ll quick_pow(ll a,ll b,ll r)//快速幂 { ll ans = 1; while(b) { if(b&1)ans=(ans*a)%r; a = (a*a)%r; b>>=1; } return ans; } bool Miller_Rabbin(int n,int a)//米勒拉宾素数测试 { int r=0,s=n-1,j; if(!(n%a)) return false; while(!(s&1)){ s>>=1; r++; } ll k=qpow(a,s,n); if(k==1) return true; for(j=0;j<r;j++,k=k*k%n) if(k==n-1) return true; return false; } bool IsPrime(int n)//判断是否是素数 { int tab[]={2,3,5,7}; for(int i=0;i<4;i++) { if(n==tab[i]) return true; if(!Miller_Rabbin(n,tab[i])) return false; } return true; }
米勒拉宾素数测试
猜你喜欢
转载自blog.csdn.net/qq_40093905/article/details/80264611
今日推荐
周排行