%I#42 2022年7月31日11:05:03

%S 2940015054475841410604171929573971767106259912825291524181，

%电话：201796324744312602013070663308555332648943931395152507，

%电话：5285767556445355752595974249617316191371623617964632676712591720477774697897977810223

%弱素数（或孤立素数）：改变任何一个十进制数字总是产生一个复合数，但第一个数字不能变为0。

%C该定义可以重述为“带d位的素数p，这样在距p的汉明距离1处（以10为基数）就不存在最多带d位数字的素数q”_N.J.A.斯隆，2019年5月6日

%C对于下列k、5、6、7、8、9、10的值，此序列中<10^k的项数为0、6、43、406、3756、37300_Jean-Marc Rebert，2015年11月10日

%H Jean-Marc Rebert，n表，n=1..3756的a（n）</a>

%H C.里维拉，<a href=“http://www.primepuzzles.net/puzzles/puzz_017.htm“>弱素数</a>

%H Eric Weisstein的数学世界，<a href=“http://mathworld.wolfram.com/WeaklyPrime.html“>弱素数</a>

%p过滤器：=程序（n）

%p局部L，i，d，ds；

%p如果不是isprime（n），则返回假fi；

%p L：=换算（n，基数，10）；

%p代表i从1到nops（L）do

%p如果i=nops（L），则ds:={$1..9}减去{L[i]}

%p elif i＝1，则ds：＝{1,3,7,9}减去{L[i]}

%p其他ds:={$0..9}减去{L[i]}

%p fi；

%ds-do中的p代表d

%p如果isprime（n+（d-L[i]）*10^（i-1）），则返回false fi；

%日期

%p od；

%p为真

%p端程序：

%p选择（过滤器，[seq（i，i=11..10^6,2）]）；#_罗伯特·伊斯雷尔，2015年12月15日

%t选择[Prime@Range[10^5]，函数[n，函数[w，Total@Map[Boole@PrimeQ@#&，DeleteCase[#，n]]&@Union@Flatten@Map[Function[d，FromDigits@ReplacePart[w，d->#]&/@If[d==1，#，Prepend[#，0]&@Range@9]，Range@Length@w]==0]@IntegerDigits@n]]（*Michael De Vlieger_，2016年12月13日*）

%o（PARI）isokp（n）={v=数字（n）；对于（k=1，#v，w=v；如果（k==1，idep=1，idp=0）；对于

%o lista（nn）={forprime（p=2，nn，if（isokp（p），print1（p，“，”））；）；}\\_Michel Marcus_，2015年12月15日

%o（Python）

%o来自sympy import isprime

%o def h1（n）：#hamming距离n的1个邻居，不是从0开始

%o s=str（n）；d=“0123456789”；L=长度

%o如果c=s[i]和not（i==0和c==“0”）

%o def ok（n）：返回isprime（n）和all（不是h1（n）中k的isprime

%o打印（[k表示范围内的k（10**6），如果ok（k）]）#_Michael S.Branicky_，2022年7月31日

%Y参见A050249、A158125（弱底漆）、A186995、A192545。

%K nonn，基础

%O 1,1号机组

%A _Eric W.Weisstein，2009年3月13日

%E由_Charles R Greathouse IV_编辑，2010年8月2日

%E a Andrew Howroyd_于2018年2月23日将缺失的a（3385）插入b文件