Select[Prime[Range[100]],IrreduciblePolynomialQ [x^5+x^3+x+1,Modulus->#]&]