【组合数学】斯特林数

斯特林数包括第一类斯特林数和第二类斯特林数

第一类斯特林数：

定义：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;  
    }  
}