pollard_rho 算法进行质因数分解

 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 }