ref: 45ce4de65064ffb0e8b1904118cea7e16528cffa
dir: /sys/src/libmp/port/mpinvert.c/
#include "os.h" #include <mp.h> // use extended gcd to find the multiplicative inverse // res = b**-1 mod m void mpinvert(mpint *b, mpint *m, mpint *res) { mpint *v; v = mpnew(0); mpextendedgcd(b, m, v, res, nil); if(mpcmp(v, mpone) != 0) abort(); mpfree(v); mpmod(res, m, res); }