标签：
以前我在判断素数上一直只会 sqrt(n) 复杂度的方法和所谓的试除法(用sqrt(n)以内的素数来除)。
现在我发现其实还有一种方法叫做费马小定理。
有关请见 http://baike.baidu.com/link?url=1BurQrmJP3j9QiD4OnA2X3TAbSSCPvTgbaqbo6qSQPVSuXLjVe-lL2SNi6N5wblwJFrIJs41pmDbCZ6z9je4h_；
代码如下：
llg ch(llg a,llg b){ llg now=a,ans; while(b){ if(b&1) ans+=now,ans%=mod; b>>=1; now<<=1; now%=mod; } return ans; } bool mi(llg a,llg b){ llg s=1; while(b){ if(b & 1) s=ch(s,a); a=ch(a,a); b>>=1; } return s==1; } bool pd(llg x){ if(x==2) return 1; mod=x; for(int i=1,y;i<=5;i++){ y=rand()+1; if(!mi(y,x-1)) return 0; } return 1; }
其中ch这个函数是快速乘法，可以在10^18大小的模数下相乘。
费马小定理的逆定理不一定成立，于是我们就需要多取几组加以实验。只要你人品不是太坏，就不会错。考验RP的时刻到来了
大概就是这样了吧。
标签：
原文地址：http://www.cnblogs.com/lcf-2000/p/5547695.html
