###methode um die wurzeln von f mod p  auszurechnen
###input: f polynom ueber Z, p prim
###output: nullstellen von f mod p
nst_mod_p := function(f,p)
local k,r;
k := PolynomialAlgebra(GF(p));
f := Coerce(k,f);
r := Roots(f);
return Apply(r, i-> Coerce(Z,i[1]));
end;