【4.5 扩展欧几里得算法详解】

$更好的阅读体验$

4.5 扩展欧几里得算法

思想

• 欧几里得算法：$gcd(a,b)=gcd(b,a\%b)$，特别的$gcd(a,0)=a$

• 裴蜀定理：对于任意正整数$a,b$，一定存在非零的$x,y$，使得$ax+by=gcd(a,b)$

• $$\{\begin{array}{l}b=0时:\{\begin{array}{l}gcd(a,b)=a\\ ax+by=gcd(a,b)\end{array}\Rightarrow \{\begin{array}{l}x=1\\ y=0\end{array}\\ \\ \\ \\ \\ b!=0时：\{\begin{array}{l}\begin{array}{rl}①& 设ax+by=gcd(a,b)=d\\ & \because 由欧几里得算法可知：gcd(a,b)=gcd(b,a\%b)=d\\ & \therefore 由裴蜀定理得：b{x}^{\prime }+(a\%b)y^{\prime }=d\\ 又& \because ax+by=d\\ & \therefore 联立\{\begin{array}{l}ax+by=d\\ b{x}^{\prime }+(a\%b)y^{\prime }=d\\ a\%b=a-\lfloor \frac{a}{b}\rfloor b\end{array}\Rightarrow \{\begin{array}{l}x=y^{\prime }\\ y=x^{\prime }-\lfloor \frac{a}{b}\rfloor y^{\prime }\end{array}\\ ②& 设a^{\prime }=b,b^{\prime }=a\%b\\ & \therefore gcd(b,a\%b)=gcd(a^{\prime },b^{\prime })=d\\ & \because gcd(a^{\prime },b^{\prime })=gcd(b^{\prime },a^{\prime }\%b^{\prime })=d\\ & \therefore b^{\prime }{x}^{\prime \prime }+a^{\prime }\%b^{\prime }{y}^{\prime \prime }=d\\ 又& \because b{x}^{\prime }+(a\%b)y^{\prime }=d\\ & \therefore 联立\{\begin{array}{l}b{x}^{\prime }+(a\%b)y^{\prime }=d\\ b^{\prime }{x}^{\prime \prime }+a^{\prime }\%b^{\prime }{y}^{\prime \prime }=d\\ a^{\prime }\%b^{\prime }=a^{\prime }-\lfloor \frac{a^{\prime }}{b^{\prime }}\rfloor b^{\prime }\end{array}\Rightarrow \{\begin{array}{l}{x}^{\prime }=y^{\prime \prime }\\ y^{\prime }=x^{\prime \prime }-\lfloor \frac{a^{\prime }}{b^{\prime }}\rfloor y^{\prime \prime }\end{array}\\ ③& 设a^{\prime \prime }=b^{\prime },b^{\prime \prime }=a^{\prime }\%b^{\prime }\\ & \dots \\ & \dots \\ & 直到b=0时，联立解得\{\begin{array}{l}{x}^i=1\\ y^i=0\end{array}\\ & 然后逐步返回每一次联立所得的结果\{\begin{array}{ll}{x}^{i-1}=y^i& \\ y^{i-1}=x^i-\lfloor \frac{a^i}{b^i}\rfloor y^i& 最后返回得到x和y的值\end{array}\end{array}\end{array}\end{array}$$

注意

• 当方程符合$ax+by=gcd(a,b)$的形式时，才可以用扩展欧几里得算法求解$(x_0,y_0)$

4.5.1 欧几里得算法的推论

• 可以进一步求解任意方程$ax+by=n$，得到一个整数解

• $$\begin{array}{r}\{\begin{array}{ll}& (1)判断方程ax+by=n是否有整数解，有解的条件为：gcd(a,b)可以整除n\\ & (2)用扩展欧几里得算法求ax+by=gcd(a,b)得到一个解(x_0,y_0)\\ & (3)在a{x}_0+b{y}_0=gcd(a,b)两边同时乘\frac{n}{gcd(a,b)}\Rightarrow \frac{a{x}_{0}n}{gcd(a,b)}+\frac{b{y}_{0}n}{gcd(a,b)}=n\\ & (4)对照ax+by=n可知该方程的一个解为(x^{\prime },y^{\prime })，其中\{\begin{array}{l}{x}^{\prime }=\frac{x_{0}n}{gcd(a,b)}\\ y^{\prime }=\frac{y_{0}n}{gcd(a,b)}\end{array}\end{array}\end{array}$$

模板

void exgcd(int a,int b,int &x,int &y){
    
    if(!b){  //若b=0时
        x=1,y=0;
        return ;
    }
    else{  //b!=0时
        exgcd(b,a%b,x,y);  //递归到下一层
        int t=x;  //返回时执行
        x=y;
        y=t-a/b*y;
    }
    
}

---

例题 877. 扩展欧几里得算法

原题链接

描述

给定 n 对正整数 ai,bi，对于每对数，求出一组 xi,yi，使其满足 ai×xi+bi×yi=gcd(ai,bi)。

输入格式
第一行包含整数 n。

接下来 n 行，每行包含两个整数 ai,bi。

输出格式
输出共 n 行，对于每组 ai,bi，求出一组满足条件的 xi,yi，每组结果占一行。

本题答案不唯一，输出任意满足条件的 xi,yi 均可。

数据范围
1≤n≤105,
1≤ai,bi≤2×109
输入样例：

2
4 6
8 18

输出样例：

-1 1
-2 1

代码

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

void exgcd(int a,int b,int &x,int &y){
    
    if(!b){
        
        x=1,y=0;
        return ;
        
    }
    else{
        
        exgcd(b,a%b,x,y);
        
        int t=x;
        
        x=y;
        
        y=t-a/b*y;
        
    }
    
}

int main(){
    
    int n;
    
    cin>>n;
    
    while(n--){
        
        int a,b,x,y;
        
        cin>>a>>b;
        
        exgcd(a,b,x,y);
        
        cout<<x<<" "<<y<<endl;
        
    }
    
    return 0;
        
}

---

4.5.2 解一元线性同余方程

---

概念

• $ax\equiv b(modm)$，即$\frac{ax}{m}$与$\frac{b}{m}$的余数相同，且$a,b,m$为整数，求$x$的值
• 该方程即为一元线性同余方程

思想

• 对$ax\equiv b(modm)$做等价变形：$ax+my=b$

• $$\begin{array}{rlrl}& \because & ax& \equiv b(modm)\\ & \therefore & ax\%m& =k(b\%m),(k\in \mathbb{Z})\\ & \therefore & ax-\lfloor \frac{ax}{m}\rfloor m& =k(b-\lfloor \frac{b}{m}\rfloor m)\\ & \therefore & ax-kb& =(\lfloor \frac{ax}{m}\rfloor -k\lfloor \frac{b}{m}\rfloor )m\\ & \because & \lfloor \frac{ax}{m}\rfloor ,& \lfloor \frac{b}{m}\rfloor ,k\in \mathbb{Z}\\ & \therefore & (\lfloor \frac{ax}{m}\rfloor & -k\lfloor \frac{b}{m}\rfloor )\in \mathbb{Z}\\ & & 设(\lfloor \frac{ax}{m}\rfloor & -k\lfloor \frac{b}{m}\rfloor )=y,(y\in \mathbb{Z})\\ & \therefore & ax-kb& =my\Rightarrow ax-my=b\\ 又& \because & y& 可以为负数\\ & \therefore & ax\equiv b(mod& m)\leftrightarrow ax+my=b\end{array}$$

• 由扩展欧几里得算法的推论可知：当且仅当$gcd(a,m)$可以整除$b$时，$ax+my=b$存在整数解

• $$由扩展欧几里得算法可知：\{\begin{array}{l}当gcd(a,m)=b时：\{\begin{array}{l}x=x_0\\ y=y_0\end{array}\\ 当gcd(a,m)为b的整数倍时：\{\begin{array}{l}{x}^{\prime }=\frac{x_{0}b}{gcd(a,m)}\\ y^{\prime }=\frac{y_{0}b}{gcd(a,m)}\end{array}\end{array}$$

---

例题 878. 线性同余方程

原题链接

描述

给定 n 组数据 ai,bi,mi，对于每组数求出一个 xi，使其满足 ai×xi≡bi(modmi)，如果无解则输出 impossible。

输入格式
第一行包含整数 n。

接下来 n 行，每行包含一组数据 ai,bi,mi。

输出格式
输出共 n 行，每组数据输出一个整数表示一个满足条件的 xi，如果无解则输出 impossible。

每组数据结果占一行，结果可能不唯一，输出任意一个满足条件的结果均可。

输出答案必须在 int 范围之内。

数据范围
1≤n≤105,
1≤ai,bi,mi≤2×109

2
2 3 6
4 3 5

输出样例：

输出样例：

impossible
-3

代码

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

typedef long long LL;

void exgcd(LL a,LL b,LL &x,LL &y){
    
    if(!b){  //若b=0时
        x=1,y=0;
        return ;
    }
    else{
        exgcd(b,a%b,x,y);
        LL t=x;
        x=y;
        y=t-a/b*y;
    }
    
}

LL gcd(LL a,LL b){
    return b ? gcd(b,a%b) : a;
}

int main(){
    
    int n;
    
    cin>>n;
    
    while(n--){
        
        LL a,b,m,x,y;
        
        cin>>a>>b>>m;
        
        LL d=gcd(a,m);
        
        exgcd(a,m,x,y);
        
        if(b%d) cout<<"impossible"<<endl;
        else cout<<b/d*x%m<<endl;
        
    }
    
    return 0;
    
}

---