拓展欧几里得—逆元—中国剩余定理—同余

快速幂

模板：

long long mul(int a,int b,int mod)
{
    long long ans=1;
    while(b)
    {
        if(b&1) ans=ans*a%mod;
        a=a*a%mod;//注意爆范围
        b>>=1;
    }
    return ans;
}

快速乘

模板：

int mul(int a,int b,int mod)
{
    int ans=0;
    while(b)
    {
        if(b&1) ans=(ans+a)%mod;//注意爆范围
        a=(a+a)%mod;
        b>>=1;
    }
    return ans;
}

欧几里得算法

欧几里得算法定理
定理: gcd(a,b) = gcd ( b,a mod b)
证明:
a可以表示成a=kb+r，则r=a mod b
1.假设d是a,b的一个公约数，则有
a|d, b|d，而r=a-kb， 因此r|d
因此d是(b,a mod b)的公约数,证明充分性
2.假设d是(b,a mod b)的公约数，则
b|d,r|d，但是a=kb +r
因此d也是(a,b)的公约数,证明必要性
因此(a,b)和(b,a mod b)的公约数是一样的，其最大公约数也必
当b=o时，(a， o)的最大公约数是a；
模板：

int gcd(int a,int b)
{
    return b?gcd(b,a%b):a;
}

拓展欧几里得

解 x，y的方法的理解: .
1、显然当b=0，gcd (a, b) =a。此时x=1，y=0;
也就是1a+0b= a;
2、设ax1+by1=gcd(a,b);
在下一个exgcd里 : bx2+(a mod b)y2=gcd(b,a mod b);
根据欧几里德原理有 : gcd(a,b)=gcd(b,a mod b);
则: ax1+by1=bx2+(a mod b)y2;
即: ax1+by1=bx2+(a-(a/b)b)y2=ay2+bx2-(a/b)by2;
根据恒等定理得:
x1=y2; y1=x2-(a/b)*y2;
这样我们就得到解x1,y1的方法: x1，y1的值基于X2，y2。

模板：

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

同余定理

a ≡ b ( mod m) <–> m | (a-b);
线性运算：如果a≡b(mod m)，c≡d(mod m)，那么a ± c≡b ± d(mod m)，且a * c≡b * d(mod m);
幂运算：如果a ≡ b (mod m)，那么a^n ≡ b^n (mod m);

逆元

当n为质数时，可以用快速幂求逆元:
a/b≡ a * x (mod n)
两边同乘b可得a=a* b* x (mod n)
即1≡ b* x (mod n)
同 b*x≡ 1 (mod n)
费马小定理可知，当n为质数时
b^(n- 1)≡ 1 (mod n)
拆一个b出来可得b *b^(n-2)= 1 (mod n)
故当n为质数时，b的乘法逆元 x=b^(n-2);
模板：

#include<iostream>
#include <cmath>
using namespace std;
typedef long long ll;
ll mul(ll a,ll b,ll c)
{
    long long ans=1;
    while(b)
    {
        if(b&1)
        {
            ans=ans*a%c;
        }
        a=a*a%c;
        b>>=1;
    }
    return ans;
}

int main()
{
    int t;
    cin >> t;
    while(t--)
    {
        ll a,b,c;
        cin >> a >> b;
        c=mul(a,b-2,b);
        if(a%b) cout << c << endl;
        else cout << "impossible" << endl;
    }
    return 0;
}

当n不是质数时，可以用扩展欧几里得算法求逆元：
a有逆元的充要条件是a与p互质，所以gcd(a, p) = 1;
假设a的逆元为x，那么有a * x ≡ 1 (mod p)
等价：ax + py = 1
exgcd(a, p, x, y)
模板：

#include <iostream>
using namespace std;
typedef long long LL;
int n;
int exgcd(int a, int b, int &x, int &y)
{
    if (!b) {
        x = 1, y = 0;
        return a;
    }
    int d = exgcd(b, a % b, y, x);
    y -= a / b * x;
    return d;
}

int main()
{
    cin >> n;
    while (n --)
    {
        int a, p, x, y;
        cin >>  a >> p;
        int d = exgcd(a, p, x, y);
        if (d == 1) cout << ((LL)x + p) % p << endl;//将结果变成正数
        else puts("impossible");

    }
    return 0;
}

中国剩余定理

x 1≡m1(%a1);
x2 ≡m2(%a2);
x3 ≡m3(%a3);
：
xk ≡mk(%ak);
先看前两项
k1a1+m1=g2a2+m2;(移项)
k1a1-g2a2=m2-m1;
先判断是否有解
gcd(a1,a2)|m2-m1；
拓展欧几里得得到的是k1 * a1-g2 * a2=gcd(a1,a2)的解k0;
k1=(m2-m1)/gcd(a1，a2) * k0；
这里用gcd (a1,-a2) 也可以，负号放g2里
k1* a1-g2* a2=m2-m1 的通解 k=k1 +t* (a2/gcd(a1,a2));
int dd=a2/gcd(a1,a2);（取绝对值）
在这里要取k1的最小正整数，还可以防止爆范围;
k1=(k1%dd+ dd)%dd (因为t是任意的，相当于放后面了)
将k带入第一个方程得x=a1*k1+m1 + t * dd *a1;
新的 m1=a1 * k1+m1,a1=dd * a1;
同时m1是将最小的k1,带入方程得到的解，是最小的;
模板：

#include <iostream>
using namespace std;
typedef long long LL;
LL x, y;

LL exgcd(LL a, LL b)
{
    if(b == 0)
    {
        x = 1, y = 0;
        return a;
    }
    LL d = exgcd(b, a % b);
    LL z = x;
    x = y;
    y = z - a / b * y;
    return d;
}

int main()
{
    int n;
    cin >> n;
    n--;
    LL a1, a2, m1, m2;
    cin >> a1 >> m1;
    LL d;
    while(n--)
    {
        cin >> a2 >> m2;
        d = exgcd(a1, a2);
        if((m2 - m1) % d)
        {
            cout << "-1" << endl;
            return 0;
        }
        LL k0 = (m2 - m1) / d * x;
        k0 = (k0 % abs(a2 / d) + abs(a2 / d)) % abs(a2 / d);
        m1 = (a1 * k0 + m1);
        a1 = abs(a2 / d * a1);
    }
    cout << m1 << endl;
    return 0;
}

最后放个好看的图片；