%I#44 2022年9月8日08:45:15

%S 1,2,3,5,6,13,14,17,26,31,38,40,46,56,60,61,66,68,72,73,80,87,89,93，

%电话：95115122126128146149156158160163180186192203206208，

%电话：220221235237238251226426628028222902943003330349350

%N个正整数k，使得k^6+k^5+k^4+k^3+k^2+k+1是素数。

%C相应的素数为A088550_Bernard Schott，2012年12月20日

%Ck=5978493*2^150006-1是这个序列中一个非常大的项的例子。生成的素数由N-1方法证明（因为k是素数，k*（k+1）是完全因子化的，这为Phi_7（k）-1提供了33.33…%的因子化）_谢尔盖·巴塔洛夫，2015年3月13日

%H Vincenzo Librandi，n的表，n=1..1000的a（n）</a>

%e2在序列中是因为2^6+2^5+2^4+2^3+2^2+2+127是素数。

%p A100330：=程序（n）

%p选项记忆；

%p局部a；

%p如果n=1，则

%第1页；

%p其他

%p表示来自procname（n-1）+1 do的a

%p如果是素数（numtheory[分圆]（7，a）），则

%p返回a；

%p end if；

%p端do：

%p end if；

%p端程序：

%p序列（A100330（n），n=1..30）；#_R.J.Mathar，2014年2月7日

%t选择[Range[350]，PrimeQ[Sum[#^i，{i，0，6}]&]（*雷·钱德勒，2004年11月17日*）

%t Do[如果[PrimeQ[n^6+n^5+n^4+n^3+n^2+n+1]，打印[n]]，{n，1，500}]（*_文森佐图书馆_，2014年2月8日*）

%o（岩浆）[n:n在[1..500]|IsPrime（n^6+n^5+n^4+n^3+n^2+n+1）]；//_文森佐·利班迪（Vincenzo Librandi），2014年2月8日

%o（PARI）is（n）=isprime（polcyclo（7，n））\\_Charles R Greathouse IV_，2015年4月28日

%Y参考A100331，A250174（Phi_14（n）=n^6-n^5+n^4-n^3+n^2-n+1素数；这两个序列也可以被认为是彼此扩展到负n值）、A250177（Phi21（n）素数）。

%K nonn公司

%O 1,2号机组

%2004年11月16日，A_Ray G.Opao