数论基础

Posted 2021-02-13 vscoder

前言：

我的数学真是太菜了，所以这里就都不给出证明了

我是一定不会告诉你其实是因为我不会证的，想看证明的可以看我们大佬的博客

正文：

GCD

$gcd$ 的求法当然要用欧几里得定理，就是辗转相除

int gcd(int a,int b)
{
    if(!b) return a;
    return gcd(b,a%b);
}

求出了 $gcd$ 之后，我们就可以求出 $lcm$（最小公倍数）

有一个性质是 $gcd(a,b) imes lcm(a,b)=a imes b$，这样就可以 $O(1)$ 求出 $lcm$

EXGCD

$exgcd$ 即扩展欧几里得定理

可以用来求解 $ax+by=gcd(a,b)$ 的一组解

具体做法就是辗转相除到底然后回带，具体怎么带请参考代码

以下是魔改董大佬的写法

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

以下是魔改 $pc$ 大佬的写法

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

快速幂

用于指数较大的情况

同时可以一边计算一边取模，防止爆 $long long$

ll qpow(ll a,ll b,ll p)
{
    ll ans=1,base=a%p;//记得先取模
    while(b)
    {
        if(b&1) ans=ans*base%p;
        base=base*base%p;
        b>>=1;
    }
    return ans;
}

龟速乘

魔改的快速幂，防止直接乘起来会爆 $long long$

ll qmul(ll a,ll b,ll p)
{
    ll ans=0,base=a%p;
    while(b)
    {
        if(b&1) ans=(ans+base)%p;
        base=(base+base)%p;
        b>>=1;
    }
    return ans;
}

乘法逆元

如果 $abequiv 1 (mod p)$ 且 $gcd(a,p)=1$ ，则称 $a$ 在 $mod p$ 意义下的乘法逆元为 $b$

我们知道在取模运算下是不能随便除的，所以我们定义了逆元，$dfrac a bequiv a imes inv(b) (mod p)$

化简即可得 $b imes inv(b)equiv 1 (mod p)$

关于它的求法，这里介绍三种

我们知道 $a imes inv(a)equiv 1 (mod p)$ 且 $gcd(a,p)=1$

所以第一种方法就是刚刚介绍的 $EXGCD$

由逆元的定义我们可以得到 $a imes inv(a) +k imes p=1$

而又因为 $gcd(a,p)=1$ ，所以 $a imes inv(a) +p imes k=gcd(a,p)$

于是我们就可以直接解出 $a$ 在 $mod p$ 意义下的逆元

exgcd(a,p,inv,k);
inv=(inv%p+p)%p;//保证逆元为正数

第二种求法是费马小定理

如果 $p$ 为质数，并且 $gcd(a,p)=1$ ，那么 $a^{p-1}equiv 1 (mod p)$

所以它可以用来解决模数 $p$ 是质数的情况

由 $a^{p-1}equiv 1 (mod p)$ 可知 $a imes a^{p-2}equiv 1 (mod p)$

所以 $a$ 在 $mod p$ 意义下的逆元即为 $a^{p-2}$，快速幂可求

第三种求法是欧拉定理

同样的，如果 $gcd(a,p)=1$，则有 $a^{φ(p)}equiv 1 (mod p)$

所以它也可以解决模数 $p$ 不是质数的情况

由 $a^{φ(p)}equiv 1 (mod p)$ 可知 $a imes a^{φ(p)-1}equiv 1 (mod p)$

所以 $a$ 在 $mod p$ 意义下的逆元即为 $a^{φ(p)-1}$，快速幂可求

另外逆元还可以线性递推以及求阶乘逆元，这里我就不写了

后序：

这些模板我仅仅理解了一部分

以后还会陆续补充一些证明及应用

如果有错误的地方欢迎指正

$ps:$ 数学真是有趣呐~~~