SHARE
TWEET

Untitled

a guest Aug 13th, 2017 47 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. CE(k1,k2,b,n1,n2)=
  2. {
  3. forstep(k=k1,k2,[1,2],
  4. for(n=n1,n2,
  5. N=k*b^n-1;
  6. my(s=Mod(2*polchebyshev(k*b/2,1,polchebyshev(b/2,1,2)),N));
  7. for(i=1,n-2,s=2*polchebyshev(b,1,s/2));
  8. if((s==0) && !ispseudoprime(N),print("k=",k," b=",b," n=",n))))
  9. }
RAW Paste Data
Top