OFFSET
1,1
COMMENTS
Phi(7,x) =1+x+x^2+x^3+x^4+x^5+x^6 =A053716(x) is the 7th cyclotomic polynomial.
EXAMPLE
616067011 = 29^6+29^5+29^4+29^3+29^2+29+1 (29 is prime) and 616067011^6+616067011^5+616067011^4+616067011^3+616067011^2+616067011+1 = 54672347801779330810964871392077416495507203132755717 is prime. Thus, 616067011 is a member of this sequence.
MAPLE
for k from 1 do
p := ithprime(k) ;
n := numtheory[cyclotomic](7, p) ;
pn := numtheory[cyclotomic](7, n) ;
if isprime( pn) then
print(n) ;
end if;
end do: # R. J. Mathar, Feb 07 2014
PROG
(Python)
import sympy
from sympy import isprime
{print(n**6+n**5+n**4+n**3+n**2+n+1) for n in range(10**5) if isprime(n) and isprime((n**6+n**5+n**4+n**3+n**2+n+1)**6+(n**6+n**5+n**4+n**3+n**2+n+1)**5+(n**6+n**5+n**4+n**3+n**2+n+1)**4+(n**6+n**5+n**4+n**3+n**2+n+1)**3+(n**6+n**5+n**4+n**3+n**2+n+1)**2+(n**6+n**5+n**4+n**3+n**2+n+1)+1)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Derek Orr, Feb 06 2014
STATUS
approved