ref: c186019e2661ec8b8c301e5d7dc1bcd99f71f48c
dir: /3rd/mp/mpinvert.c/
#include "platform.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){ fprintf(stderr, "mpinvert: impossible\n"); exit(2); } mpfree(v); mpmod(res, m, res); }