同余方程(扩展欧几里得)

已知整数 a 和 $b$，求关于 x 的同余方程 ax≡1(mod b) 的最小正整数解。

输入格式

输入一行，输入两个整数 a,b（ 2≤a,b≤2×10的9次方）。

输出格式

输出一行，输出一个整数，即同余方程的最小正整数解。输入数据保证一定有解。

样例输入

3 8

样例输出

3

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

首先要知道 ax≡1(mod b)的意思是 (ax-1)%b=1，那么就可以转换方程ax+bk=1;

并且a关于p的逆元存在的条件是gcd(a,p)=1,即a与p互质；

扩展欧几里得是用来算ax+by=gcd(a,b)的解的，所以可以把bk看成by，方程为ax+by=1；

所以直接用扩展欧几里得算，注意！！！但最后求出的x有可能为负数，所以输出时把x换成(x+b)%b;若x是正数，则不变，若为负数，则变成最小正数解；