斯特林数包括第一类斯特林数和第二类斯特林数
第一类斯特林数:
定义:n个不同的元素构成m个圆排列(循环排列)的方法数
分成两种情况:
1· 前n-1个元素已经构成了m个圆排列,只要将第n个元素插入1~n-1个元素中任意一个元素的左边即可((n-1)*s(n-1,m))
2· 前n-1个元素已经构成了m-1个圆排列,那么第n个元素只能单独构成一个圆排列才能构成m个圆排列(s(n-1,m-1))
递推式: s(n,m)=s(n-1,m-1)+(n-1)*s(n-1,m);
long long s[MAXN][MAXN];//第一类Stirling数
#define MAX 1e9+7
void init()//预处理
{
memset(s,0,sizeof(s));
s[1][1]=1;
for(int i=2;i<=MAXN-1;i++)
for(int j=1;j<=i;j++)
{
s[i][j]=s[i-1][j-1]+(i-1)*s[i-1][j];
if(s[i][j]>=MAX)
s[i][j]%=MAX;
}
}
第二类斯特林数:
定义:n个不同的元素拆分成m个集合(不考虑次序)的方法数
分成两种情况:
1·前n-1个元素已经拆分成了m个集合,第n个元素和任意一个集合合并(m*s(n-1,m))
2·前n-1个元素拆分成了m-1个集合,第n个元素单独成为一个集合(s(n-1,m-1))
递推式: s(n,m)=m*s(n-1,m)+s(n-1,m-1);
long long s[MAXN][MAXN];//第一类斯特林数
#define MAX 1e9+7//取模
void init()//预处理
{
memset(s,0,sizeof(s));
s[1][1]=1;
for(int i=2;i<=MAXN-1;i++)
for(int j=1;j<=i;j++)
{
s[i][j]=s[i-1][j-1]+j*s[i-1][j];
if(s[i][j]>=MAX)
s[i][j]%=MAX;
}
}