%I#23 2022年9月8日08:46:00

%S 5,11111511051115115115011501015501101111015110501，

%电话：1155115101155111555150051501015011150551510015115151511，

%电话：5155155001550155501555111001511005111011011115011100511110501115001115101115001151500011150011550151051

%N只包含数字“0”、“1”和“5”的素数。

%H Robert Israel，n的表，n的a（n）=1..10000</a>

%H<a href=“/index/Pri#PrimesWithDigits”>为给定集合中具有数字的素数项建立索引</a>

%p N:=10000:#获得前N个项

%p计数：=0：

%允许p：={0,1,5}：

%允许：=nops（允许）：

%p subst:=seq（i=允许[i+1]，i=0..允许-1）：

%当计数<N do时，从0开始p代表d

%当计数<N do时，x1的p从1降到1

%当计数<N do时，t的p从0到允许的^d-1

%p L：=subs（subst，convert（x1*nallowed^d+t，base，nallowed））；

%pX:=加（L[i]*10^（i-1），i=1..d+1）；

%p如果是素数（X），则

%p计数：=计数+1；

%p A[计数]：=X；

%功率因数

%日期日期：

%p序列（A[n]，n=1..n）；#_罗伯特·伊斯雷尔，2014年4月20日

%t选择[FromDigits/@Tuples[{0,1,5}，6]，PrimeQ]（*Harvey P.Dale_，2021年7月23日*）

%o（PARI）L=[0,1,5]；对于（d=1,6，u=向量（d，i，10^（d-i））~；对于vec（v=向量（d，i，[1+（i==1&！L[1]），#L]），ispseudoprime（t=向量（d，i，L[v[i]]）*u）&print1（t“，”））/*函数a（n）请参见A199327*/

%o（岩浆）[0，1，5]]；//_文森佐·利班迪（Vincenzo Librandi），2014年4月22日

%Y参见A020449-A020472、A199326-A199329、A061247、A199340-A199349。

%K nonn，基础

%O 1，1

%A _M.F.Hasler，2011年11月5日