斯特林数 & 贝尔数学习笔记

背景：

好久没有更$\text{blog}$，最近都在准备模拟（水 ）赛。

第一类斯特林数：

第一类斯特林数第$n$行$m$列，记做${[}_n^m]$。

$\text{que}$：

题目传送门：https://www.luogu.org/problem/P5408。
现在你有$n$个珠子，每一个珠子的颜色各不相同，求这些珠子组成$m$个项链的方案数。
现在求$1,2,3,...m$时方案数。

$\text{sol}$：

设$S_{n,m}$表示上面的方案数。
那么有：
$$S_{n,m}=S_{n-1,m-1}+(n-1)S_{n-1,m}$$

解释一下：对于当前的珠子来说，有两种选择：
$[1]$：新开一个项链，那么方案数就要继承$S_{n-1,m-1}$；
$[2]$：在原来的项链里加入，考虑它会加在某一个珠子的后面，由于之前已经有$(n-1)$个珠子，可以加在任意一个后面，因此方案数为$(n-1)S_{n-1,m}$。

可是这是$\Theta (nm)$的大暴力啊。
第一类斯特林数有一个性质就是它可以表示下面的形式：
$$S_{n,m}=[x^m]x(x+1)(x+2)+...+(x+n-1)$$

现在这个式子就可以分治$\text{NTT}$解决了。
时间复杂度：$\Theta (n{\log }^{2}n)$。
而模板题卡这种做法，因此不能通过。

有$\Theta (n\log n)$的做法，可我太菜了，不想学。

$\text{code}$：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
const LL mod=167772161,G=3,inv_G=55924054;
using namespace std;
	int r[800010],f[800010],a[20][800010],b[20][800010];
	int n,A,B,l,limit;
int dg(int x,int k)
{
	if(!k) return 1;
	int op=dg(x,k>>1);
	if(k&1) return (LL)op*op%mod*x%mod; else return (LL)op*op%mod;
}
int get_inv(int x)
{
	return dg(x,mod-2);
}
void init(int n)
{
	limit=1,l=0;
	while(limit<(n<<1))
		limit<<=1,l++;
	for(int i=1;i<limit;i++)
		r[i]=((r[i>>1]>>1)|((i&1)<<(l-1)));
}
void NTT(int *now,int limit,int op)
{
	for(int i=0;i<limit;i++)
		if(i<r[i]) swap(now[i],now[r[i]]);
	for(int mid=1;mid<limit;mid<<=1)
	{
		int wn=dg(op==1?G:inv_G,(mod-1)/(mid<<1));
		for(int j=0;j<limit;j+=(mid<<1))
		{
			int w=1;
			for(int k=0;k<mid;k++,w=((LL)w*wn)%mod)
			{
				int x=now[j+k],y=(LL)w*now[j+k+mid]%mod;
				now[j+k]=(x+y)%mod;
				now[j+k+mid]=(x-y+mod)%mod;
			}
		}
	}
	if(op==1) return;
	int INV=get_inv(limit);
	for(int i=0;i<limit;i++)
		now[i]=(LL)now[i]*INV%mod;
}
void poly_CDQ(int *f,int dep,int l,int r)
{
	if(l==r)
	{
		f[0]=l;f[1]=1;
		return;
	}
	int mid=(l+r)>>1;
	init(r-l+1+1);
	for(int i=0;i<limit;i++)
		a[dep][i]=b[dep][i]=0;
	poly_CDQ(a[dep],dep+1,l,mid),poly_CDQ(b[dep],dep+1,mid+1,r);
	init(r-l+1+1);
	NTT(a[dep],limit,1),NTT(b[dep],limit,1);
	for(int i=0;i<limit;i++)
		f[i]=(LL)a[dep][i]*b[dep][i]%mod;
	NTT(f,limit,-1);
}
int main()
{
	scanf("%d",&n);
	poly_CDQ(f,0,0,n-1);
	for(int i=0;i<=n;i++)
		printf("%d ",f[i]);
}

第二类斯特林数：

第二类斯特林数第$n$行$m$列，记做$\left\{{}_n^m\right\}$。
这一个更加广泛应用。

$\text{que}$：

题目传送门：https://www.luogu.org/problem/P5395 。
现在你有$n$个珠子，每一个珠子的颜色各不相同，求这些珠子放进$m$个盒子的方案数。
现在求$1,2,3,...m$时方案数。

$\text{sol}$：

设$S_{n,m}$表示上面的方案数。
那么有：
$$S_{n,m}=S_{n-1,m-1}+m{S}_{n-1,m}$$

解释一下：对于当前的珠子来说，有两种选择：
$[1]$：新开一个盒子，那么方案数就要继承$S_{n-1,m-1}$；
$[2]$：在原来的盒子里加入，考虑它会加在某一个盒子里，由于之前已经有$m$个盒子，可以加在任意一个里面，因此方案数为$(n-1)S_{n-1,m}$。

