|
（PARI）eva（n）=subst（Pol（n），x，10）
insert_zeros（num，len）=我的（d=数字（num），v=[]）；对于（k=1，#d-1，v=concat（v，concat，[d[k]]，vector（len）））；v=连接（v，d[#d]）；eva（v）
a（n）=素数（p=10，对于（k=1，n，如果（！ispseudoprime（eva（insert_zeros（p，k））），中断，如果（k==n，return（p））））
（Python）
从sympy导入isprime，nextprime
def insert_zeros（n，k）：返回int（（“0”*k）.join（list（str（n）））
def-ok（p，n）：返回范围（1，n+1）中k的全部（isprime（insert_zeros（p，k）））
定义a（n，startat=11）：
p=启动
为True时：
如果正常（p，n）：返回p
p=下一素数（p）
打印（[a（n）代表范围（1，6）中的n]）#迈克尔·布拉尼基2021年6月3日
| 