|
（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
| 