逆元模板

一、作用

(a +  b) % p = (a%p +  b%p) %p 
(a  -  b) % p = (a%p  -  b%p) %p 
(a  *  b) % p = (a%p *  b%p) %p 
(a  /  b) % p = (a%p  /  b%p) %p 

因为上三条是同余公式，正确，但是第四条满足左右恒等关系，但是又存在计算过程中爆精度的问题，所以有了逆元的思想

a关于p的逆元表示为inv（a），同时a与p互质

(a  /  b) % p = (a * inv(a) ) % p = (a % p * inv(a) % p) % p

便可以转换成乘法计算；

二、代码

int extgcd(int a, int b, int& x, int& y)//欧几里得算法
{
    int d = a;
    if(b != 0){
        d = extgcd(b, a % b, y, x);
        y -= (a / b) * x;
    }else {
        x = 1;
        y = 0;
    }
    return d;
}
int mod_inverse(int a, int m)
{
    int x, y;
    extgcd(a, m, x, y);
    return (m + x % m) % m;
}
 

一般m选择mod，且为质数

三、证明：

（一）、扩展欧几里得算法

存在整数对（x，y），使得ax+by=gcd（a，b）

若a，b互质，gcd（a,b）=1;

（二）、证明x是a关于b的逆元

a*x % b + b*y % b = 1 % b

a*x % b = 1 % b

a*x = 1 (mod b)

是不是类似于a*x=1,x=1/a。

所以扩展欧几里得算法x的值便是表示a关于b的逆元

（三）、至于为啥要返回（m+x%m）%m我不是很清楚

由于x%m已经为负数（即溢出了）时，再直接对m取模肯定是错的，所以通过加一个m改变符号。

四、应用

卡塔兰大数

ll extgcd(ll a, ll b, ll& x, ll& y)
{
    ll d = a;
    if(b != 0){
        d = extgcd(b, a % b, y, x);
        y -= (a / b) * x;
    }else {
        x = 1;
        y = 0;
    }
    return d;
}

void pre()
{
    int i;
    ll x, y;
    Catelan[0] = 1, Catelan[1] = 1;
    for(i = 2; i < N-5; i++)
    {
        Catelan[i] = Catelan[i-1]*(4*i-2) % mod;
        extgcd(i+1, mod, x, y);
        Catelan[i] = (Catelan[i]*((x+mod)%mod)) % mod;
    }
}