%I#27 2022年3月23日14:40:20

%S 1,2,3,5,6,8,9,11,15,21,30,39,50,63,83,95,991561733508541308，

%电话：17692903525053455639619572392136841669476845861963515，

%电话：6946870539133508134993187160493095

%N数N，使得1+Sum_{i=1..N}2^（2i-1）是素数。

%如果这个序列是无限的，那么A124401也是无限的。

%C等于A127965（n）/2。

%C和具有简单闭式1+2/3*（4^n-1）_Stefan Steinerberger_，2007年11月24日

%C超过a（30）的项对应于可能素数，参见A000978_M.F.Hasler，2008年8月29日

%F a（n）=地板（A000978（n）/2）=天花板（对数（4）（A000799（n）））；A000978（n）=2a（n）+1；A000979（n）=（2*4^a（n）+1）/3.-_M.F.Hasler，2008年8月29日

%e a（1）=1，因为1+2=3是质数；

%e a（2）=2，因为1+2+2^3=11是素数；

%e a（3）=3，因为1+2+2^3+2^5=43是质数；

%e a（4）=5，因为1+2+2^3+2^5+2^7+2^9=683是质数；

%e。。。

%t a＝｛｝；Do[If[PrimeQ[1+Sum[2^（2n-1），{n，1，x}]]，AppendTo[a，x]]，{x，1，1000}]；一

%tb={}；做[c=1+和[2^（2n-1），{n，1，x}]；如果[PrimeQ[c]，AppendTo[b，c]]，{x，0，1000}]；a={}；Do[AppendTo[a，FromDigits[IntegerDigits[b[x]]，2]]，{x，1，Length[b]}]；d={}；Do[AppendTo[d，（1/2）（数字计数[a[x]]，10，0]+DigitCount[a[x]]，10,1]]），{x，1，长度[a]}]；d日

%t位置[累计[2^（2*范围[1000]-1）]，_？（PrimeQ[#+1]&）]//Flatten（*程序生成序列的前21个项。要生成更多项，请增加Range常量。*）（*哈维·P·戴尔，2022年3月23日*）

%o（PARI）for（n=1999，ispseudoprime（2^（2*n+1）\3+1）&print1（n“，”））\\_M.F.Hasler_，2008年8月29日

%o（哈斯克尔）

%o导入数据。列表（findIndices）

%o a127936 n=a127936_列表！！（n-1）

%o a127936_list=查找索引（（==1）。a010051“）a007583_列表

%o---Reinhard Zumkeller，2013年3月24日

%o（Python）

%o来自sympy import isprime

%o A127936=[i表示i在范围（1,10**3）中，如果是i素数（int（'01'*i+'1'，2））]

%o#_Chai Wah Wu_，2014年9月5日

%Y请参阅A127962、A127963、A1279964、A1279605、A12796、A000979、A00097、A124400、A126614、A179955、A12.7956、A127957、A127988、A127936、A179936、A1 24401、A010051、A007583。

%K nonn，更多

%O 1,2号机组

%A _Artur Jasinski，2007年2月8日，2007年02月09日

%E根据Andrew S.Plewe的建议，由N.J.A.Sloane编辑，2007年6月11日

%2007年11月24日，斯泰芬·施泰因贝谢（_Stefan Steinerberger_）的E 2条新条款

%E 2008年7月12日来自_Dmitry Kamenetsky的其他6个术语

%E a（30）-a（40）由M.F.Hasler_于2008年8月29日根据A000978计算得出