Revision history [back]

My collaborator finally figured out how to do this. For your interest, here is how we do it. Let K be the number field, P the prime ideal, and f the polynomial that we want to factor modulo P.

G = K.factor(P)[0][0].residue_field()
f.change_ring(G).factor()

Hope this helps if someone arrives at the same problem in the future.