OFFSET
1,1
COMMENTS
Previous name: Primes of the form q = p^6 + p^3 + 1 = sigma_3(p^2), where p is prime.
LINKS
Harry J. Smith, Table of n, a(n) for n = 1..200
EXAMPLE
73 = 64 + 8 + 1, 757 = 729 + 27 + 1, etc.
PROG
(PARI) { n=0; for (m=1, 10^9, p=prime(m); q=p^6 + p^3 + 1; if (isprime(q), write("b066101.txt", n++, " ", q); if (n==200, return)) ) } \\ Harry J. Smith, Nov 13 2009
CROSSREFS
KEYWORD
dead
AUTHOR
Labos Elemer, Dec 04 2001
STATUS
approved