斯特林数学习笔记

第一类斯特林数

定义

$s(i,j)$表示把$i$个物品分成$j$个环的方案，或者记做$\left[\begin{array}{c}i\\ j\end{array}\right]$

递推式

考虑第$i$个物品的划分
$s(i,j)=s(i-1,j-1)+s(i-1,j)\ast (i-1)$

枚举第$1$个物品所在环的大小

$s(i,j)={\sum }_{k=1}^n\left(\genfrac{}{}{0ex}{}{n-1}{k-1}\right)(k-1)!s(i-k,j-1)$

快速求$s(n,i)$

考虑构造生成函数$F_n={\sum }_{i=0}^{\infty }s(n,i)x^i$
由递推式可得
$F_n=x{F}_{n-1}+(n-1)F(n-1)=(x+n-1)F(n-1)$

可得$F_n=x^{\overline{n}}$
即${\sum }_{i=0}^{n}s(n,i)x^i=x^{\overline{n}}$
这是就是第一类斯特林数和上升幂的关系

而且也有关于下降幂，由于
$x^{\underline{n}}=(-1{)}^n(-x{)}^{\overline{n}}$

所以${\sum }_{i=0}^n(-1{)}^{n-i}s(n,i)x^i=x^{\underline{n}}=(\genfrac{}{}{0ex}{}{x}{n})n!$

---

对和上升幂的式子
分治$NTT$可以做到$O(nlo{g}^{2}n)$

考虑倍增

如果$n$为奇数，递归求解$n-1$,乘一个$x+n-1$

如果$n$为偶数
有$x^{\overline{2n}}=x^{\overline{n}}(x+n{)}^{\overline{n}}$
假设已经求出了$x^{\overline{n}}=f1(x)={\sum }_{i=0}^n{a}_i{x}^i$
考虑如何求出
$f2(x)={\sum }_{i=0}^n{a}_i(x+n{)}^i$

$f2(x)={\sum }_{i=0}^n{a}_i{\sum }_{j=0}^i{x}^j{n}^{i-j}\left(\genfrac{}{}{0ex}{}{i}{j}\right)$

$={\sum }_{j=0}^n\frac{x^j}{j!}{\sum }_{i=j}^n{a}_{i}i!\frac{n^{i-j}}{(i-j)!}$

$={\sum }_{j=0}^n\frac{x^j}{j!}{\sum }_{i=0}^{n-j}{a}_{i+j}(i+j)!\frac{n^i}{i!}$

令$A(x)={\sum }_{i=0}^n{a}_{i}i!x^i,B(x)={\sum }_{i=0}^n\frac{n^{n-i}}{(n-i)!}{x}^i$

把$(A\ast B)(x)$乘出来左移$n$位即可

复杂度$T(n)=T(\frac{n}{2})+O(nlogn)=O(nlogn)$

代码：

inline poly calc_s(int n){
	poly res;
	if(n==1){res.pb(0),res.pb(1);return res;}
	if(n&1){
		res=calc_s(n-1);
		res.resize(n+1);
		for(int i=n;i;i--)res[i]=add(res[i-1],mul(res[i],n-1));
		return res;
	}
	int mid=n>>1;
	poly a=calc_s(mid),b(mid+1),c(mid+1);
	for(int i=0;i<=mid;i++)c[i]=mul(a[i],fac[i]);
	for(int i=0,p=1;i<=mid;i++,Mul(p,mid))b[i]=mul(p,ifac[i]);
	reverse(b.bg(),b.bg()+mid+1);
	c=c*b;for(int i=0;i<=mid;i++)c[i]=mul(c[i+mid],ifac[i]);
	c.resize(mid+1),res=a*c;return res;
}

例题：传送门

第一类斯特林数快速求行
第一类斯特林数快速求列

第二类斯特林数

定义

$S(i,j)$表示把$i$个物品分成$j$个集合的方案数，也记作$\left\{\begin{array}{c}n\\ m\end{array}\right\}$

递推式

考虑第$i$个数的划分

$S(i,j)=S(i-1,j-1)+j\ast S(i-1,j)$

考虑第一个数所在的集合大小

$S(i,j)={\sum }_{k=1}^i\left(\genfrac{}{}{0ex}{}{i-1}{k-1}\right)S(i-k,j-1)$

快速求解$S(n,x)$

考虑$x^n$是用$x$种颜色去染$n$个格子，枚举最后用了几种颜色，可得

$x^n={\sum }_{i=0}^{min(x,n)}\left(\genfrac{}{}{0ex}{}{x}{i}\right)i!S(n,i)$

实际上由于$\left(\genfrac{}{}{0ex}{}{i}{j}\right)$和$S(i,j)$在$i<j$的时候都为0

所以写作

$x^n={\sum }_{i=0}^n\left(\genfrac{}{}{0ex}{}{x}{i}\right)i!S(n,i)$

和

$x^n={\sum }_{i=0}^x\left(\genfrac{}{}{0ex}{}{x}{i}\right)i!S(n,i)$

都可以
在实际题目中可能会根据不同情况选择变量

对第二个式子二项式反演
得

$x!S(n,x)={\sum }_{i=0}^x(-1{)}^{x-i}(\genfrac{}{}{0ex}{}{x}{i})i^n$

得

$S(n,x)={\sum }_{i=0}^x\frac{(-1{)}^{x-i}}{(x-i)!}\frac{i^n}{i!}$

这是一个卷积的形式，$NTT$可以$O(nlogn)$求出

poly f,g;
for(int i=0;i<=n;i++)f.pb((i&1)?mod-ifac[i]:ifac[i]);
for(int i=0;i<=n;i++)g.pb(mul(ksm(i,n),ifac[i]));
poly S=f*g;

第二类斯特林数快速求行
第二类斯特林数快速求列

补充一个是$s(n,k)=S(-k,-n)$（虽然不知道这个会有什么鬼用）

斯特林反演：

${\sum }_{i=1}^n(-1{)}^{i-1}{S}_{n,i}(i-1)!=[n=1]$
或者写作

若$f(n)={\sum }_{i=1}^{n}S(n,i)g(i)$
$g(n)={\sum }_{i=1}^n(-1{)}^{n-i}s(n,i)f(i)$

看到一个简单清楚的证明：
设$f_n={\sum }_{i=0}^{n}S(n,i)g_i$
那么写成生成函数是
$f(x)=g(e^x-1)$,即$[x^n]f={\sum }_{i=0}^n{g}_i({\sum }_{j=1}^{\infty }\frac{x^j}{j!}{)}^i$
$e^x-1$是集合的$egf$，所以后面乘上$n!$就是$n$分成$i$个环的方案$=S(n,i)$
换元得$f(\ln (x+1))=g(x)$
$[x^n]g={\sum }_{i=0}^n{f}_i({\sum }_{j=1}^{\infty }\frac{x^j(-1{)}^{j-1}}{j}{)}^i$
后面${\sum }_j\frac{x^j}{j}$是环的$egf$，不过有个$(-1{)}^{j-1}$，总共是$n$
可以直接得到总共一定是$(-1{)}^{n-i}$，而环对应就是$s(n,i)$
所以得证