数论基础知识（数论第一波，欧拉筛骚操作）

---

前言

数论中所有的数通常都是正整数，在出现的字母在未经说明的情况下为正整数

模运算

有些时候答案非常大，题目会要求你输出答案取模一个数的值
对于取模，有以下定理
1、$(a+b)\%c=(a\%c+b\%c)\%c$
2、$(a-b)\%c=(a\%c-b\%c)\%c$
3、$(a\ast b)\%c=(a\%c)\ast (b\%c)\%c$
对于符号，我们可以简化成语言
和的余数等于余数的和，差的余数等于余数的差，积的余数等于余数的积。
肯定有人问：“为什么没有商的余数等于余数的商”
对于除法而言，及时除数与被除数都为整数，但商也不一定是整数，而加减乘在$a,b$均为整数的情况下一定是整数

整除

$a\ne 0$，若存在整数$q\ne 0$使得$aq=b$，那么称a是b的因数（约数），b是a的倍数，a能整除b，b能被a整除，记为$a|b$
对于整除，有以下定理
1、$a|b，b|c\Rightarrow a|c$，a是b的因子，b是c的因子，那么a也就是c的因子
2、$a|b，c|d\Rightarrow ac|bd$，a是b的因子，c是d的因子，那么ab也就是cd的因子
3、$a|b，a|c\Rightarrow a|pb+pc，过程：b=na，c=ma，原式等于a|pna+pma，等于a|a(pn+pm)，\because a|a，\therefore a|pb+pc$

质数（prime number，又称素数）

• 质数定义：不含有除了1和自身以外的其他因子的数称为质数（素数）
• 特殊：2是最小的质数，也是唯一的偶质数
• 质数判定定理：若一个数n找不到小于等于$\sqrt{n}$的非1因数，则这个数是质数
• 素数定理：n充分大时，n以内的素数个数约为n/log₂n个
  关于质数的东西有很多呀，不慌，慢慢来

整数的标准分解

一个大于1的整数一定可以唯一地写成若干个质数的幂的积
就像这样 $n=p_1^{w1}\ast p_2^{w2}\ast p_3^{w3}\ast p_4^{w4}\ast \dots \dots \ast p_m^{wm}$

• 互质（coprime）：两个整数除1以外再没有其他因数，那么这两个数互质
• 两个互质的数的唯一分解中没有相同的质数

代码

void factorize(int n){
    for(int i = 2; n != 1; i ++){
        if( n%i == 0 ){
            p[++k] = i;
            while( n%i == 0 ){
                n /= i;
                w[k]++
            }
        }
    }
}

最大公约数与最小公倍数

• 两个数$a,b$的最大公约数记为$gcd(a,b)$或者$(a,b)$（记住共产党）
• 两个数$a,b$的最小公倍数记为$lcm(a,b)$或者$[a,b]$

我们用整数的标准分解来深度了解一下$gcd$和$lcm$
$x=p_1^{a1}\ast p_2^{a2}\ast p_3^{a3}\ast p_4^{a4}\ast \dots \dots \ast p_m^{am}$
$y=p_1^{b1}\ast p_2^{b2}\ast p_3^{b3}\ast p_4^{b4}\ast \dots \dots \ast p_n^{bn}$
其中$p_n$和$p_m$足够大，$a_i，b_i$可以为0
我们知道两个数的$gcd$是两个数都有的最大因子，注意，是最大，说明$p_x$一定取$p_x^{ax}$或$p_x^{bx}$，而两个数都有，说明取$p_x^{min(ax,bx)}$
设$o=max(n,m)$
∴$gcd(x,y)=p_1^{min(a1,b1)}\ast p_2^{min(a2,b2)}\ast p_3^{min(a3,b3)}\ast \dots \dots \ast p_o^{min(ao,bo)}$
而两个数的$lcm$是既有因子$x$，又有因子$y$， 说明$lcm$的$p{1}_x$一定满足$p_x^{max(ax,bx)}|p{1}_x$，但又要求最小，所以$p{1}_x$取$p_x^{max(ax,bx)}$
∴$lcm[x,y]=p_1^{max(a1,b1)}\ast p_2^{max(a2,b2)}\ast p_3^{max(a3,b3)}\ast \dots \dots \ast p_o^{max(ao,bo)}$
而显然$(x,y)\ast [x,y]=x\ast y$

