同余方程(extend_gcd&&线性同余方程)

求关于x的同余方程 ax ≡ 1(mod b) 的最小正整数解。

输入格式

输入只有一行，包含两个正整数a,b,用一个空格隔开。

输出格式

输出只有一行,包含一个正整数x，表示最小正整数解。

输入数据保证一定有解。

数据范围

2≤a,b≤2∗109

输入样例：

3 10

输出样例：

7

思路：因为  ax ≡ 1(mod b)的含义为ax%b=1%b，移项化简得(ax-1)%b==0,则ax-1=bt，即：ax-bt=1，令y为一负数，使得ax+by=gcd（a，b），这样就得到了extend_gcd的形式，用extend_gcd可以求出一个特解x，因为b%b=0，所以(x%b+b)%b后对(ax-1)%b==0无影响，即不影响原同余方程的正确性，因为用extend_gcd求出的一个特解x可能是负值，所以为了得到最小正整数解最后用extend_gcd可以求出一个特解x后让(x%b+b)%b即可得到答案。

完整代码：

#include <iostream>
#define int long long 

using namespace std;

int a,b,x,y;

inline int extend_gcd(int a,int b,int &x,int &y)
{
    if(b==0){
        x=1,y=0;
        return a;
    }
    else{
        int res=extend_gcd(b,a%b,y,x);
        y-=x*(a/b);
        return res;
    }
}

signed main()
{
    cin>>a>>b;
    extend_gcd(a,b,x,y);
    cout<<(x%b+b)%b<<endl;
    return 0;
}