一只蒟蒻的数论学习笔记（约数篇二）

更多请见DUMBLOG

一只蒟蒻的数论学习笔记（约数篇二）

欧拉函数

定义

1 到 $N$ 中与 $N$ 互质的数的个数被称为欧拉函数，记为 $\phi (N)$。

性质

显然，如果 $N$ 为质数，则 $\phi (N)=N-1$

欧拉函数为 积性函数： 若 $\gcd (a,b)=1$，则 $\phi (ab)=\phi (a)\phi (b)$。

若 $n=p^k$，$p$ 为质数，则 $\phi (n)=n(1-\frac{1}{p})$。

证明：
若 $x$ 与 $p^k$ 不互质，则有 $p|x$ 成立。
$x$ 共有 $\frac{n}{p}$ 个，因此 $\phi (n)=n-\frac{n}{p}=n(1-\frac{1}{p})$

若 $N=p_1^{c_1}{p}_2^{c_2}...p_m^{c_m}$，则：
$$\phi (N)=N\times \frac{p_1-1}{p_1}\times \frac{p_2-1}{p_2}\times ...\times \frac{p_m-1}{p_m}$$
$$\phi (N)=N\times \left(1-\frac{1}{p_1}\right)\times \left(1-\frac{1}{p_2}\right)\times ...\times \left(1-\frac{1}{p_m}\right)$$
即：
$$\phi \left(n\right)=n\times \prod _{i=1}^m\left(1-\frac{1}{p_i}\right)$$

证明：$\phi$ 是积性函数。把 $n$ 拆成 $p_i^{a_i}$ 的乘积形式可立即得证。

代码

我们只要在分解质因数时顺便将欧拉函数求出来就好了。

int phi(int n) 
{
	int ans = n;
	for (int i = 2; i <= sqrt(n); 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;
}

同余

定义

若 $a\mathrm{mod}m=b\mathrm{mod}m$ ，则称 $a$ 与 $b$ 同余，写作 $a\equiv b(\mathrm{mod}m)$

同余类、剩余系

对于一个集合 { x ∣ x = a + k m , a ∈ [ 0 , m − 1 ] , k ∈ Z } \{ x|x=a+km,a \in \left [ 0,m-1 \right ] , k \in Z \} {x∣x=a+km,a∈[0,m−1],k∈Z} 的所有数模 $m$ 同余，余数都是 $a$ 。该集合被称为一个模 $m$ 的同余类，简记为 ā。
模 $m$ 的同余类一共 $m$ 个，它们构成 $m$ 的完全剩余系。
1 到 $m$ 中与 $m$ 互质的数代表的同余类共有 $\phi (m)$ 个，它们构成 $m$ 的简化剩余系

欧拉定理

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

证明

设 S = { a 1 , a 2 , . . . , a φ ( n ) } S=\{ a_{1},a_{2},...,a_{\varphi(n)} \} S={a1​,a2​,...,aφ(n)​} 为 $n$ 的简化剩余系。
设 $a$ 与 $n$ 互质。
$\because a\equiv 1(\mathrm{mod}n)$，$a{a}_1\equiv 1(\mathrm{mod}n)$
$\therefore a{a}_1\equiv 1(\mathrm{mod}n)$
同理，$a{a}_2,a{a}_3,...,a{a}_{\phi (n)}\equiv 1(\mathrm{mod}n)$
$\therefore a{a}_1\times a{a}_2\times ...\times a{a}_{\phi (n)}\equiv 1(\mathrm{mod}n)$
$\therefore a^{\phi (n)}\left(a_1{a}_2...a_{\phi (n)}\right)\equiv a_1{a}_2...a_{\phi (n)}(\mathrm{mod}n)$
$\therefore a^{\phi (n)}\equiv 1(\mathrm{mod}n)$

费马小定理

若 $p$ 为质数，则对于任意整数 $a$，有 $a^p\equiv a(\mathrm{mod}p)$，即 $a^{p-1}\equiv 1(\mathrm{mod}p)$
证明太简单了，懒得证（用欧拉函数和欧拉定理证）

扩展欧几里得算法

贝祖定理（裴蜀定理）

定义

对于任意整数 $a$，$b$，存在一对整数 $x$，$y$，满足 $\mathbf{ax}+\mathbf{by}=\gcd (\mathbf{a},\mathbf{b})$。特别地，一定存在 $x$，$y$ 满足 $ax+by=d$。
等价的表述：不定方程 $ax+by=c$ （$a$，$b$，$c$ 为整数）有解的充要条件为 $\gcd (a,b)|c$
推论：$a$，$b$ 互质等价于 $ax+by=1$ 有解

原理

考虑如何求得 $ax+by=d$ 的一个解。这里 $d=\gcd (a,b)$

考虑使用欧几里德算法的思想，令 $a=bq+r$，其中 $r=a\mathrm{mod}b$；递归求出 $bx+ry=d$ 的一个解。

设求出 $bx+ry=d$ 的一个解为 $x=x_0$，$y=y_0$，考虑如何把它变形成 $ax+by=d$ 的解。

将 $a=bq+r$ 代入 $ax+by=d$、化简得 $b(xq+y)+rx=d$

我们令 $xq+y=x_0$，$x=y_0$，则上式成立

故 $x=y_0$，$y=x_0-y_{0}q$ 为 $ax+by=d$ 的解

边界情况：$b=0$ 时，令 $x=1$，$y=0$

代码

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;
}

