数学：乘法逆元-拓展GCD

乘法逆元应用在组合数学取模问题中，这里给出的实现不见得好用

给出拓展GCD算法：

扩展欧几里得算法是指对于两个数a,b
一定能找到x,y（均为整数，但不满足一定是正数）
满足x*a+y*b=gcd(a,b)
gcd（x,y)是指x 与 y的最大公约数

有啥用呢？求解形如 a*x +b*y = c 的通解

然后我们先介绍同余方程，再介绍乘法逆元

同余方程
a≡b(mod m) 等价于小学的运算式 b÷m 余数为a
也就是a mod m=b

其实介绍这个就是看怎么把≡拿掉

乘法逆元
ax ≡ 1 (mod m)
我们称 x 是 a 关于 m 的乘法逆元
可以等价于这样的表达式： a*x + m*y = 1

当满足这个式子的时候：a*x + b*y = c 有解的充要条件： c % gcd(a , b) == 0

一般，我们能够找到无数组解满足条件，但是一般是让你求解出最小的那组解

我们求解出来了一个特殊的解 x0 ，我们用 x0 % m其实就得到了最小的解了

#include<cstdio>
using namespace std;
inline long long read()
{
    long long x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9') {if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int a,b;
void exgcd(int a,int b,int &x,int &y)
{
    if(b==0) {x=1;y=0;return;}
    exgcd(b,a%b,x,y);
    int t=x;x=y;y=t-a/b*y;
}
//ax ≡ 1 (mod b)
//-> a*x + b*y = 1
//->求出x和y后让x%b就是最小解了 
int main()
{
    a=read();b=read();
    int x,y;
    exgcd(a,b,x,y);
    x=(x%b+b)%b;
    printf("%d",x);
    return 0;
}