标签：a% gcd int 欧几里得 扩展 算法 tag end aligned
Bézout’s定理
∀
a
,
b
∈
Z
\forall a, b \in Z
∀a,b∈Z，
∃
x
,
y
∈
Z
+
\exists x,y \in Z^+
∃x,y∈Z+，使得
a
x
+
b
y
=
g
c
d
(
a
,
b
)
(0)
ax + by = gcd(a,b) \tag{0}
ax+by=gcd(a,b)(0)
扩展欧几里得算法求该定理的解
已知
g
c
d
(
a
,
b
)
=
g
c
d
(
b
,
a
%
b
)
(1)
\begin{aligned} gcd(a, b) = gcd(b, a \% b) \end{aligned} \tag{1}
gcd(a,b)=gcd(b,a%b)(1)
根据
B
e
ˊ
z
o
u
t
′
s
Bézout's
Beˊzout′s定理，
∃
x
′
,
y
′
∈
Z
+
\exists x^{'}, y^{'} \in Z^+
∃x′,y′∈Z+，使得：
b
x
′
+
(
a
%
b
)
y
′
=
g
c
d
(
b
,
a
%
b
)
(2)
\begin{aligned} bx^{'} + (a \% b)y^{'} = gcd(b, a \% b) \end{aligned} \tag{2}
bx′+(a%b)y′=gcd(b,a%b)(2)
又
a
%
b
=
a
−
b
∗
⌊
a
/
b
⌋
(3)
\begin{aligned} a \% b = a - b * \lfloor a / b \rfloor \end{aligned} \tag{3}
a%b=a−b∗⌊a/b⌋(3)
将
(
3
)
(3)
(3)带入
(
2
)
(2)
(2)，得：
b
x
′
+
(
a
−
b
∗
⌊
a
/
b
⌋
)
y
′
=
g
c
d
(
b
,
a
%
b
)
a
y
′
+
b
(
x
′
−
⌊
a
/
b
⌋
y
′
)
=
g
c
d
(
b
,
a
%
b
)
(4)
\begin{aligned} bx^{'} + (a - b * \lfloor a / b \rfloor) y^{'} &= gcd(b, a \% b) \\ ay^{'} + b(x^{'} - \lfloor a / b \rfloor y^{'}) &= gcd(b, a \% b) \end{aligned} \tag{4}
bx′+(a−b∗⌊a/b⌋)y′ay′+b(x′−⌊a/b⌋y′)=gcd(b,a%b)=gcd(b,a%b)(4)
令
x
=
y
′
x = y^{'}
x=y′，
y
=
x
′
−
⌊
a
/
b
⌋
y
′
y= x^{'} - \lfloor a / b \rfloor y^{'}
y=x′−⌊a/b⌋y′，将
(
1
)
(1)
(1)代入
(
4
)
(4)
(4)，整理得：
a
x
+
b
y
=
g
c
d
(
b
,
a
%
b
)
=
g
c
d
(
a
,
b
)
(5)
\begin{aligned} ax + by = gcd(b, a \% b) = gcd(a,b) \end{aligned} \tag{5}
ax+by=gcd(b,a%b)=gcd(a,b)(5)
(
5
)
(5)
(5)说明，知道
(
2
)
(2)
(2)的解就可以推出定理的解。公式(0)到公式(2)明显是欧几里得算法求解最大公约数的过程。公式(3)，(4)，(5)是说明：公式(0)的解与公式(2)的解的关系。将上述过程归纳成代码形式：
void exgcd(int a, int b, int& x, int& y) { //公式(0)
if (b == 0) {
x = 1;
y = 0;
return a;
}
exgcd(b, a % b, x, y); //公式(2)
int tp = x;
x = y;
y = tp - a / b * y;
}
参考
- https://www.cnblogs.com/fusiwei/p/11775503.html
标签：a%,gcd,int,欧几里得,扩展,算法,tag,end,aligned 来源： https://blog.csdn.net/oooooorz/article/details/122326348
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。