OFFSET
1,1
PROG
(Python)
import sympy
from sympy import isprime
for n in range(10000):
sum = 0
i = k = 1
while i<=n:
product = n
for x in range(k):
product *= i
i += 1
if i>n: break
sum += product
k += 1
if isprime(sum): print str(n)+', ',
CROSSREFS
KEYWORD
nonn
AUTHOR
Alex Ratushnyak, Mar 21 2014
STATUS
approved