(PARI)rot(n)=如果(#Str(n)==1,v=向量(1),v=矢量(#n-1));对于(i=2,#n,v[i-1]=n[i]);u=矢量(#n);对于(i=1,#n,u[i]=n[i]);v=连接(v,u[1]);v(v)
eva(n)=x=0;对于(k=1,#n,x=x+(n[k]*10^(#n-k)));x个
a(n)=i=0;r=rot(数字(n));当(r!=数字(n)时,如果(ispseudoprime(eva(r)),i++);r=rot(r));如果(ispseudoprime(eva(r)),i++);我
|