扩展欧几里得

最新推荐文章于 2024-06-02 14:55:04 发布

wtglndpyyzqq

最新推荐文章于 2024-06-02 14:55:04 发布

阅读量135

点赞数

本文链接：https://blog.csdn.net/weixin_48700632/article/details/112973863

版权

GCD算法：求最大公约数
重要等式gcd(a,b)=gcd(b,a mod b)
扩展的欧几里得算法：扩展欧几里得算法本身过程求的是 ax + by = gcd(a,b)的解。通过辗转相除的方法向下递进，不断缩小系数，保证会出现有确定解的最后一层。
具体步骤如下
gcd（a，b）=gcd(b,a mod b)
有x1，y1使 bx1+(a mod b)y2=gcd（b, a mod b）
让 ax+by=bx1+（a mod b）y2
又因为 a mod b=a-(a/b)b
所以ax+by=bx1+(a-(a/b)b)y1=bx1+ay1+a/bby1
=ay1+b*(x1-(a/b)y1)
解 x=y1 , y=x1-a/by1
递归直到b=0

落谷p1028 同余方程
代码如下

#include <bits/stdc++.h>
using namespace std;
int a,b,k,x,y;
void exgcd(int a,int b){
	if(b==0){ x=1; y=0; return;}
	exgcd(b,a%b);
	k=x;
	x=y;
	y=k-a/b*y;
	return;
}
int main(){
	scanf("%d%d",&a,&b);
	exgcd(a,b);
	printf("%d",(x+b)%b);
}