OFFSET
1,1
COMMENTS
LINKS
EXAMPLE
10778947368421 is in the sequence since 10778947368421 = 20^10 + 20^9 + 20^8 + 20^7 + 20^6 + 20^5 + 20^4 + 20^3 + 20^2 + 20 + 1, 20 is not prime, and 10778947368421 is prime.
MAPLE
f:= proc(n)
local p, j;
if isprime(n) then return NULL fi;
p:= add(n^j, j=0..10);
if isprime(p) then p else NULL fi
end proc:
map(f, [$1..1000]); # Robert Israel, Nov 19 2014
PROG
(Python)
from sympy import isprime
A198244_list, m = [], [3628800, -15966720, 28828800, -27442800, 14707440, -4379760, 665808, -42240, 682, 0, 1]
for n in range(1, 10**4):
for i in range(10):
m[i+1]+= m[i]
if not isprime(n) and isprime(m[-1]):
A198244_list.append(m[-1]) # Chai Wah Wu, Nov 09 2014
(Magma) [a: n in [0..500] | not IsPrime(n) and IsPrime(a) where a is (n^10+n^9+n^8+n^7+n^6+n^5+n^4+n^3+n^2+n+1)]; // Vincenzo Librandi, Nov 09 2014
(PARI) forcomposite(n=0, 10^3, my(t=sum(k=0, 10, n^k)); if(isprime(t), print1(t, ", "))); \\ Joerg Arndt, Nov 10 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Jonathan Vos Post, Dec 21 2012
EXTENSIONS
a(5)-a(6) from Robert G. Wilson v, Dec 21 2012
a(7) from Michael B. Porter, Dec 27 2012
Corrected terms a(6)-a(7) and added terms by Chai Wah Wu, Nov 09 2014
STATUS
approved