OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
P:= x^5 - x^4 - x^3 - x^2 - x - 1:
select(p -> Irreduc(P) mod p, [seq(ithprime(i), i=1..1000)]); # Robert Israel, Mar 13 2024
MATHEMATICA
P = x^5 - x^4 - x^3 - x^2 - x - 1;
Select[Prime[Range[1000]], IrreduciblePolynomialQ[P, Modulus -> #]&] (* Jean-François Alcover, Mar 24 2024, after Robert Israel *)
PROG
(Python)
from itertools import islice
from sympy import Poly, nextprime
from sympy.abc import x
def A371566_gen(): # generator of terms
p = 2
while True:
if Poly(x*(x*(x*(x*(x-1)-1)-1)-1)-1, x, modulus=p).is_irreducible:
yield p
p = nextprime(p)
(PARI) a371566(upto) = forprime (p=2, upto, my(f=factormod(x^5 - x^4 - x^3 - x^2 - x - 1, p)); if(#f[, 1]==1, print1(p, ", "))) \\ Hugo Pfoertner, Mar 22 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert Israel, Mar 27 2024
STATUS
approved