另一个Pollard Rho实施

时间:2015-07-30 21:18:07

标签: algorithm primes factorization

为了解决项目Euler(https://projecteuler.net/problem=3)上的第3个问题,我决定实施Pollard的Rho算法(至少部分算法,我计划包括骑行后来)。奇怪的是它适用于如下数字:82123(因子= 41)和16843009(因子257)。然而,当我尝试项目欧拉数:600851475143时,当最大素因子为6857时,我最终获得71.这是我的实现(对不起代码墙和缺少类型转换):

#include <iostream>
#include <math.h>
#include <vector>

using namespace std;

long long int gcd(long long int a,long long int b);
long long int f(long long int x);

int main()
{
long long int i, x, y, N, factor, iterations = 0, counter = 0;
vector<long long int>factors;

factor = 1;
x = 631;
N = 600851475143;
factors.push_back(x);

while (factor == 1)
{
    y = f(x);
    y = y % N;
    factors.push_back(y);
    cout << "\niteration" << iterations << ":\t";
    i = 0;
    while (factor == 1 && (i < factors.size() - 1))
    {
        factor = gcd(abs(factors.back() - factors[i]), N);
        cout << factor << " ";
        i++;
    }
    x = y;
    //factor = 2;
    iterations++;
}

system("PAUSE");
return 0;
}

long long int gcd(long long int a, long long int b)
{
    long long int remainder;
    do
    {
    remainder = a % b;
    a = b;
    b = remainder;
    } while (remainder != 0);
    return a;
}

long long int f(long long int x)
{
//x = x*x * 1024 + 32767;
x = x*x + 1;
return x;
}

1 个答案:

答案 0 :(得分:2)

Pollard的rho算法不保证任何东西。它不能保证找到最大的因素。它并不保证它找到的任何因素都是素数。它甚至不能保证找到一个因素。 rho算法是概率性的;它可能会找到一个因素,但不一定。由于你的函数返回一个因子,它可以工作。

那就是说,你的实施并不是很好。没有必要存储函数的所有先前值,并且每次通过循环计算gcd。这是伪代码,用于更好的函数版本:

function rho(n)
    for c from 1 to infinity
        h, t := 1, 1
        repeat
            h := (h*h+c) % n # the hare runs ...
            h := (h*h+c) % n # ... twice as fast
            t := (t*t+c) % n # as the tortoise
            g := gcd(t-h, n)
        while g == 1
        if g < n then return g

此函数返回单个因子n,可以是素数也可以是复合因子。它只存储随机序列的两个值,并在找到一个循环时停止(当g == n时),用不同的随机序列重新启动(通过递增c)。否则它会继续运行,直到找到一个因子,只要将输入限制为64位整数,就不会花费太长时间。通过将rho应用于剩余的辅因子来查找更多因子,或者如果找到的因子是复合因子,则在找到所有素因子时停止。

顺便说一句,你不需要Pollard的rho算法来解决Project Euler#3;简单的试验分工就足够了。该算法查找数字的所有素因子,从中可以提取最大值:

function factors(n)
    f := 2
    while f * f <= n
        while n % f == 0
            print f
            n := n / f
        f := f + 1
    if n > 1 then print n