欧几里得算法与 EX

合集 - 【算法学习】(35)

1.【算法学习】排序2024-08-272.【算法学习】Manacher 马拉车2024-08-293.【算法学习】KMP 算法2024-09-014.LCA 最近公共祖先（树链和倍增）这次真有树链了！！！2024-08-305.线段覆盖问题2024-08-296.【算法学习】学换根dp有感2024-09-207.二分图最大匹配2024-09-178.【算法学习】01BFS2024-09-179.洛谷 P1892 [BOI2003] 团伙 种类并查集 扩展域并查集2024-09-1110.【算法学习】高斯消元法2024-09-1011.贝叶斯公式2024-09-0712.背包2024-09-0513.【算法学习】模拟退火2024-09-0514.【算法学习】基环树2024-10-2915.【算法学习】树链部分2024-10-2916.【算法学习】莫队2024-10-2917.【算法学习】分块九讲2024-10-2818.平衡树2024-10-2319.圆方树2024-10-2220.【算法学习】点分治2024-10-1121.公式2024-10-0922.【算法学习】笛卡尔树2024-10-0423.【算法学习】悬线法2024-10-04

24.欧几里得算法与 EX2024-09-23

25.【算法学习】逆元与求解2024-09-2326.【算法学习】费马定理2024-09-2327.三分2024-09-2328.裴蜀定理2024-09-2329.【算法学习】欧拉函数φ2024-09-2330.【算法学习】二维转一维问题2024-11-0331.【算法学习】扫描线2024-11-0132.【算法学习】矩阵乘法2024-10-3133.【算法学习】同余最短路2024-11-0434.【算法学习】组合数学2024-11-0535.【算法学习】反悔贪心2024-11-06

收起

欧几里得算法

欧几里得算法又称辗转相除法，用来求两个数的最大公约数的算法。

省流：\(gcd(a,b)=gcd(b,a\mod b)\)

学习笔记

点击查看代码

#include <bits/stdc++.h>
using namespace std;
int a,b;
int gcd(int x,int y){
	return y==0?x:gcd(y,x%y);
}
int main(){
    cin>>a>>b; 
    cout<<gcd(a,b);
    return 0;
}

拓展欧几里得算法 exgcd

学习笔记

学习笔记

常用于求解 \(ax+by=gcd(a,b)\) 的一组可行解。

接下来就是推例子了：

\[\because gcd(a,b)=gcd(b,a\mod b) \]

\[\therefore ax+by=bx_2+(a \mod b)y_2 \]

\[\because a \mod b=a-b\times \lfloor \frac{a}{b} \rfloor \]

\[\therefore ax+by=bx_2+( a-b\times \lfloor \frac{a}{b} \rfloor )y_2 \]

\[\therefore ax+by=bx_2+ay_2-by_2\lfloor \frac{a}{b} \rfloor \]

\[\therefore ax+by=ay_2+b(x_2-y_2\lfloor \frac{a}{b} \rfloor) \]

然后我们就发现 \(x=y_2,y=x_2-y_2\lfloor \frac{a}{b} \rfloor\)

\(x_2,y_2\) 又可以从 \(x_3,y_3\) 推过来。

我们就需要不断做欧几里得法，当 \(a=gcd(a,b),b=0\) 时，我们这时取 \(x_n=1,y_n=0\)，因为 \(b=0\) 所以 \(y_n\) 取任何数都没问题，但取 \(0\) 的时候 \(x\) 和 \(y\) 增长速度更慢一些，还有注意解可能为负数的情况，所以要取模变正数。

点击查看代码

#include <bits/stdc++.h>
using namespace std;
#define ll long long
void exgcd(ll a,ll b,ll &x,ll &y){
	if(b==0){
		x=1,y=0;
		return;
	}
	exgcd(b,a%b,y,x);
	y-=a/b*x;
}
int main(){
	ios::sync_with_stdio(false);
	ll a,b;
	cin>>a>>b;
	ll d=__gcd(a,b);
	if(b%d!=0){
		cout<<"impossible\n";
	}
	else{
		ll x,y;
		exgcd(a,b,x,y);
		x=(x*(b/d)%b+b)%b;//注意这要调整解
		cout<<x<<"\n";
	}
    return 0;
}