%I#37 2019年1月23日03:06:23

%S 1,1,1,2,2,3,7,8,13,23,41,67111193360630109119835586448，

%电话：1187621649401517365813571125178646867887524716340693060794，

%电话：574624510806204203569213843339872656139137562095

%N Sophie-Germain素数介于2^N和2^（N+1）之间。

%精确地说，Sophie-Germain素数p，2^n<p<=2^（n+1）。因为2是一个Sophie-Germain素数，所以这个精确的定义只对确定a（0）和a（1）很重要。替代定义（2^n<=p<2^（n+1））将给出序列0，2，1，1，2，2，3，7，8，13，23，41，67，111，193。。。

%C Sophie Germain素数p在A005384中。相应的素数s=2p+1称为安全素数，在A005385中。2^（n+1）到2^。

%H Paul D.Beale，<a href=“http://arxiv.org/abs/1411.2484“>基于Pohlig-Hellman指数密码的一类新的可扩展并行伪随机数生成器，arXiv预印本arXiv:1411.24842014-2015。

%H Paul D.Beale，Jetanat Datephanyawat，<a href=“https://arxiv.org/abs/1811.11629“>基于非加密RSA指数密码的可扩展并行和矢量化伪随机数生成器类，arXiv:1811.11629[cs.CR]，2018。

%F a（n）=A211397（n+1）-A211397（n）.-_米歇尔·马库斯，2014年9月22日

%t nmax=36；rtable=表[0，{nmax}]；

%t做[r=0；

%t做[If[And[PrimeQ[i]，PrimeQ[2i+1]]，r++]，{i，1+2^n，

%t2^（n+1）}]；打印[n，“”，r]；

%可传输[[n+1]=r，{n，0，nmax-1}]；

%可移植（*Paul D.Beale_，2014年9月19日*）

%o（PARI）a211395（n）={局部（r，i）；r=0；对于（i=2^n+1，2^（n+1），如果（isprime（i）&&isprim（2*i+1），r=r+1））；r}\\_Michael B.Porter_，2013年2月8日

%Y参见A005384、A005385、A211397。

%K nonn公司

%0、5

%A _Brad Clardy，2013年2月8日

%E a（29）-a（36）摘自Paul D.Beale，2014年9月19日