欧几里得算法和扩展欧几里得算法
欧几里得算法
说白了就是辗转相除算法
我们在求a和b的gcd的时候
首先我们会发现
gcd(a,b)=gcd(a−b,b) (a>=b)
g
c
d
(
a
,
b
)
=
g
c
d
(
a
−
b
,
b
)
(
a
>=
b
)
那么我们不停地迭代就发现
gcd(a,b)=gcd(a mod b,b)
g
c
d
(
a
,
b
)
=
g
c
d
(
a
m
o
d
b
,
b
)
然后我们就可以不停地迭代下去
当小的数等于0的时候答案就是大的那个数
最后写出来就是
int gcd(int a,int b){
if(!b)return a;
return gcd(b,a%b);
}
扩展欧几里得
然后我们需要求的从
gcd(x,y)
g
c
d
(
x
,
y
)
变成了
ax+by=c
a
x
+
b
y
=
c
也就是不定方程求解
然后考虑一个不定方程
ax+by=c
a
x
+
b
y
=
c
有解的充要条件就是
gcd(a,b)|c
g
c
d
(
a
,
b
)
|
c
即,我们想求出满足
ax+by=gcd(a,b)
a
x
+
b
y
=
g
c
d
(
a
,
b
)
的解
然后惊讶的发现
gcd(a,b)=gcd(b,a mod b)
g
c
d
(
a
,
b
)
=
g
c
d
(
b
,
a
m
o
d
b
)
然后我们用
b
b
和代替
a
a
和
发现
gcd(a,b)=bx+(a mod b)y
g
c
d
(
a
,
b
)
=
b
x
+
(
a
m
o
d
b
)
y
又因为
a mod b=a−(a/b)∗b
a
m
o
d
b
=
a
−
(
a
/
b
)
∗
b
(这里是整除)
所以
gcd(a,b)=bx+(a−(a/b)∗b)y
g
c
d
(
a
,
b
)
=
b
x
+
(
a
−
(
a
/
b
)
∗
b
)
y
整理一下
gcd(a,b)=ay+b(x−(a/b)∗y)
g
c
d
(
a
,
b
)
=
a
y
+
b
(
x
−
(
a
/
b
)
∗
y
)
又惊讶地发现:
x−>y
x
−
>
y
y−>x−(a/b)∗y
y
−
>
x
−
(
a
/
b
)
∗
y
然后就可以利用exgcd的过程进行求解了
void exgcd(int a,int b,int &x,int &y){
if(!b)x=1,y=0;
else{
exgcd(b,a%b,y,x);
y-=(a/b)*x;
}
}
但是问题又来了
我们怎么求出最小的正整数解呢
首先我们还原一下式子,把
x,y
x
,
y
同时乘上
c/gcd(a,b)
c
/
g
c
d
(
a
,
b
)
又因为通解可以表示成:
x=x′+bgcd(a,b)∗k
x
=
x
′
+
b
g
c
d
(
a
,
b
)
∗
k
y=y′−agcd(a,b)∗k
y
=
y
′
−
a
g
c
d
(
a
,
b
)
∗
k
所以x的最小正整数解就是
(x%bgcd(a,b)+bgcd(a,b))%bgcd(a,b)
(
x
%
b
g
c
d
(
a
,
b
)
+
b
g
c
d
(
a
,
b
)
)
%
b
g
c
d
(
a
,
b
)