算法学习笔记【数学知识】

该笔记中主要讲述算法竞赛中常用的数学知识，主要为数论，线性代数，组合数学，概率论，博弈论五个部分。
该笔记中全程用 C++ 语言进行代码讲解。
这里给的例题都会给一道题解，剩下几道练习
制作不易，点个赞吧

• 内容没更新完，敬请期待

质数

定义： 若一个正整数无法被（除了1和它本身）以外的任何自然数除整，则称该数为质数（或素数），否则称该正整数为合数。
素数的分配： 在整个自然数中，质数的数量还是比较少的，对于一个足够大的整数 $N$，不超过 $N$ 的质数大约有 $N/lnN$ 个，即每个 $lnN$ 个数中大约有1个质数。

素数的判定：

试除法：
因为一个素数不能被一个除了1和它本身的数除整，那么我们就直接一个一个除，又因为，一个整数 $n$ 是由一个大于等于 $\sqrt{n}$和一个大于等于$\sqrt{n}$ 的数组成的，那么我们就可以直接进行时间复杂度为$O(\sqrt{n})$ 的计算。

代码实现：

bool isprime(int n)
{
	for(int i = 2; i <= n/i; i++)
	{
		if(n % i== 0) return false;
	}
	return true;
}

素数的筛选：

素数筛（欧拉筛）：
该筛法是算法竞赛中最常用的素数筛法，时间复杂度直接降到了 $O(n)$ 级别，但空间复杂度就会很高，也是O(n)。
欧拉筛的基本思路是：

1. 我们从第一个数字开始找，因为第一个数字2是一个素数，然后就将后面的2的倍数给标记一下，以后将不再遍历，然后再遍历到3 的时候也是一样。
2. 因为遍历的时候

代码实现：

bool used[N];
int primes[N];
void isprime()
{
    int idx = 0;
    for(int i =2 ; i <= N; i++)
    {
        if(!used[i]) primes[++idx] = i;
        for(int j = 1; j <= idx && i * primes[j] <= N; j++)
        {
            used[i * primes[j]] = true;
            if(i % primes[j] == 0) break;
        }
    }
}

质因数的分解：

算术基本定理：

任何一个大于1的正整数都能唯一分解为有限个质数的乘积，可写作：
$$N=p_1^{c_1}{p}_2^{c_2}\cdot \cdot \cdot p_m^{c_m}$$
其中 $c_i$ 都是正整数，$p_i$ 都是质数，且满足 $p_1<p_2<\cdot \cdot \cdot p_m$。
实现方法：
试除法：
依照素数判定中的试除法进行计算，也是遇到一个数进行相除，直到除到除不整为止。
还有一个比这时间复杂度更少的算法，Pollard-Roh，但是这个比较难，这里就不解释了。
下面就直接看代码
代码实现：

int num[N];
int primes[N];
void divide(int x)
{
    int idx = 0;
    for(int i = 2; i <= x/i; i++)
    {
        if(x % i == 0)
            primes[++idx] = i,num[idx] = 0;
        while(x % i == 0)
        {
            x /= i;
            num[idx] ++;
        }
    }
    if(x > 1)
        primes[++idx] = x,num[idx] = 1;
    for(int i = 1; i <= idx; i++)
        cout << primes[i] << " " << num[i] << endl;
}

例题（题解）：

质数距离：
Acwing 196. 质数距离题解

约数

