Miller Robin大素数判定

Miller Robin算法 当要判断的数过大,以至于根n的算法不可行时,可以采用这种方法来判定素数。 用于判断大于2的奇数(2和偶数需要手动判断),是概率意义上的判定,因此需要做多次来减少出错概率。 Template: ```C++ typedef long long ll; ll kmul(ll a,ll b,ll mod) { ll res=0; while (b) { if (b&1) res=(res+a)%mod; a=(a+a)%mod; b>>=1; } return res; } ll kpow(ll a,ll b,ll mod) { ll res=1; while (b) { if (b&1) res=kmul(res,a,mod)%mod; a=kmul(a,a,mod)%mod; b>>=1; } return res; } bool Mil_Rb(ll n,ll a) { ll d=n-1,s=0,i; while (!(d&1)) { d>>=1; s++; } ll t=kpow(a,d,n); if (t==1||t==-1) return 1; for (i=0;ia[i]&&!Mil_Rb(n,a[i])) return 0; } return 1; } ```
posted @ 2018-07-25 00:48  __orange  阅读(2057)  评论(0编辑  收藏  举报