扩展欧几里得算法

基本知识

裴蜀定理：对于任意一对正整数ab，则一定存在xy，是的ax+by=gcd(a,b)
扩展欧几里得算法就是一种构造xy的算法

例题

扩展欧几里得算法

#include <bits/stdc++.h>
using namespace std;
int exgcd(int a, int b, int &x, int &y) //在欧几里得算法的基础上进行扩展，算出xy
{
    if (!b)
    {
        x = 1, y = 0;
        return a;
    }
    int d = exgcd(b, a % b, y, x);
    y -= a / b * x;
    return d;
}
int main()
{
    int n;
    cin >> n;
    while (n--)
    {
        int a, b;
        int x, y;
        cin >> a >> b;
        int l = exgcd(a, b, x, y);
        cout << x << " " << y << endl;
    }
    return 0;
}

由上述代码可以求出一个xy的组合
由公式可以推知所有的解：x=x0-b/ak，y=y0+b/ak这样就可以得到所有的解

线性同余方程
由题目内容做等价变换可以得到ax+my=b，其中y为待定系数，根据题中条件如果b不是gcd(a,m)的倍数，则一定没有解，否则，一定有解

#include <bits/stdc++.h>
using namespace std;
int exgcd(int a, int b, int &x, int &y) //在欧几里得算法的基础上进行扩展，算出xy
{
    if (!b)
    {
        x = 1, y = 0;
        return a;
    }
    int d = exgcd(b, a % b, y, x);
    y -= a / b * x;
    return d;
}
int main()
{
    int n;
    cin >> n;
    while (n--)
    {
        int a, b, m;
        cin>>a>>b>>m;
        int x,y;
        if(b%d) cout<<"impossible"<<endl;
        else cout<<(long long)x*(b/d)%m<<endl;
    }
    return 0;
}