定义：
若一个整数 $n$ 除以整数 $d$ 的余数都为 $0$，即 $d$ 能整除 $n$，则称 $d$ 是 $n$ 的约数，$n$ 是 $d$ 的倍数，记为 $d|n$。
算术基本定理：
在前面的推论里面解释到，一个正整数 $N$ 能被唯一分解为 $N=p_1^{c_1}{p}_2^{c_2}\cdot \cdot \cdot p_m^{c_m}$，其中 $c_i$ 都是正整数，$p_i$ 都是质数，且都满足 $p_1<p_2<p_3<\cdot \cdot \cdot <p_m$，则 $N$ 的正约数集合可写作：
$$\left\{p_1^{b_1}{p}_2^{b_2}\cdot \cdot \cdot p_m^{b_m}\right\}，其中0\le b_i\le c_i$$
所以 $N$ 的正约数的个数为：($\prod$表示连乘)
$$(c_1+1)\times (c_2+1)\times \cdot \cdot \cdot \times (c_m+1)=\prod (c_i+1)$$
那么约数之和为：
$$(1+p_1+p_1^2+\cdot \cdot \cdot +p_1^{c_1})\times \cdot \cdot \cdot \times (1+p_m+p_m^2+\cdot \cdot \cdot +p_m^{c_m})=\prod _{i=1}^m(\sum _{j=0}^{c_i}(p_i^j))$$

求$N$的正约数集合的方法：

1. 试除法：

因为一个数是由小于等于$\sqrt{n}$ 和大于等于$\sqrt{n}$ 的数组成的，那么我们只用遍历前面 $2\sim \sqrt{n}$的数就行了。
推论：
一个整数 $n$ 的约数的个数上界为 $2\sqrt{n}$。
其他和质因数分解的方法一样，这里就不多说了，直接上代码：

int factor[N];
void divide(int x)
{
    int idx = 0;
    for(int i = 2; i <= x/i; i++)
    {
        if(x % i == 0)
        {
            factor[++idx] = i;
            if(i != x / i) factor[++idx] = x/i;
        }
    }
    sort(factor+1,factor + 1+idx);
    for(int i = 1; i <= idx; i++)
        cout << factor[i]  << " ";
}

2. 倍数法：

如果用 “试除法” 求 $1\sim n$中每个数的正约数集合，那么时间复杂度就为$O(N\sqrt{N})$,时间很长，那么反过来考虑的话，我们可以用倍增的方法进行求每个数 $d$ ，$1\sim N$ 中以 $d$ 为约数的数就是 $d$ 的倍数，那么我们就考虑，用倍增的方法进行求倍数，可以减少很多重复的操作。
时间复杂度为：
$O(N+N/2+N/3+\cdot \cdot \cdot +N/N)=O(N\log N)$
代码实现：

void divide(int x)
{
    vector<int> factor[50005];
    for(int i =1 ; i <= x; i++)
        for(int j = 1; j <= x / i; j++)
            factor[i * j].push_back(i);

    for(int i =1 ; i <= x; i++)
    {
        for (int j : factor[i]) {
            cout << j << " ";
        }
        cout << endl;
    }
}

最大公约数：

定义：
最大公约数：

若自然数 $d$ 同时是自然数 $a$ 和 $b$ 的约数，则称 $d$ 是 $a$ 和 $b$ 的公约数，在所有 $a$ 和 $b$ 的公约数中最大的一个则称为 $a$ 和 $b$ 的最大公约数，记为 $gcd(a,b)$。

最小公倍数：

若自然数 $m$ 同时是自然数 $a$ 和 $b$ 的倍数，则称 $m$ 是$a$和$b$的公倍数，在所有 $a$ 和 $b$的公倍数中最小的那个就叫做最小公倍数， 记为 $lcm(a,b)$.

同理，我们可以定义三个数以及更多个数的最大公约数、最小公倍数。
定理：
$$\forall a,b\in N,gcd(a,b)\times lcm(a,b)=a\times b$$

计算最大公约数：

九章算术·更相减损术：

• $\forall a,b\in N,a\ge b,有gcd(a,b)=gcd(b,a-b)=gcd(a,a-b)$
• $\forall a,b\in N,有gcd(2a,2b)=2gcd(a,b)$]

代码实现：

int gcd(int a,int b)
{
    int t = abs(a-b);
    while(t)
    {
        a = b;
        b = t;
        t = abs(a - b);
    }
    return b;
}

欧几里得算法：

$$\forall a,b\in N,b\ne 0,gcd(a,b)=gcd(b,a\mathrm{mod}b)$$
代码实现：

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

