-- Function: unsigned long int mpz_remove (mpz_t ROP, mpz_t OP, mpz_t
F)
Remove all occurrences of the factor F from OP and store the
result in ROP. The return value is how many such occurrences were
removed.
Not entirely clear whether ROP=F^k or ROP=OP/F^k