meditationatae

Just another WordPress.com site

Some Rabin-Miller pseudoprime primality tests

? m=55;s=m;for(X=1,6,s=s^4-s^2+1);print(m,” “,round(log(s)/log(10)))
55 7128

? ispseudoprime(s,30)
%39 = 1

? m=1948;s=m;for(X=1,6,s=s^4-s^2+1);print(m,” “,round(log(s)/log(10)))
1948 13474

? ispseudoprime(s,30)
%40 = 1

? m=3269;s=m;for(X=1,6,s=s^4-s^2+1);print(m,” “,round(log(s)/log(10)))
3269 14395

? ispseudoprime(s,30)

 

Advertisements

Written by meditationatae

May 30, 2013 at 9:36 pm

Posted in History

Tagged with , ,

%d bloggers like this: