factor X^3+X+1 mod 31