%I#46 2022年2月21日01:12:37

%S 1,7,13,17,19,31,37,43,47,49,59,61,67,71,73,79,91,97101103107109，

%电话127133137139149151157163167169181193199211217223，

%电话：227229241247259263269271278283289

%N伯努利数B_{2n}的分母为6。

%C A053176中的所有素数都在序列中。如果n在序列中，它的因式分解只包含A053176中的素数_Benoit Cloitre_，2002年10月19日

%C B（2n）有分母6，当（n^2-1）*B（2n）是整数时_Benoit Cloitre_，2004年2月15日

%C A156543的后续。-_Reinhard Zumkeller，2009年2月10日

%D B.C.Berndt，Ramanujan笔记本第四部分，Springer-Verlag，见第76页。

%H T.D.Noe，n的表格，n=1..1000的a（n）</a>

%H<a href=“/index/Be#Bernoulli”>为与伯努利数相关的序列的条目建立索引</a>

%F a（n）似乎对c*n渐近，5<c<6。-_Benoit Cloitre_，2002年10月19日

%t ok[n_]：=整数Q[（n^2-1）*BernoulliB[2n]]；选择[范围[300]，确定]（*_Jean-François Alcover_，2012年6月27日，在_Benoit Cloitre_*之后）

%t结果={}；做[count=0；

%t Do[If[Not[PrimeQ[2*Divisors[n][[i]]+1]]，计数++]，

%t{i，2，除数Sigma[0，n]}]；如果[count==DivisorSigma[0，n]-1，AppendTo[result，n]，{n，1，10000}]；结果（*Reichard R.Forberg，2016年8月6日*）

%t位置[BernoulliB[2范围[300]]，_？（分母[#]==6&）//压扁（*哈维·P·戴尔，2017年1月28日*）

%o（PARI）isok（n）=分母（bernfrac（2*n））==6；\\_米歇尔·马库斯，2016年2月6日

%o（Magma）[n:n in[0..400]|分母（伯努利（2*n））eq 6]；//_文森佐·利班迪（Vincenzo Librandi），2016年2月6日

%Y参考A051222、A053176、A156543。

%K nonn，很好，很容易

%O 1,2号机组

%A _N.J.A.斯隆_