快速分解算法

算法思路

太懒了不写
参见JZOJ4458下的pdf

代码

    while (1){
        a=rand()%1000+2;
        r1=r2=2;
        do{
            r1=f(r1);
            r2=f(f(r2));
            p=gcd(abs(r2-r1),n);
            if (p>1) break;
        }while (r1!=r2);
        if (r1!=r2) break;
    }
    printf("%d\n",p);