shithub: riscv

ref: b4a9468fdd4674b838ece2c448a7f81a9c330899
dir: /sys/src/libmp/port/mpinvert.c/

View raw version
#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);
}