GF(b,n)=
{
my(s=Mod(2*polchebyshev(b/2,1,polchebyshev(b/2,1,4)),b^2^n+1));
for(i=1,2^n-2, s=2*polchebyshev(b,1,s/2));
s==0
}
R0YoYixuKT0KeyAKbXkocz1Nb2QoMipwb2xjaGVieXNoZXYoYi8yLDEscG9sY2hlYnlzaGV2KGIvMiwxLDQpKSxiXjJebisxKSk7CmZvcihpPTEsMl5uLTIsIHM9Mipwb2xjaGVieXNoZXYoYiwxLHMvMikpOwpzPT0wIAp9