注意事项：
对于欧几里得算法的时间复杂度为：$O(\log (a+b))$
欧几里得算法是最常用的求最大公约数的方法，不过对于高精度除法(取模)不容易实现，需要高精度的时候可以考虑用更相减损术代替欧几里得算法。

互质与欧拉函数：

定义：
$$\forall a,b\in N,若gcd(a,b)=1，则称a,b互质$$

对于三个数或者更多个数的情况，我们把 $gcd(a,b,c)=1$ 的情况称为 $a,b,c$ 互质数。把 $gcd(a,b)=gcd(b,c)=gcd(a,c)=1$ 称为两两互质。后者显然是一个更强的条件。

欧拉函数：
定义：
$$1\sim N中与N互质的数的个数被称为欧拉函数，记为\phi (N)。$$
在算术基本定理中，$N=p_1^{c_1}{p}_2^{c_2}\cdot \cdot \cdot p_m^{c_m}$，则：
$$\phi (N)=N\ast \frac{p_1-1}{p_1}\ast \frac{p_2-1}{p_2}\ast \cdot \cdot \cdot \ast \frac{p_m-1}{p_m}=N\ast \prod _{质数p|N}(1-\frac{1}{p})$$
证明：

设 $p$ 是 $N$ 的质因子， $1\sim N$ 中 $p$ 的倍数有 $p,2p,3p,\cdot \cdot \cdot ,(N/p)\ast p$，共 $N/p$ 个。设 $q$ 是 $N$ 的质因子， $1\sim N$ 中 $q$ 的倍数有 $q,2q,3q,\cdot \cdot \cdot ,(N/q)\ast q$，共 $N/q$ 个。如果我们将 $N/p+N/q$ 的个数去掉，那么 $p\ast q$ 的倍数就被排除了两次，需要再加回来。因此，$1\sim N$ 中与 $N$ 含有共同质因子 $p$ 或 $q$ 的个数就为：
$$N-\frac{N}{p}-\frac{N}{q}+\frac{N}{pq}=N\ast (1-\frac{1}{p}-\frac{1}{q}+\frac{1}{pq})=N\ast (1-\frac{1}{p})(1-\frac{1}{q})$$

实际上，上述思想被称为容斥原理，类似的，可以再 $N$ 的全部质因子上使用容斥原理，即可得到像上述的公式。

根据欧拉函数的计算式，我们只需要分解质因数，即可顺便求出欧拉函数。
代码实现：
试除法求一个整数$N$ 的欧拉函数值：

int phi(int n)
{
    int ans = n;
    for(int i = 2; i <= n/i; i++)
    {
        if(n % i == 0)
        {
            ans = ans/ i * (i - 1);
            while(n % i == 0)
                n /= i;
        }
    }
    if(n > 1) ans = ans / n * (n - 1);
    return ans;
}

也可以用线性筛的方法求 $1\sim N$之间的每个数的欧拉函数：
讲解写在另一个博客中了，不懂的可以去翻看
线性筛求欧拉函数

void phi(int x)
{
    cnt = 0;
    phis[1] = 1;
    for(int i = 2; i <= x; i++)
    {
        if(!used[i])
        {
            phis[i] = i-1;
            primes[++cnt] = i;
        }
        for(int j = 1; i * primes[j] <= x; j++)
        {
            used[i * primes[j]] = true;
            if(i % primes[j] == 0)
            {
                phis[i * primes[j]] = primes[j] * phis[i];
                break;
            }
            phis[i * primes[j]] = (primes[j] - 1) * phis[i];
        }
    }
}

性质：

1. $\forall n>1，1\sim n$ 中与 $n$ 互质的数的和为 $n\ast \phi (n)/2$。

证明：因为 $gcd(n,x)=gcd(n,n-x)$，所以与 $n$ 不互质的数 $x,n-x$都是成对存在的，平均值为 n/2。因此，与 $n$ 互质的数的平均值也是 $n/2$，进而得出性质1。

