%I#18 2019年3月5日16:24:36

%S 491617647115112171811184729383340034657495167637687，

%电话：8831901104631058912073132171453314737149571528715787，

%电话：158231600717681178631871318869205332093924019246592515326561

%N不规则素数，不规则指数为3。

%H Amiram Eldar，n的表，n=1..100000时的a（n）（由Hart等人计算，T.D.Noe中的术语1..9824，由Buhler等人计算）

%H J.Buhler、R.Crandall、R.Ernvall、T.Metsankyla和M.A.Shokrollahi，<A href=“http://dx.doi.org/10.1006/jsco.1999.1011“>不规则素数和分圆不变量到1200万</a>，《符号计算杂志》31，2001，89-96。

%H贝努利数，素数的不规则指数</a>

%H William Hart、David Harvey和Wilson Ong，<a href=“https://doi.org/10.1090/mcom/3211“>不规则素数到20亿</a>，《计算数学》，第86卷，第308期（2017年），第3031-3049页；也可在<a href=”https://arxiv.org/abs/1605.02398“>arXiv:1605.02398[math.NT]</a>，2016年。

%H David Harvey，<a href=“https://web.maths.unsw.edu.au/~davidharvey/papers/twobillion/“>不规则素数到20亿（包括所有小于2^31的素数的列表）。

%t Do[p=素数[n]；k=1；c=0；当[2*k<p-3时，如果[Mod[Numerator[BernoulliB[2*k]]，p]==0，c++]；k++]；如果[c==3，打印[p]]，{n，3，1000}]

%t做[p=素数@n；k=1；c=0；当[2*k<p-3时，如果[Mod[Numerator[BernoulliB[2*k]]，p]==0，c++]；k++]；如果[c==3，打印@p]，{n，313887}]

%Y参见A000928、A000367、A060974、A060995、A073276和A073277。

%K nonn公司

%O 1,1号机组

%A _Robert G.Wilson v_，2002年7月22日

%E由_Robert G.Wilson v_扩展，2006年9月20日