%I#27 2019年1月17日13:44:05

%序号：0,1,3,5,931593591469289730933111156971795542261111031

%N序列中素数的指数由A（0）=11定义，当N>0时，A（N）=10*A（N-1）+21。

%C数字n，使得（120*10^n-21）/9是素数。

%C数字n，使数字1后跟n>=0，数字3后跟数字1是质数。

%C与<=3111项相对应的数字是认证素数。有关较大的数字，请参阅P.De Geest，PDP参考表。

%D Klaus Brockhaus和Walter Oberschelp，Zahlenfolgen mit hominem Ziffernkern，MNU 59/8（2006），第462-467页。

%H Patrick De Geest，<a href=“http://www.worldofnumbers.com/deplat.htm#pdp131“>PDP参考表-131。

%H Makoto Kamada，<a href=“https://stdkmd.net/nrr/1/13331.htm#prime“>133…331形式的质数。

%H<a href=“/index/Pri#Pri_rep”>包含重单位的素数的索引条目。

%如果n>1，F a（n）=A082697（n-1）-2。

%e 131是质数，因此1是项。

%t做[If[PrimeQ[（1*10^n+3*（10^n-1）/9）*10+1]，打印[n]]，{n，1，2500}]

%t选择[Range[0，2000]，PrimeQ[（120 10^#-21）/9]&]（*文森佐·Librandi_，2014年11月3日*）

%o（PARI）a=11；对于（n=01500，如果（i素数（a），打印1（n，“，”））；a=10*a+21）

%o（PARI）表示（n=01500，如果是（i素数（（120*10^n-21）/9），打印1（n，“，”））

%Y参考A000533、A002275、A068645、A082697。

%K nonn，难，更多

%氧1,3

%A _Robert G.Wilson v_，2000年8月18日

%E Klaus Brockhaus_和Walter Oberschelp（Oberschelp（AT）informatik.rwth-aachen.de）的更多条款和附加评论，2004年12月28日

%E编辑：N.J.A.Sloane，2007年6月15日

%Herman Jamke（hermanjamke（AT）fastmail.fm）于2008年1月1日在De Geest网站上发布的E更新

%E a（15）=111031，来自雷·钱德勒，2011年4月14日

%E更新评论部分和链接，作者：Patrick De Geest_，2014年11月2日

%E编辑：雷·钱德勒，2014年11月4日