1. 程式人生 > >擴充套件gcd模板

擴充套件gcd模板

void exgcd(ll a, ll b, ll &x, ll &y){
    if(b == 0){
        x = 1;
        y = 0;
        return;
    }
    ll x1, y1;
    exgcd(b, a%b, x1, y1);
    x = y1;
    y = x1-(a/b)*y1;
}

求的是ax+by=gcd(a,b),對於ax+by=c。要先判斷c是否能整除gcd(a,b),否則無解。 例如求逆元。 ax≡1(mod)->ax +b*mod = 1

//求a的逆元
	exgcd(a, mod, x,
y); while(x<0) x += mod; printf("%lld\n", x%mod);//因為mod是質數,一定與 ``