扩展欧几里得

扩展欧几里得

食用范围

给定 $n$ 对正整数$a_i,b_i$，对于每对数，求出一组 $x_i,y_i$ ，使其满足 $a_i\ast x_i+b_i\ast y_i=gcd(a_i,b_i)$

此等式存在的条件：$gcd(a,m)|b$ 成立

裴蜀定理

对于任意正整数 $a,b$ ，那么一定存在非零整数 $x,y$ 使得 $ax+by=min(gcd(a,b))$

求这个 $x,y$ 的过程，叫扩展欧几里得算法

欧几里得算法的基本原理：$gcd(a,b)=gcd(b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$

推导过程

注：$a\mathrm{m}\mathrm{o}\mathrm{d}b=a-\lfloor \frac{a}{b}\rfloor \times b$

Code

ll exgcd(ll a, ll b, ll &x, ll &y) {
    //return gcd(a,b) and get x,y for ai*x+bi*y=gcd(ai,bi)
	if(b == 0) {
		x = 1, y = 0;
		return a;
	}
	int d = exgcd(b, a % b, y, x);
	y -= a / b * x;
	return d;
}

使用扩展欧几里得求解线性同余方程

问题

求 $x$ $s.t.(使得)$ $ax\equiv b\mathrm{m}\mathrm{o}\mathrm{d}m$

$ax\equiv b\mathrm{m}\mathrm{o}\mathrm{d}m\iff ax=my+b$

$\iff ax-my=b$

令 $y=y^{\prime }$

$\iff ax-m{y}^{\prime }=b$

使用扩展欧几里得 （成立条件：$gcd(a,m)|b$）求出 $x$

再把 $x$ 放大 $b/d$ 倍就得到了答案 $d=gcd(a,m)$

Code

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;

ll exgcd(ll a, ll b, ll &x, ll &y) {//return gcd(a,b) and get x,y for ai*x+bi*y=gcd(ai,bi)。
	if(b == 0) {
		x = 1, y = 0;
		return a;
	}
	int d = exgcd(b, a % b, y, x);
	y -= a / b * x;
	return d;
}

int main() {
	int n;
	scanf("%d", &n);
	while(n--) {
		ll a, b, m;
		scanf("%d%d%d", &a, &b, &m);
		ll x, y;
		ll d = exgcd(a, m, x, y);
		if(b % d) puts("impossible");
		//因为b是d的b/d倍，所以最终答案还要转化为在b意义下的答案，也就是把d扩大b/d倍
		else printf("%lld\n", x * (b / d) % m);

	}

	return 0;
}

总结

求逆元的两种方式

• 如果模数是质数 ，快速幂 + 费马小定理
• 如果模数不是质数，$exgcd$