%I#9 2019年2月21日23:14:55

%S 101111111311111111911131111611117111311114111611122229，

%电话：233333333333033334333533333833344449477774999959999，

%电话：67777717777777777377777777477777773777779777777797977779998883949997777

%至少有一个数字在十进制展开式中出现四次的N素数。

%C序列可能是无限的。

%H Harvey P.Dale，n的表格，n=1..2000的a（n）</a>

%e 101117正好有四个1。

%p isdgctm:=proc（n，d）局部dgs，a，i；dgs:=转换（n，基数，10）；a:=[seq（0，j=0..9）]；

%dgs中i的p为a:=下层（i+1=op（i+1，a）+1，a；od：如果convert（a，set）与｛d｝＜＞｛｝相交，则为true；否则为假；fi；结束：

%对于从1到10000的n，p表示p:=ithprime（n）；如果是isdgctm（p，4），则打印f（“%d，”，p）；fi；od:#R.J.Mathar_，2009年6月21日

%t选择[Prime[Range[18610000]]、MemberQ[DigitCount[#]、4]和]（*哈维·P·戴尔，2015年4月12日*）

%K nonn，基础

%O 1，1

%A _ Ki Punches，2009年6月19日

%E由R.J.Mathar_编辑和更正，2009年6月21日