约数个数与约数和

我们从整数的唯一分解入手
$$x=p_1^{a_1}\ast p_2^{a_2}\ast p_3^{a_3}\ast p_4^{a_4}\ast \dots \dots \ast p_m^{a_m}$$
我们可以知道，既然$p_1^{a_1}|x$，说明$p_1^0|x$，$p_1^1|x$，$p_1^2|x$……$p_1^{a_1}|x$，而这里就有$a_1+1$个约数
而既然$p_2^{a_2}|x$，说明就有$a_2+1$个约数，后面以此类推，而前面$a_n+1$约数可以与后面$a_m+1$个约数逐个组合，根据乘法原理，就有$(a_n+1)\ast (a_m+1)$个约数
而根据以上结论，就可以推出$x$的约数个数为$$d[x]=(a_1+1)\ast (a_2+1)\ast (a_3+1)\dots \dots \ast (a_m+1)$$
而我们知道$x$有$p_1^0|x$，$p_1^1|x$，$p_1^2|x$……$p_1^{a_1}|x$这a_1+1个约数，而且对于每一个$x$的质因子都有$a_i+1$个约数，而这$a_1+a_2+\dots \dots a_m+m$个约数可以多个组合，因式分解一下就可得到$$s[x]=(1+p_1^1+p_1^2\dots +p_1^{a_1})\ast (1+p_2^1+p_2^2\dots +p_2^{a_2})\dots \ast (1+p_m^1+p_m^2\dots +p_m^{a_m})$$
我们设$(1+p_i^1+p_i^2\dots +p_i^{a_i})=Q$为$(1)$式，则$p_i\ast Q=(p_i^1+p_i^2\dots +p_i^{a_i+1})$为$(2)$式，
$$\therefore (2)-(1)=p_i\ast Q-Q=(p_i-1)\ast Q=p_i^{a_i+1}-1$$
$$\therefore Q=\frac{p_i^{a_i+1}-1}{p_i-1}$$

欧拉筛法

埃拉托尼斯筛法相对于欧拉筛法而言太不实用，请读者自主学习
欧拉筛法如下：

void sieve(int x){
    for(int i = 2; i <= n; i ++){
        if( !v[i] ){
            p[++cnt] = i;
            v[i] = 1;
        }
        for(int j = 1; j <= cnt&&i*p[j] <= x; j ++){
            v[i*p[j]] = 1;
            if( i % p[i] ==  0)
                break;
        }
    }
}

我相信有很多读者不明白第二个循环，其实无论此时$i$是否为质数，我们都需要把之前筛出来的素数乘$i$来筛掉后面的素数，而欧拉筛的难点就在于对if (i % prime[j] == 0)这步的理解，当i是prime[j]的整数倍时，记 m = i / prime[j]，那么 i * prime[j+1] 就可以变为 (m * prime[j+1]) * prime[j]，这说明 i * prime[j+1] 是 prime[j] 的整数倍，不需要现在筛出，因为在之后筛除过程中i * prime[j+1] 这个合数一定会被prime[j]筛除，prime[j]之后的所有素数同理，所以break跳出循环。

欧拉函数求约数个数d[n]

