OFFSET
1,1
COMMENTS
Any subsequent terms are > 10^5. - Lucas A. Brown, Dec 11 2022
MATHEMATICA
Reap[ Do[ If[ Factor[x^n - x - 1, Modulus -> 3][[0]] =!= Times, Print[n]; Sow[n]], {n, 2, 3000}]][[2, 1]] (* Jean-François Alcover, Apr 03 2013 *)
Select[Range[1000], IrreduciblePolynomialQ[x^# - x - 1, Modulus -> 3] &] (* Robert Price, Sep 19 2018 *)
PROG
(Sage)
P.<x> = GF(3)[]
for n in range(10^6):
if (x^n-x-1).is_irreducible():
print(n)
(PARI)
for (n=1, 10^6, if ( polisirreducible(Mod(1, 3)*(x^n-x-1)), print1(n, ", ") ) );
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Joerg Arndt, Mar 29 2013
EXTENSIONS
a(35)-a(43) from Lucas A. Brown, Dec 11 2022
STATUS
approved