OFFSET
1,2
COMMENTS
EXAMPLE
a(1) = 1 because 2^2 - 2 + 1 = 3 is prime.
a(2) = 2 because 6^2 - 6 + 1 = 31 is prime.
a(3) = 4 because 210^2 - 210 + 1 = 43891 is prime.
MATHEMATICA
Select[Range@ 400, PrimeQ[#^2 - # + 1 &@ Product[Prime@ k, {k, #}]] &] (* Michael De Vlieger, Nov 28 2015 *)
PROG
(PARI) a002110(n) = prod(k=1, n, prime(k));
for(n=0, 1e3, if(ispseudoprime(a002110(n)^2 - a002110(n) + 1), print1(n, ", ")))
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Altug Alkan, Nov 26 2015
STATUS
approved