2. 若 $a,b$ 互质，则 $\phi (ab)=\phi (a)\phi (b)$。

证明：根据欧拉函数的计算式，对 $a$，$b$ 分解质因数，直接可得到性质2，把性质2推广到一般的函数中，就可以得到 积性函数

3. 积性函数：如果当 $a,b$ 互质时，有$f(ab)=f(a)\ast f(b)$,那么称函数 $f$ 为积性函数。
4. 若 $f$ 是积性函数，且在算术基本定理中 $n={\prod }_{i=1}^m{p}_i^{c_i}$，则 $f(n)={\prod }_{i=1}^{m}f(p_i^{c_i})$。

证明：把 $n$ 分解质因子，按照积性函数的定义，性质4显然是正确的。

5. 设 $p$ 为质数，若 $p|n$ 且 $p^2\mid n$，则 $\phi (n)=\phi (n/p)\ast p$。

证明：若 $p\mid n$ 且 $p^2\mid n$，则 $n,n/p$ 包含相同的质因子，只是 $p$ 的指数不同。直接按照欧拉函数将 $\phi (n)$ 和 $\phi (n/p)$ 计算出来，然后进行相除，得到的就是 p。

6. 设 $p$ 为质数，若 $p|n$ 且 $p^2\nmid n$，则 $\phi (n)=\phi (n/p)\ast (p-1)$。

证明：若 $p\mid n$ 且 $p^2\nmid n$，说明 $p,n/p$ 互质，因为$\phi$是积性函数，所以 $\phi (n)=\phi (n/p)\ast \phi (n)$，而 $\phi (p)=p-1$。

7. ${\sum }_{d\mid n}\phi (d)=n$

证明：设 $f(n)={\sum }_{d\mid n}\phi (d)$。用乘法分配律展开比较，再利用 $\phi$ 是积性函数，得到：若 $n,m$ 互质，那么 $f(nm)={\sum }_{d\mid n}\phi (d)=({\sum }_{d\mid n}\phi (d))\ast ({\sum }_{d\mid n}\phi (n))=f(n)\ast f(m)$。即 ${\sum }_{d\mid n}\phi (d)$ 是积性函数。对于单个质因子，$f(p^m)={\sum }_{d\mid p^m}\phi (d)=\phi (1)+\phi (p)+\phi (p^2)+\cdot \cdot \cdot +\phi (p^m)$ 是一个等比数列求和数列求和再加1，结果为 $p^m$。所以 $f(n)={\sum }_{i=1}^{m}f(p_i^{c_i})={\prod }_{i=1}^m{p}_i^{c_i}=n$。

例题（题解）：

反素数(约数)：
AcWing 198. 反素数题解
Hankson的趣味题(最大公约数)：
Acwing 200. Hankson的趣味题题解
可见的点（互质与欧拉函数）：
Acwing 201. 可见的点题解
扩展
蓝桥杯-3522 互质数的个数

同余

定义：
  若整数 $a$ 和整数 $b$ 除以正整数 $m$的余数相等，则称 $a,b$ 模 $m$ 同余，则记为 $a\equiv b(\mathrm{mod}m)$。
同余与剩余系：
  对于 $\forall a\in [0,m-1]$，集合$a+km(k\in Z)$ 的所有数模 $m$ 同余，余数都是 $a$。该集合称为一个模 $m$ 的同余类，简称为 $\bar{a}$
  模 $m$ 的同余类一共有 $m$ 个，分别为 $\bar{0},\bar{1},\bar{2}...\overline{m-1}$。它们构成 $m$ 的完全剩余系。
  模 $1\sim m$ 中与$m$ 互质的数代表的同余类共有 $\phi (m)$个，它们构成 $m$ 的简化剩余系，例如：模 10 的简化剩余系为 { 1 ˉ , 3 ˉ , 7 ˉ , 9 ˉ } \left\{\bar{1},\bar{3},\bar{7},\bar{9}\right\} {1ˉ,3ˉ,7ˉ,9ˉ}。
  简化剩余系关于模 $m$ 乘法封闭(一个集合中任意两个元素相乘还是在这个集合中)。这是因为若 $a,b(1\le a,b\le m)$ 与 $m$ 互质，则 $a\ast b$也不可能与 $m$ 含有相同的质因子，即 $a\ast b$ 也与 $m$ 互质。再由余数的定义即可得到 $a\ast b\mathrm{mod}m$ 也与 $m$ 互质。既 $a\ast b\mathrm{mod}m$ 也属于 $m$ 的简化剩余系。
