OFFSET
1,1
COMMENTS
n such that n+1 is a prime with primitive root 3 (A019334 except for the first term). [Joerg Arndt, Jul 05 2011]
PROG
(PARI) for(n=2, 1000, if(polisirreducible(Mod(1, 3)*sum(e=0, n, x^e)), print1(n+1, ", "))) /* Joerg Arndt, Jul 05 2011 */
(PARI) forprime(p=5, 1000, if(znorder(Mod(3, p))==p-1, print1(p-1, ", "))) /* much faster */ /* Joerg Arndt, Jul 05 2011 */
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Robert G. Wilson v, Jun 22 2002
STATUS
approved