可是这是$nm$的大暴力啊。
第二类斯特林数有一个性质就是它可以表示下面的形式：
$$S_{n,m}=\frac{1}{m!}\sum _{k=0}^m(-1{)}^k{C}_m^k(m-k{)}^n$$

$$S_{n,m}=\frac{1}{m!}\sum _{k=0}^m(-1{)}^k\frac{m!}{k!(m-k)!}(m-k{)}^n$$

$$S_{n,m}=\sum _{k=0}^m(-1{)}^k\frac{1}{k!(m-k)!}(m-k{)}^n$$

$$S_{n,m}=\sum _{k=0}^m\frac{(-1{)}^k}{k!}\frac{(m-k{)}^n}{(m-k)!}$$

为什么呢？
假设每一个盒子都不同，并且我们允许空盒的存在。
但是我们当然就是不允许空盒存在啊。
所以考虑容斥。
枚举我当前有$k$个空盒子。
那么先把这$k$个盒子选出来，也就是$C_m^k$。
然后剩下$m-k$个盒子，$n$个球可以随便放，也就是$(m-k{)}^n$。
考虑到我们盒子是没有区别的，所以除以$m!$。

现在这个式子是一个卷积的形式，就可以$\text{NTT}$解决了。
时间复杂度：$\Theta (n\log n)$。

$\text{code}$：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
const int mod=167772161,G=3,inv_G=55924054;
using namespace std;
	int a[1000010],b[1000010],f[1000010],g[1000010],inv[1000010],Inv[1000010];
	int limit,n,l,r[1000010];
int dg(int x,int k)
{
	if(!k) return 1;
	int op=dg(x,k>>1);
	if(k&1) return (LL)op*op%mod*x%mod; else return (LL)op*op%mod;
}
int get_inv(int x)
{
	return dg(x,mod-2);
}
void init(int n)
{
	limit=1,l=0;
	while(limit<(n<<1))
		limit<<=1,l++;
	for(int i=1;i<limit;i++)
		r[i]=((r[i>>1]>>1)|((i&1)<<(l-1)));
}
void NTT(int *now,int limit,int op)
{
	for(int i=0;i<limit;i++)
		if(i<r[i]) swap(now[i],now[r[i]]);
	for(int mid=1;mid<limit;mid<<=1)
	{
		int wn=dg(op==1?G:inv_G,(mod-1)/(mid<<1));
		for(int j=0;j<limit;j+=(mid<<1))
		{
			int w=1;
			for(int k=0;k<mid;k++,w=((LL)w*wn)%mod)
			{
				int x=now[j+k],y=(LL)w*now[j+k+mid]%mod;
				now[j+k]=(x+y)%mod;
				now[j+k+mid]=(x-y+mod)%mod;
			}
		}
	}
}
void dft(int *f,int n,int limit)
{
	NTT(f,limit,-1);
	int INV=get_inv(limit);
	for(int i=0;i<n;i++)
		f[i]=(LL)f[i]*INV%mod;
}
void Init()
{
	inv[0]=inv[1]=1;
	for(int i=2;i<=n;i++)
		inv[i]=((LL)mod-mod/i)*inv[mod%i]%mod;

	Inv[0]=Inv[1]=1;
	for(int i=2;i<=n;i++)
		Inv[i]=(LL)Inv[i-1]*inv[i]%mod;
}
int main()
{
	scanf("%d",&n);
	Init();
	for(int i=0;i<=n;i++)
	{
		f[i]=(((i&1)?-1ll:1ll)*Inv[i]+mod)%mod;
		g[i]=(LL)dg(i,n)*Inv[i]%mod;
	}
	init(n+1);
	NTT(f,limit,1),NTT(g,limit,1);
	for(int i=0;i<limit;i++)
		f[i]=(LL)f[i]*g[i]%mod;
	dft(f,n+1,limit);
	for(int i=0;i<=n;i++)
		printf("%d ",f[i]);
}

贝尔数：

$\text{que}$：

现在你有$n$个珠子，每一个珠子的颜色各不相同，求这些珠子放进若干个盒子的方案数。
现在求$1,2,3,...m$时方案数。

$\text{sol}$：

枚举盒子数$i$，显然我们有：
$$B_{n,m}=\sum _{i=1}^m\left\{{}_n^i\right\}$$

可是这样的时间复杂度是$\Theta (nm)$的，而且还要预处理第二类斯特林数，极其复杂。
考虑它的递推式：
$$B_n=\sum _{k=1}^{n-1}{C}_{n-1}^k{B}_k$$

理解起来也不难，枚举珠子数$k$，考虑从$n-1$个珠子里选出$k$个珠子放在一个盒子里，那么
未完待续。

后记：

有什么好的套路或性质以后再补充吧，毕竟还没有做题。
欢迎大佬指出错误。
一些性质：
$$n^k=\sum _{i=0}^k\left\{{}_k^i\right\}{C}_n^{i}i!$$

持续更哦。