费马小定理：

  若 $p$ 是质数，则对于任何整数 $a$，有 $a^p\equiv 1(\mathrm{mod}n)$，其中 $\phi (n)$为欧拉函数。

欧拉定理：

  若正整数 $a,n$ 互质，则 $a^{\phi (n)}=1(\mathrm{mod}n)$，其中 $\phi (n)$ 为欧拉函数。

证明：

  设 $n$ 的简化剩余系为 { a 1 ‾ , a 2 ‾ . . . , a φ n ‾ } \left\{\overline{a_1},\overline{a_2}...,\overline{a_{\varphi{n}}}\right\} {a1​​,a2​​...,aφn​​}。对 $\forall a_i,a_j$，若$a\ast a_i\equiv a\ast a_j(\mathrm{mod}n)$，则 $a\ast (a_i-a_j)\equiv 0$。因为 $a,n$ 互为质数，所以 $a_i-a_j\equiv 0$，即 $a_i\equiv a_j$。故当 $a_i\ne a_j$ 时，$a{a}_i,a{a}_j$ 也代表不同的同余类。
  又因为简化剩余系关于模 $n$ 乘法封闭，故 $\overline{a{a}_i}$也在简化剩余系集合中，因此，集合 $\left\{\overline{a_1},\overline{a_2},...,\overline{a_{\phi (n)}}\right\}$ 与集合 $\left\{\overline{a{a}_1},\overline{a{a}_2},...,\overline{a{a}_{\phi (n)}}\right\}$都能表示 $n$ 的简化剩余系。综上所诉：
$$a^{\phi (n)}{a}_1{a}_2...a_{\phi n}\equiv (a{a}_1)(a{a}_2)...(a{a}_{\phi (n)})\equiv a_1{a}_2...a_{\phi (n)}(\mathrm{mod}n)$$

  因此 $a^{\phi (n)}\equiv 1(\mathrm{mod}n)$。
  当 $p$ 是质数时，$\phi (p)=p-1$，并且只有 $p$ 的倍数与 $p$ 不互质。所以，只要 $a$ 不是 $p$ 的倍数，就有 $a^{p-1}\equiv 1(\mathrm{mod}n)$，两边同乘 $a$ 就是费马小定理。另外，若 $a$ 是 $p$ 的倍数，费马小定理显然成立。
  证毕。
证明：

  设 $b=q\ast \phi (n)+r$，其中 $0\le r<\phi (n)$，即 $r=b\mathrm{mod}\phi (n)$。于是：
$$a^b\equiv a^{p\ast \phi (n)+r}\equiv (a^{\phi (n)}{)}^q\ast a^r\equiv 1^q\ast a^r\equiv a^r\equiv a^{b\mathrm{mod}\phi (n)}(\mathrm{mod}n)$$

证毕。
   许多计数类的题目要求我们把答案对一个质数 $p$ 取模后输出。面对 $a+b,a-b,a\ast b$ 这样的算式，可以在计算前把 $a,b$对 $p$取模，面对乘方算式，根据欧拉定理的推论，可以先把底数对p取模、指数对$\phi (p)$取模，再计算乘方。

• 当$a,n$不一定互质且$b>\phi (n)时，有a^b\equiv a^{bmod\phi (n)+\phi (n)}modn$

矩阵乘法

组合计数

概率和数学期望

0/1分数规划

博弈论与SG函数