（岩浆）IsA118118:=函数（n）；D： =整数（n）；返回[1..#D]中的所有{k，j>：k，[0..9]|j eq D[k]或不是IsPrime（Seqint（S）），其中S:=插入（D，k，k，[j]）}；端函数；[p:p在PrimesUpTo（8000000）|IsA118118（p）中]//克劳斯·布罗克斯2011年2月28日
（PARI）isokp（n）={v=数字（n）；for（k=1，#v，w=v；for（j=0，9，if（j！=v[k]，w[k]=j；ntest=subst（Pol（w），x，10）；if（isprime（ntest），return（0）；）；）
列表（nn）={forprime（p=2，nn，if（isokp（p），print1（p，“，”））；}\\米歇尔·马库斯2015年12月15日
（Python）
从sympy导入isprime
def h1（n）：#n的hamming距离1邻居
s=str（n）；d=“0123456789”；L=长度
如果c=第[i]节）
def ok（n）：如果k=n）
打印（[k代表范围内的k（10**6），如果正常（k）]）#迈克尔·布拉尼基2022年6月19日
