(PARI)isokp(n)={v=数字(n);对于(k=1,#v,w=v;如果(k==1,idep=1,idp=0);对于
列表(nn)={forprime(p=2,nn,if(isokp(p),print1(p,“,”));}\\米歇尔·马库斯2015年12月15日
(Python)
从sympy导入isprime
def h1(n):#hamming distance 1是n的邻居,不以0开头
s=str(n);d=“0123456789”;L=长度
如果c=s[i]和not(i==0和c==“0”)
def ok(n):返回isprime(n)和all(对于h1(n)中的k,不是isprim(k))
打印([k代表范围内的k(10**6),如果正常(k)])#迈克尔·布拉尼基2022年7月31日