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;
}