Pollard_rho算法

思路代码
int find_factorplus(int N) {
    a = 2;
    for（ int i= 1; i <= 1000000; i++ ) {
        b = f(a);
        p = GCD( abs( b - a ) , N);
        if( p > 1 ) return p;//Found factor: p
        a = b;
    }
    return 0;//Failed. :-(
}

floyd判圆代码
int find_factorplus(int N) {
    a = 2;
    b = a;
    do {
        a = f(a);//a runs once
        b = f(f(b));//b runs twice as fast
        p = GCD( abs( b - a ) , N);
        if( p > 1 ) return p;//Found factor: p
    } while( b != a );
    return 0;//Failed. :-(
}

1.需要一个随机数函数f（x）=(x^2+1)%N
2一般设置一个初始变量a=2,另外需要一个随机变量b用于判环，其值为f(f(x))
3.当a==b的时候就是出现循环，这时退出，当p>1时就是找到一个因子，可以继续递归循环找另一个因子。