http://benji3up2kxewkqfcq7buxk2xd6zwy3zggnurkrm3l4cvwy2iipvyyad.onion/mirrors/gmpdoc/Low_002dlevel-Functions.html
Function: void mpn_tdiv_qr (mp_limb_t * qp , mp_limb_t * rp , mp_size_t qxn , const mp_limb_t * np , mp_size_t nn , const mp_limb_t * dp , mp_size_t dn ) ¶ Divide { np , nn } by { dp , dn } and put the quotient
at { qp , nn − dn +1} and the remainder at { rp , dn }.