Allow storing factorization in residue ring contexts (fmpz_mod, nmod, mpn_mod, ...)

[fredrik-johansson]

One of many ways that this would be useful: we can compute inverses modulo prime powers using Hensel lifting instead of computing a full-precision XGCD.