• 若$n$是素数，约数只有$1,n$，所以$d[n]=2$
• 若$i\%p[j]\ne 0$，$i$不存在$p[j]$这个因子，说明$i\ast p[j]$只包含$p[j]$的一次方，没问题吧，而我们知道约数个数的公式是乘上每个质因子的次数+1，因为$i\ast p[j]$只包含$p[j]$的一次方，$p[j]$的次数也就是一，乘上的也就是$(1+1)=2$了
• 若$i\%p[j]=0$，$i$至少存在$p[j]$的一次方，说明$i\ast p[j]$至少包含$p[j]$的二次方，没问题吧，我们用$num[i]$表示$i$的最小的质因子的次数（因为我们需要通过次数来计算约数个数），则$d[i]$包含了$（1+num[i]）$这个因式，因为$d[i]$肯定也是从之前的某个数乘上$(1+num[i])$得到的，而现在要求$d[i\ast p[j]$，可以肯定的是$d[i\ast p[j]]$肯定存在$num[i]+2$这个因式，而$d[i]$存在$num[i]+1$这个因式，那么如何转换呢，就是先把$num[i]+1$除了之后，再乘上$num[i]+2$，就可以推出$d[i\ast p[j]]$了
  因为$i\ast p[j]$是在$i$的基础上增加了$p$的一次方，所以$num[i\ast p[j]]=num[i]+1$了
  所以$$\begin{gathered} num[i\ast p[j]]=num[i]+1 \\ d[i\ast p[j]]=d[i]/(num[i]+1)\ast (num[i\ast p[j]]+1) \end{gathered}$$

代码

void sieve(int x){
    d[1] = 1;
    for(int i = 2; i <= n; i ++){
        if( !v[i] ){
            p[++cnt] = i;
            v[i] = 1;
            d[i] = 2;
        }
        for(int j = 1; j <= cnt&&i*p[j] <= x; j ++){
            v[i*p[j]] = 1;
            if( i % p[i] ==  0){
                num[i*p[j]] = num[i]+1;
                d[i*p[j]] = d[i]/(num[i]+1)*(num[i]+2);
                break;
            }
            d[i*p[j]] = d[i]*2;
            num[i*p[j]] = num[i]+1;
        }
    }
}

欧拉函数求约数和s[n]

• 若$i$为质数，约数只有$i,1$，所以$s[i]=i+1$
• 若$i\%p[j]\ne 0$，说明$i\ast p[j]$只含有$p[j]$的一次方，没问题吧。而因为我们知道
  求约数和，就是求每个质因子等比数列的乘积，因为$i\ast p[j]$只含有$p[j]$的一次方，所以要在$s[i]$基础上乘的就是$p[j{]}^1$的等比数列，也就是$p[j{]}^1+p[j{]}^0=p[j]+1$
• 若$i\%p[j]=0$，说明$i\ast p[j]$至少含有$p[j]$的二次方，我们设$psum[i]$表示关于$i$的最小质因子$p$的等比数列求和$$1+p+p^2+\dots \dots +p^w$$因为约数和就是这东西的乘积，可以肯定的是，$s[i]=psum[i]$乘上另外质因子的等比数列求和。由于$i\ast p[j]$和$i$均以$p[j]$作为最小的质因子（欧拉筛的性质，自己看代码体会吧），所以$i\ast p[j]$一定包含$psum[i]\ast p+1$这个因式（因为$i\ast p[j]$比$p[j]$多一个$p[j]$的1次方，$psum[i]=1+p+p^2+\dots \dots +p^w$，所以$psum[i\ast p[j]]=1+p+p^2+\dots \dots +p^{w+1}$，所以$psum[i\ast p[j]]$其实就等于$psum[i]\ast p[j]+1$，看的出来吧），那么$s[i\ast p[j]]$到$s[i]$的转换就和$d[i\ast p[j]]$到$d[i]$的转换差不多了
  $$psum[i\ast p[j]]=psum[i]\ast p[j]+1$$
  $$s[i\ast p[j]]=\frac{s[i]}{psum[i]}\ast (psum[i\ast p[j]])$$
  代码

void sieve(int x){
    s[1] = 1;
    for(int i = 2; i <= x; i ++){
        if( !v[i] ){
            p[++cnt] = i;
            v[i] = 1;
            psum[i] = s[i] = i+1;
        }
        for(int j = 1; j <= cnt&&p[j]*i <= x; j ++){
            v[i*p[j]] = 1;
            if( i%p[j] == 0 ){
                psum[i*p[j]] = psum[i]*p[j]+1;
                s[i*p[j]] = s[i]/psum[i]*psum[i*p[j]];
                break;
            }
            s[i*p[j]] = s[i]*(p[j]+1);
        	psum[i*p[j]] = 1+p[j];
        }
    }
}