(PARI)GT_Trunc1(nmax,prop,b=10)={\\有关详细信息,请参阅链接
my(n=0,v=矢量(nmax),g=1,lgs=1,Ige,an,c);
对于(k=1,b-1,如果(prop(k),v[n++]=k));
lge=n;c=lge-lgs+1;
而(c,g++;对于(k=lgs,lge,对于(m=0,b-1,an=b*v[k]+m;
如果(prop(an),v[n++]=an;如果(n>=nmax,返回(v));););
lgs=lge+1;lge=n;c=lge-lgs+1;);
如果(n,返回(v[1..n]);
打印(“无解决方案”);}
v=GT_Trunc1(1000000,isprime,16)
(PARI)isok(n)={while(n,if(!isprime(n),return(0));n\=16);1}\\乔格·阿恩特2015年3月7日
(PARI)my(A=素数([0,15]),i=1);直到(#A<i+=1,A=concat(A,素数([A[i],A[i]+1]*16))#A237600型=A\\M.F.哈斯勒2018年11月7日
(Python)
从gmpy2导入is_prime
对于范围(1,10**9)中的n:
如果is_prime(n):
s=格式(n,'x')
对于范围(1,len(s))中的i:
如果不是is_prime(int(s[:-i],16)):
打破
其他:
(Python)
从sympy导入primerange
|