%I #21 Jul 13 2022 14:08:00
%S 2,3,5,13,17,31,61,73,89,149,163,251,349,353,461,523,599,647,863,941,
%T 947,1087,1117,1229,1277,1291,1297,1409,1439,1489,1567,1579,1609,1627,
%U 1753,1783,1831,2039,2131,2293,2531,2609,2753,2861,3037,3163,3167,3299
%N Primes p such that 1 + p + p^2 + p^3 + p^4 + p^5 + p^6 is prime.
%C Primes in A100330. The generated prime numbers are exactly A194257. [_Bernard Schott_, Dec 21 2012]
%H Robert Israel, <a href="/A163268/b163268.txt">Table of n, a(n) for n = 1..10000</a> (first 102 terms from Zak Seidov)
%p select(p -> isprime(p) and isprime(1+p+p^2+p^3+p^4+p^5+p^6), [2,seq(i,i=3..10000,2)]); # _Robert Israel_, May 05 2017
%t f[n_]:=1+n+n^2+n^3+n^4+n^5+n^6; lst={}; Do[p=Prime[n]; If[PrimeQ[f[p]], AppendTo[lst,p]], {n,7!}]; lst
%t Select[Prime[Range[500]],PrimeQ[Total[#^Range[0,6]]]&] (* _Harvey P. Dale_, Jul 13 2022 *)
%o (PARI) n=0;forprime(p=2,10000,isprime((p^7-1)/(p-1))&&print(n++" "p))\\ _Zak Seidov_, Mar 09 2013
%Y Cf. A053182, A065509, A088550, A100330, A163268, A194257.
%K nonn
%O 1,1
%A _Vladimir Joseph Stephan Orlovsky_, Jul 24 2009
%E Edited (but not checked) by _N. J. A. Sloane_, Jul 25 2009