FactorSquareFree[(x^17) + 1, Modulus -> 17]

a = (x^5) + 4  (x^2) + 9  x + 2

FactorSquareFree[a, Modulus -> 17]