算法竞赛入门经典计数基础

有重复元素的全排列（可重集）

问题描述：有k个元素，其中第i个元素有n_i个，求全排列个数。

这个问题，在组合数学中，就是可重集的全排列。

可重集在组合数学上表示为S={a₁n₁, a₂n₂, a₃n₃…a_kn_k}

则S的全排列 = $$\frac{n!}{n_1!n_2!n_3!...n_k!}$$

结论推导过程

推导过程：一共有n个位置
首先考虑第a₁号元素，有n₁个。
然后所以有C(n, n₁)
然后考虑第a₂号元素，有n₂个。
所以由C(n - n₁, n₂)
由乘法原理
设X为S的全排列个数
则由乘法原理
X = C(n, n₁)C(n - n₁, n₂)…C(n - n₁ - n₂…n_k-1, n_k)
又知
C(n, k) = p(n, k) / k! = n! / n!(m - n)!
代入X等式的右边
得
$$\frac{n!}{n_1!(n-n_1)!}\frac{(n-n_1)!}{n_1!(n-n_1-n_2)!}\dots \frac{(n-n_1-n_2-\cdots -n_k)!}{n_k!(n-n_1-n_2-\dots n_k)!}$$
约去公因式
即得$$X=\frac{n!}{n_1!n_2!n_3!...n_k!}$$

可重复选择的组合

将求的组合数转换成方程的解，这个转换就很有灵性了。
问题描述：有n个不同元素，每个元素元素可以选多次，选k个元素，有多少种方法

方程1 x₁ + x₂ + x₃ + … + x_n = k
方程2 y₁ + y₂ + y₃ + … + y_n = k + n
其中y_i = x_i + 1
为什么方程1要转换成2
假设不转换成2，那么存在零解，对于k个1，并不知道取成多少份。
转换成了方程2，每一个y_i都必须大于1.。那么k + n必须要分成n个部分，就是在k + n - 1个候选分割线选n - 1,即C(k + n - 1, n - 1) = C(n + k - 1, k).

组合数推导（杨辉三角）

$$(a+b{)}^n=\sum C_n^i{a}^{n-k}{b}^k$$
组合数推导
$$C_n^i=\frac{n-k+1}{k}{C}_n^{n-i}$$

c[0] = 1;
for(int i = 1; i <= n; ++i) 
	c[i] = c[i - 1] * (n - i + 1) / i;\

重头戏

约数的个数

Q:给出正整数n的唯一分解式$$n=p_1^{a_1}{p}_2^{a_2}\dots p_k^{a_k}$$,求n的正约数的个数。

n的正约数个数为：
$$\prod (a_i+1)=(a_1+1)(a_2+1)(a_3+1)\dots (a_k+1)$$
个人对这条结论的简易理解：
无论p_i中的a_i取多少，都是n的因数。
对于某个p_i，可以取0，1，2…a₁
那根据乘法原理
即得上式。

小于n且与n互素的整数个数

$$\psi (x)=\sum (-1{)}^{|s|}\frac{n}{\prod p_i}$$
为什么可以用容斥原理可以得到上面那条公式呢
可以知道，对于唯一分解式$$n=p_1^{a_1}{p}_2^{a_2}\dots p_k^{a_k}$$其中p_i必然都是n的素因子，那首先第一步就是加上n为因子的$$-\frac{n}{p_1}$$仔细体会一下这个式子就可以了。
那根据容斥，必然有多减的了，所以又要加上$$\frac{n}{p_1{p}_2}$$
依次加$$(-1{)}^{|s|}\frac{n}{p_1{p}_2\dots p_k}$$,|s|指的是分母k的大小

下一步式子我自己顺着推不出了，但是逆着推是可以理解的

$$\psi (n)=n(1-\frac{1}{p_1})(1-\frac{1}{p_2})\dots (1-\frac{1}{p_k})$$
就是每个括号都选一个，乘起来恰好就是第一条式子。

最重要的就是做代码解释的笔记了

代码解释

假设不给出唯一分解式，就要找出素因子。
至于为什么只要找到$$\sqrt{n}$$
发现对于唯一分解式$$n=p_1^{a_1}{p}_2^{a_2}\dots p_k^{a_k}$$
如果n不为质数，那么$$p_k<=\sqrt{n}$$
如果n为质数，那没得好说。$$p_k=n$$
上述好像是对的，刘汝佳的代码也只是找到了根号n，书上也说要依法判断根号n的就行了。
$$ans=ans/i\ast (i-1)$$这个式子是$$1-\frac{1}{p_1}$$的变形
“除干净”是
while(n % i == 0) n /= i;
这个操作是除以 i，将因数是i的合数同时自己本身又是n的合数的因数剔除。下一个n % i == 0，i就肯定是素因子了。

int euler_phi(int n) {
	int m = sqrt(n + 0.5);
	for(int i = 2; i <= m;++i) {
		if(n % i == 0) {
			ans = ans /i * (i - 1);
			while(n % i == 0) n /= i;
		}
	}
	if(n > 1) ans = ans / n * (n -  1);  //这条语句就是如果n是质数的话，那么for循环不会对结果有作用
}		

1~n的所有欧拉函数值

void phi_table(int n, int *phi) {
	for(int i = 2; i <= n; ++i) p[i] = 0;
	p[1] = 1;
	for(int i = 2; i <= n; ++i) {//从2开始，2是素数
		if(!phi[i]) {  //这句话就判定了i必然是素数，因为嵌套循环里将 <i 的合数都筛了
			for(int j = i; j <= n; j += i)  {
				if(!phi[j]) phi[j] = j;//将因子有i的都筛了
				phi[j] = phi[j] / i * (i - 1);  //i必然是素数，所以可以相除
			}
		}
	}