%I#34 2019年11月11日03:51:46

%编号：1824147616322462248425203064340238204026664070267158，

%电话：907012202125481299413042153581764617670183361956420624，

%电话：225002412626132261882624029074296583077831126324304434016

%N数N，使N^1024+1是素数（广义费马素数）。

%C根据Bunyakovsky的猜想，这个序列是无限的_Charles R Greathouse IV，2012年4月26日

%H T.D.Noe，n表，n=1..1000的a（n）（摘自Yves Gallot）

%H Yves Gallot，<a href=“http://yves.gallot.pagesperso-orange.fr/primes网站/“>广义费马素数搜索</a>

%H Jeppe Stig Nielsen，<a href=“http://jeppesn.dk/generalized-fermat.html“>按基数排序的广义费马素数。

%H Eric Weisstein的数学世界，<a href=“http://mathworld.wolfram.com/GeneralizedFermatNumber.html“>广义费马数</a>

%t执行[k=1；而[PowerMod[n，1024，2*k*1024+1]！=2*k*1024&&k<2*10^6，k++]；如果[k==2*10^6&&PrimeQ[n^1024+1]，打印[n]]，{n，2，13954，2}]

%t Do[If[PrimeQ[n^1024+1]，Print[n]，##&[]]，{n，1，100}]（*包括第一学期，跑得更快，_Daniel Jolly_，2014年11月4日*）

%o（PARI）is A057002（n）=isprime（n^1024+1）\\_Michael B.Porter_，2010年4月3日

%Y其他数字序列n，其中n^（2^k）+1是固定k的素数：A005574、A000068、A006314、A00631.3、A0063.15、A006316、A056994、A0569965、A057465、A088361、A0883652、A226528、A22652、A226530、A251597、A253854、A244150、A243959、A321323。

%Y参考A006093。

%K nonn公司

%O 1,2号机组

%A _Robert G.Wilson v_，2000年9月9日

%E更多来自Jeppe Stig Nielsen的条款，2003年9月27日

%E根据T.D.Noe的建议由N.J.A.Sloane编辑，2008年5月14日