%I#23 2012年7月23日19:32:39

%S 1,2,3,4,5,6,7,8,9,13,14,15,19,21,23,26,40,43,45,52,54,55,69,77,90,99，

%电话：1061281471762022673314585125559089421004112313741386，

%U 1467型

%N个数b，这样至少有一个数c和一个单位数d，这样（10^c-d）*10^b-1和（10^c-d）*10 ^b+1是0<c<2*b的双素数。

%C在这种格式中，一位数d总是1、4或7，因为（10^C-d）*10^b+-1中的一个是3的倍数。

%C对于一个b，可能有多个匹配项（C，d）。

%C相关最小C值的顺序为：1、1、1，2、3、6、1、4、2、11、9、4、7、12、9、9、42、62、5、31、2、72、88、141、119、181、6、38、164、132、53、293、150、704、557、980、952、1596、529、2221、200、169、1371，。。。其相关d值为4、4、1、1、1,1、7、4,7、4,1、1、7,1、4、1,7、1,7,7、4、1,1、7,4,4、4,7,1、1,7、7,4、1。。。。

%e（10^1-7）*10^1-1=29素数31双素数，因此a（1）=1。

%e（10^1-4）*10^2-1=599素数601双素数，因此a（2）=2。

%e（10^1-1）*10^3-1=8999素数9001双素数，因此a（3）=3。

%e（10^2-1）*10^4-1=989999素数990001双素数，因此a（4）=4。

%e（10^3-1）*10^5-1=99899999素数。

%e（10^3-1）*10^5+1=99900001双素数，所以a（5）=5。

%p为A213882:=过程（b）

%p局部c，d，p；

%p代表从1到2的c*b-1 do

%p代表d从0到9 do

%p p：=（10^c-d）*10^b-1；

%p如果isprime（p）和isprim（p+2），则

%p返回true；

%p结束if；

%p端do：

%p端do：

%p返回false；

%p端程序：

%p代表n从1到2000 do

%p如果是A213882（n），则

%p printf（“%d，\n”，n）；

%p end if；

%p结束do；#_R.J.Mathar，2012年7月21日

%Y参见A213883、A213884。

%K nonn公司

%O 1,2号机组

%A Pierre CAMI，2012年6月26日