码蹄集盒子与球 第二类Stirling数（斯特林数）

n个不同的球放入r个不同的盒子+盒子不能为空，很容易联想到第二类Stirling数（n个不同的球放入r个相同的盒子+盒子不能为空）。

第二类Stirling数：

S(n,m)--将n个不同的球放到m个相同的盒子的方案数

S(n,m)可以由两个部分推出：

• S(n-1,m-1)，此时也就是将第n个球放到第m个盒子中
• S(n-1,m)*m，此时前n-1个球放到m个盒子中，那么第n个球可以有n种放法（在放之前，每个盒子中的球都不相同，所以每个盒子都不相同，就有n种放法）。

即S(n,m)=S(n-1,m-1)+S(n-1,m)*m

初始化：（右面的矩阵g[i][j]代表将i个不同的小球放到j个相同的盒子中，盒子不能为空）

代码如下：

long long int g[20][20];//i个球放入j个盒子
int n,r;

void initial(){
    int min_temp=min(r,n);
    for(int i=0;i<=min_temp;i++){
        g[i][i]=1;
    }
    for(int i=1;i<=n;i++){
        g[i][0]=0;
    }
}

void update(){
    for(int i=1;i<=n;i++){
        for(int j=1;j<i;j++){
            g[i][j]=g[i-1][j-1]+g[i-1][j]*j;
        }
    }
}

将n个不同的小球放入r个不同的盒子中只需要将盒子的摆放考虑进去，也就是S(n,r)∗r!

将n个不同的小球放入r个相同的盒子中，允许为空：

将n个不同的小球放入r个不同的盒子中，允许为空：

参考博客：第二类Stirling数（第二类斯特林数）-CSDN博客

本题代码如下：

#include<bits/stdc++.h> 

using namespace std;
long long int g[20][20];//i个球放入j个盒子
int n,r;

void initial(){
    int min_temp=min(r,n);
    for(int i=0;i<=min_temp;i++){
        g[i][i]=1;
    }
    for(int i=1;i<=n;i++){
        g[i][0]=0;
    }
}

void update(){
    for(int i=1;i<=n;i++){
        for(int j=1;j<i;j++){
            g[i][j]=g[i-1][j-1]+g[i-1][j]*j;
        }
    }
}

long long int multiply(int num){
    long long int res=1;
    for(int i=2;i<=num;i++){
        res*=i;
    }
    return res;
}

int main( )
{
    cin>>n>>r;
    //n个球放入r个相同的盒子
    initial();
    update();
    //n个球放入r个不同的盒子
    cout<<g[n][r]*multiply(r);
    return 0;
}