meditationatae

Just another WordPress.com site

Some Rabin-Miller pseudoprime primality tests (updated)

? 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)
%41 = 1

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

? ispseudoprime(s,30)
%42 = 1

Advertisements

Written by meditationatae

May 30, 2013 at 10:18 pm

Posted in History

Tagged with

%d bloggers like this: