1 //************************************************
2 //pollard_rho 算法进行质因数分解
3 //************************************************
4
5 LL factor[100];//质因数分解结果（刚返回时是无序的）
6 int tol;////质因数的个数。数组小标从0开始
7
8 LL gcd(LL a,LL b)
9 {
10 if(a==0) return 1;// !!!!
11 if(a<0) return gcd(-a,b);
12 while(b)
13 {
14 LL t=a%b;
15 a=b;
16 b=t;
17 }
18 return a;
19 }
20
21 LL Pollard_rho(LL x,LL c)
22 {
23 LL i=1,k=2;
24 LL x0=rand()%x;
25 LL y=x0;
26 while(1)
27 {
28 i++;
29 x0=(mult_mod(x0,x0,x)+c)%x;
30 LL d=gcd(y-x0,x);
31 if(d!=1 && d!=x) return d;
32 if(y==x0) return x;
33 if(i==k) {y=x0;k+=k;}
34 }
35 }
36
37 //对n进行素因子分解
38
39 void findfac(LL n)
40 {
41 if(Miller_Rabin(n))
42 {
43 factor[tol++]=n;
44 return;
45 }
46 LL p=n;
47 while(p>=n)
48 p=Pollard_rho(p,rand()%(n-1)+1);
49 findfac(p);
50 findfac(n/p);
51 }