求出 ax + by = c 的所有解

先用 exgcd 求出任意一个解 $x_0$ 和 $y_0$。
则 $a{x}_0+b{y}_0=c$
再求出 $ax+by=0$ 的最小（绝对值）解 $x^{\prime }=\frac{b}{\gcd (a,b)}$ 和 $y^{\prime }=\frac{-a}{\gcd (a,b)}$
则 $(a{x}_0+b{y}_0)+k(a{x}^{\prime }+b{y}^{\prime })=c$
即 $a(x_0+k{x}^{\prime })+b(y_0+k{y}^{\prime })=c$
所有解就是 $x=x_0+k{x}^{\prime }$，$y=y_0+k{y}^{\prime }$，$k\in Z$

乘法逆元

若整数 b，m 互质，并且 $b|a$，则存在一个整数 $x$，使得 $\frac{\mathbf{a}}{\mathbf{b}}\equiv \mathbf{a}\times \mathbf{x}(\mathrm{mod}\mathbf{m})$。
称 $x$ 为 $b$ 的模 $m$ 乘法逆元，记为 $b^{-1}(\mathrm{mod}m)$。
可以推得：当 $m$ 为质数时，$b^{m-2}$ 即为 $b$ 的乘法逆元（即 $x$ ）。
那么咋求逆元捏？

快速幂求逆元

快速幂

快速幂是啥？
可以看成是二进制优化的用来求 ${\mathbf{a}}^{\mathbf{b}}\mathrm{mod}\mathbf{p}$ 的东西。

代码

求逆元相当于求 $b^{p-2}\mathrm{mod}p$，所以能用快速幂求。

#include<iostream>
#include<cstdio>
#include<algorithm>
#define int unsigned long long
using namespace std;
const int N = 3000101;
int n, m;
inline int qm(int a, int k, int p) 
{
	int ans = 1;
	while (k) 
	{
		if (k & 1)ans = ans * a % p;
		k >>= 1;
		a = a * a % p;
	}
	return ans;
}
signed main() {
	scanf("%lld%lld", &n, &m);
	for (int i = 1; i <= n; i++) printf("%lld\n", qm(i, m - 2, m));
	return 0;
}

扩展欧几里得算法（贝祖定理）求逆元

求逆元又等价于求 $ax+by=1$
所以可以直接上板子。

代码

int inv(int a, int b) 
{
	int x, y;
	exgcd(a, b, x, y);
	return x;
}

线性求逆元

在求 $i$ 的逆元时，

$\because p=\lfloor \frac{p}{i}\rfloor \times i+p\mathrm{mod}i$（⌊ ⌋为向下取整），

令 $a=\lfloor \frac{p}{i}\rfloor$，$b=p\mathrm{mod}i$，

$\therefore p=a\times i+b$

$\because a\times i+b=p$

$\therefore a\times i+b=0(\mathrm{mod}p)$

$\therefore a\times i=-b(\mathrm{mod}p)$

$\therefore i=-\frac{b}{a}(\mathrm{mod}p)$

$\therefore i^{-1}=-\frac{a}{b}(\mathrm{mod}p)$

综上，$i$ 的逆元为 $\frac{-\lfloor \frac{p}{i}\rfloor \times i}{p\mathrm{mod}i}$

设 $inv[i]$ 为 $i$ 的逆元，

$\because i^{-1}=-\frac{a}{b}(\mathrm{mod}p)$

$\therefore inv[i]=-a\times b^{-1}(\mathrm{mod}p)$

$\because b^{-1}=inv[b]$

$\therefore inv[i]=-\frac{a}{inv[b]}(\mathrm{mod}p)$

$\therefore inv[i]=-\frac{\lfloor \frac{p}{i}\rfloor }{inv[p\mathrm{mod}i]}(\mathrm{mod}p)$

$\therefore inv[i]=\frac{p-\lfloor \frac{p}{i}\rfloor }{inv[p\mathrm{mod}i]\mathrm{mod}p}(\because x\mathrm{mod}p=(x+p)\mathrm{mod}p)$

有了公式后就能用线性的复杂度求出每个数的逆元了。

代码

for (inv[1] = 1, i = 2; i <= n; ++i)
	inv[i] = (p - p / i) * inv[p % i] % p;

线性同余方程

形如 $ax\equiv c(\mathrm{mod}m)$ 的方程（关于 $x$），称为线性同余方程。
该方程等价于 $ax+my=c$，有解条件为 $\gcd (a,m)|c$ 。
所以可以用 拓展$\gcd$ 来求。
求出 $x_0$ 和 $y_0$，满足 $a{x}_0+m{y}_0=\gcd (a,m)$。
然后 $x=x_0\times \frac{b}{\gcd (a,m)}$
$x$ 即为解。




剩下的下次再说吧。