EXGCD (拓展欧几里得)
引入
解不定方程Ax+By=K(得到的x和y只是其中一组解),已知gcd(a,b)=K,
给出A、B、K,求出x和y,满足Ax+By=K。
暴力CODEING
int x,y;
void solve (int a, int b, int c) {
for (int i = -maxn; i <= maxn; i ++) {
for (int j = -maxn; j <= maxn; j ++) {
if (a * i + b * j == c) {
x = i;
y = j;
break;
}
}
}
}
很好打,但是,时间它不允许我这么打。时间复杂度:O( 4 × m a x n 2 4 \times maxn^2 4×maxn2)。
证明思路
正解
二元一次不定方程求解和欧几里得算法到底有什么关系呢?
推导过程:
题面 : 求解 a × x + b × y = c a \times x + b \times y = c a×x+b×y=c。
∵ \because ∵ g c d ( a , b ) = c gcd(a,b) = c gcd(a,b)=c
∴ \therefore ∴ 设 a × x 1 + b × y 1 = g c d ( a , b ) a \times x1 + b \times y1 = gcd ( a,b ) a×x1+b×y1=gcd(a,b)
∵ \because ∵ x = x 1 × c g c d ( a , b ) , y = y 1 × c g c d ( a , b ) x = x1 \times \displaystyle\frac{c}{gcd(a,b)},y = y1 \times \displaystyle\frac{c}{gcd(a,b)} x=x1×gcd(a,b)c,y=y1×gcd(a,b)c
如果 g c d ( a , b ) gcd ( a,b ) gcd(a,b) 不能整除 c 那么此不定方程无整数解
我们需要求出不定方程 a × x 1 + b × y 1 = g c d ( a , b ) a \times x1 + b \times y1 = gcd ( a,b ) a×x1+b×y1=gcd(a,b