扩展欧几里得算法 AcWing 878. 线性同余方程

扩展欧几里得算法 AcWing 878. 线性同余方程

原题链接

AcWing 878. 线性同余方程

算法标签

数学知识 同余 扩展欧几里得算法

思路

该图转自该题解
住： 1 Bezout即裴蜀定理
2 联立$a\ast x+m\ast y=b$与$a\ast x_0+m\ast y_0=gcd(a,m)$可得x

代码

#include<bits/stdc++.h>
#define int long long
#define rep(i, a, b) for(int i=a;i<b;++i)
#define Rep(i, a, b) for(int i=a;i>b;--i)
using namespace std;
int x, y;
inline int read(){
   int s=0,w=1;
   char ch=getchar();
   while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
   while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar();
   return s*w;
}
void put(int x) {
    if(x<0) putchar('-'),x=-x;
    if(x>=10) put(x/10);
    putchar(x%10^48);
}
int exgcd(int a, int b, int &x, int &y){
    if(!b){
        x=1, y=0;
        return a;
    }
    int d=exgcd(b, a%b, y, x);
    y-=a/b*x;
    return d;

}
signed main(){
	ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	int n=read();
	while(n--){
	    int a=read(), b=read(), m=read();
	    int d=exgcd(a, m, x, y);
	    if(b%d){
	        puts("impossible");
	    }else{
	        printf("%lld\n", x*b/d%m);
	    }
	}
	return 0;
}

模板

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


原创不易
转载请标明出处
如果对你有所帮助 别忘啦点赞支持哈
![在这里插入图片描述](https://img-blog.csdnimg.cn/9cb4cb80a18345158f8a69e2f196e5b9.png)