OFFSET
1,2
COMMENTS
All terms are odd. - Jon Perry, Sep 11 2014
LINKS
Nathaniel Johnston, Table of n, a(n) for n = 1..10000
EXAMPLE
a(5)=31 since 31^3+(31+1)^3+(31+3)^3=101863 is prime.
MAPLE
isA246666 := proc (n) return isprime(n^3+(n+1)^3+(n+3)^3) end proc; seq(`if`(isA246666(2*n-1), 2*n-1, NULL), n = 1 .. 400); # Nathaniel Johnston, Sep 09 2014
PROG
(PARI) for(n=0, 10^3, if(isprime(n^3+(n+1)^3+(n+3)^3), print1(n, ", "))); \\ Joerg Arndt, Sep 09 2014
(Python)
from sympy import isprime
A246666_list = [n for n in range(1, 10**5) if isprime(3*n*(n*(n+4)+10)+28)]
# Chai Wah Wu, Sep 09 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Carmine Suriano, Sep 01 2014
STATUS
approved