正文开始前,感谢以下几篇博客:
https://www.cnblogs.com/y2823774827y/p/10700231.html(不认识的)
https://www.cnblogs.com/p-b-p-b/p/10943453.html(PB的)
https://www.cnblogs.com/AThousandMoons/p/10949178.html(CDW的)
第一类斯特林数的定义:$egin{bmatrix}n\mend{bmatrix}$ 表示把 $n$ 个不同的物品组成恰好 $m$ 个没有区别的圆排列的方案数。
那么有递推公式 $egin{bmatrix}n\mend{bmatrix}=(n-1)egin{bmatrix}n-1\mend{bmatrix}+egin{bmatrix}n-1\m-1end{bmatrix}$。前面是把新的物品塞在一个原来的物品后边,后面是新开一个圆排列。
初始 $egin{bmatrix}1\1end{bmatrix}=1,egin{bmatrix}n\mend{bmatrix}=0(n<m)$。
这毒瘤玩意居然没有单个的通项公式
有个式子:$n^{overline{m}}=sumlimits^{m}_{i=0}egin{bmatrix}m\iend{bmatrix}n^i$。
组合意义?不会
理性证明?看PB的去(算了还是搬过来吧)
$$n^{overline{m+1}}$$
$$(n+m)n^{overline{m}}$$
$$(n+m)sumlimits^{m}_{i=0}egin{bmatrix}m\iend{bmatrix}n^i$$
$$sumlimits^{m+1}_{i=0}(megin{bmatrix}m\iend{bmatrix}+egin{bmatrix}m\i-1end{bmatrix})n^i$$
$$sumlimits^{m+1}_{i=0}egin{bmatrix}m+1\iend{bmatrix}n^i$$
求一行的:$x^{overline{n}}=sumlimits^{n}_{i=0}egin{bmatrix}n\iend{bmatrix}x^i$。
所以只要求出 $x^{overline{n}}$ 这个多项式即可。
设 $F_n(x)=x^{overline{n}}$。
当 $n$ 是奇数时 $F_n(x)=(x+n-1)F_{n-1}(x)$,不管。
当 $n$ 是偶数时 $F_n(x)=F_{frac{n}{2}}(x)F_{frac{n}{2}}(x+frac{n}{2})$。
(为方便将式子变为 $F_{2n}(x)=F_n(x)F_n(x+n)$)
假设已经求出 $F_n(x)=sumlimits_{i=0}^na_ix^i$。
$$F_n(x+n)$$
$$sumlimits_{i=0}^na_i(x+n)^i$$
$$sumlimits_{i=0}^na_isumlimits_{j=0}^idbinom{i}{j}x^jn^{i-j}$$
$$sumlimits_{i=0}^na_ii!sumlimits_{j=0}^idfrac{x^j}{j!}dfrac{n^{i-j}}{(i-j)!}$$
$$sumlimits_{j=0}^ndfrac{x^j}{j!}sumlimits_{i=j}^n(a_ii!)dfrac{n^{i-j}}{(i-j)!}$$
令 $A_i=a_ii!,B_i=dfrac{n^{n-i}}{(n-i)!}$,那么 $x^j$ 的系数是 $dfrac{A imes B}{j!}[x^{j+n}]$。
时间复杂度 $O(nlog n)$。
求一列的:
考虑生成函数 $F(x)=sumlimits_{i=1}^{+infty}(i-1)!dfrac{x^i}{i!}=sumlimits_{i=1}^{+infty}dfrac{x^i}{i}$(注意下界是 $1$)。表示把 $n$ 个球形成圆排列(非空)的指数型生成函数(斯特林数中集合没有区别,所以要指数型)。那么 $dfrac{F^m(x)}{m!}$ 就是把 $n$ 个球形成 $m$ 个圆排列且每个圆排列非空的方案数的指数型生成函数。那么 $dfrac{egin{bmatrix}i\mend{bmatrix}}{i!}=dfrac{F^m(x)}{m!}[x^i]$。$O(nlog n)$。
第二类斯特林数的定义:$egin{Bmatrix}n\mend{Bmatrix}$ 表示把 $n$ 个不同的物品放进 $m$ 个没有区别的集合,且每个集合都非空的方案数。
那么有递推公式 $egin{Bmatrix}n\mend{Bmatrix}=megin{Bmatrix}n-1\mend{Bmatrix}+egin{Bmatrix}n-1\m-1end{Bmatrix}$。前面是把新的物品随便找个集合塞下,后面是新开一个集合。
初始 $egin{Bmatrix}1\1end{Bmatrix}=1,egin{Bmatrix}n\mend{Bmatrix}=0(n<m)$。
有个式子:$n^m=sumlimits^{m}_{i=0}egin{Bmatrix}m\iend{Bmatrix}dbinom{n}{i}i!=sumlimits^m_{i=0}egin{Bmatrix}m\iend{Bmatrix}n^{underline{i}}$。
组合意义理解一下,左边是 $m$ 个不同物品扔进 $n$ 个不同集合,且可以有空集合的方案数。右边枚举有多少个非空集合,然后把 $m$ 个物品放进这些集合里。注意第二类斯特林数中集合没有区别,所以要乘回一个 $i!$ 使他们有区别。
二项式反演后可得通项公式 $egin{Bmatrix}n\mend{Bmatrix}=dfrac{1}{m!}sumlimits^m_{k=0}(-1)^kdbinom{m}{k}(m-k)^n$。
求一行的:
发现 $egin{Bmatrix}n\iend{Bmatrix}=sumlimits^i_{k=0}dfrac{(-1)^k}{k!}dfrac{(i-k)^n}{(i-k)!}$,是个卷积的形式,$O(nlog n)$ 求解。
求一列的:
考虑生成函数 $F(x)=sumlimits_{i=1}^{+infty}dfrac{x^i}{i!}=e^x-1$(注意下界是 $1$)。表示把 $n$ 个球丢入一个集合(非空)的指数型生成函数(斯特林数中集合没有区别,所以要指数型)。类似上面的做法即可。$O(nlog n)$。
还有一种方法。考虑 $F_m(x)=sumlimits_{i=0}^{+infty}egin{Bmatrix}i\mend{Bmatrix}x^i$。那么显然有 $F_0(x)=1$。
$$F_m(x)=sumlimits_{i=0}^{+infty}egin{Bmatrix}i\mend{Bmatrix}x^i$$
$$F_m(x)=sumlimits_{i=0}^{+infty}(megin{Bmatrix}i-1\mend{Bmatrix}+egin{Bmatrix}i-1\m-1end{Bmatrix})x^i$$
$$F_m(x)=msumlimits_{i=0}^{+infty}egin{Bmatrix}i-1\mend{Bmatrix}x^i+sumlimits_{i=0}^{+infty}egin{Bmatrix}i-1\m-1end{Bmatrix}x^i$$
$$F_m(x)=msumlimits_{i=0}^{+infty}egin{Bmatrix}i\mend{Bmatrix}x^{i+1}+sumlimits_{i=0}^{+infty}egin{Bmatrix}i\m-1end{Bmatrix}x^{i+1}$$
$$F_m(x)=mxF_m(x)+xF_{m-1}(x)$$
$$F_m(x)=dfrac{xF_{m-1}(x)}{1-mx}$$
$$F_m(x)=dfrac{x^m}{prodlimits_{i=1}^m(1-ix)}$$
分母实际上是 $prodlimits_{i=1}^m(x-i)$ 的系数倒置。可以用上面的方法。也可以 $O(nlog n)$,常数更小。
首先由两个很显然的式子:
$n^{overline{m}}=(-1)^{n+1}(-n)^{underline{m}}$
$n^{underline{m}}=(-1)^{n+1}(-n)^{overline{m}}$
总结一下上面提到的性质:(顺便引入一些新的)
$n!=sumlimits^n_{i=0}egin{bmatrix}n\iend{bmatrix}$(其实就是下面的式子令 $n=1$,$m=n$,也可以用抽象代数的意义证)
$n^{overline{m}}=sumlimits^m_{i=0}egin{bmatrix}m\iend{bmatrix}n^i$
$n^{underline{m}}=sumlimits^m_{i=0}(-1)^{m-i}egin{bmatrix}m\iend{bmatrix}n^i$(也可以用数学归纳法证)
$n^m=sumlimits^m_{i=0}egin{Bmatrix}m\iend{Bmatrix}n^{underline{i}}=sumlimits^m_{i=0}egin{Bmatrix}m\iend{Bmatrix}dbinom{n}{i}i!$
$n^m=sumlimits^m_{i=0}(-1)^{m-i}egin{Bmatrix}m\iend{Bmatrix}n^{overline{i}}$(下面会讲到怎么证)
$sumlimits^n_{i=0}i^k=frac{(n+1)^{underline{k+1}}}{k+1}-sumlimits_{i=0}^{k-1}(-1)^{k-i}egin{bmatrix}k\iend{bmatrix}sumlimits_{j=0}^nj^i$
$sumlimits^n_{i=0}i^k=sumlimits^k_{i=0}egin{Bmatrix}k\iend{Bmatrix}dfrac{(n+1)^{underline{i+1}}}{i+1}$
关于最后一个式子的证明:(倒数第二个自己来咕了)
$$sumlimits^n_{i=0}sumlimits^k_{j=0}egin{Bmatrix}k\jend{Bmatrix}dbinom{i}{j}j!$$
$$sumlimits^k_{j=0}j!egin{Bmatrix}k\jend{Bmatrix}sumlimits^n_{i=0}dbinom{i}{j}$$
$$sumlimits^k_{j=0}j!egin{Bmatrix}k\jend{Bmatrix}dbinom{n+1}{j+1}$$
$$sumlimits^k_{j=0}j!egin{Bmatrix}k\jend{Bmatrix}dfrac{(n+1)^{underline{j+1}}}{j+1}$$
中间 $sumlimits^n_{i=0}dbinom{i}{j}=dbinom{n+1}{j+1}$ 就是传说中的杨辉恒等式。右边是从 $n+1$ 个位置中选 $j+1$ 个位置,左边是枚举最右边的位置($i+1$),然后从前 $i$ 个位置中选剩下的 $j$ 个。
有些题目中会用到普通幂和上升/下降幂转换的方法来做。根据PB的总结,当普通幂转上升/下降幂时用第二类,当上升/下降幂转普通幂时用第一类。
斯特林反演:
长得和二项式反演很像。
$$f_n=sumlimits^n_{k=0}(-1)^kegin{bmatrix}n\kend{bmatrix}g_kiff g_n=sumlimits^n_{k=0}(-1)^kegin{Bmatrix}n\kend{Bmatrix}f_k$$
$$f_n=sumlimits^n_{k=0}egin{bmatrix}n\kend{bmatrix}g_kiff g_n=sumlimits^n_{k=0}(-1)^{n-k}egin{Bmatrix}n\kend{Bmatrix}f_k$$
用这个式子也可以证明上面的一些式子。
证明的话需要两个反转公式:
$$sum_{i=m}^n(-1)^{n-i}egin{bmatrix}n\iend{bmatrix}egin{Bmatrix}i\mend{Bmatrix}=[n=m]$$
$$sum_{i=m}^n(-1)^{n-i}egin{Bmatrix}n\iend{Bmatrix}egin{bmatrix}i\mend{bmatrix}=[n=m]$$
两个奇怪的柿子:
$$egin{Bmatrix}n+1\m+1end{Bmatrix}=sum_{k}inom{n}{k}egin{Bmatrix}k\mend{Bmatrix}$$
证明:后面那个相当于从 $n$ 个物品中选出一些,并分成 $m$ 个无序非空集合。
前面那个是把 $n+1$ 个物品分成 $m+1$ 个无序非空集合。
钦定这 $n+1$ 个物品中有一个关键物品,分成的 $m+1$ 个集合中恰好一个含有关键物品,删掉这个集合,显然和后半部分等价。
$$egin{bmatrix}n+1\m+1end{bmatrix}=sum_{k}egin{bmatrix}n\kend{bmatrix}inom{k}{m}$$
证明:后面那个相当于把 $n$ 个物品分成若干个圆排列,再选出 $m$ 个圆排列。
前面那个是把 $n+1$ 个物品分成 $m+1$ 个圆排列。
钦定这 $n+1$ 个物品中有一个关键物品,分成的 $m+1$ 个圆排列中恰好一个含有关键物品,将包含关键物品的那个环从关键物品那里断开。
断开之后会变成一个排列,注意到排列是可以和大小之和一定的圆排列集合一一对应的。
将这些圆排列和原来没被断掉的 $m$ 个圆排列合在一起看,显然和后半部分等价。
暂时先咕到这里。