5777同时满足变参数的的lucas伪素数与76为底的强伪素数

mathematica

http://ntheory.org/data/slpsps.txt
上面的是lucas伪素数，变参数的

1. MillerRabin[n0_,a0_]:=Module[{n=n0,a=a0,s,m,t1,k},
   
2.     s=0;m=n-1;While[Mod[m,2]==0,m=m/2;s=s+1];
   
3.     t1=PowerMod[a,m,n];
   
4.     If[t1==1,Return[True]];
   
5.     k=0;While[k<s-1&&t1!=n-1,k=k+1;t1=Mod[t1^2,n]];
   
6.     If[t1==n-1,Return[True],Return[False]]
   
7. ]
   
8. 
   

复制代码

mathematica

5701
是素数，且5777通过了这个SPRP测试，
为什么呢？
因为5701-5777=-76