题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2512
题目大意:
因为长期钻研算法, 无暇顾及个人问题,BUAA ACM/ICPC 训练小组的帅哥们大部分都是单身。某天,他们在机房商量一个绝妙的计划"一卡通大冒险"。这个计划是由wf最先提出来的,计划的内容是,把自己的联系方式写在校园一卡通的背面,然后故意将自己的卡"遗失"在某处(如水房,TD,食堂,主M。。。。)他们希望能有MM看到他们遗失卡,能主动跟他们联系,这样就有机会请MM吃饭了。他们决定将自己的一卡通夹在基本相同的书里,然后再将书遗失到校园的各个角落。正当大家为这个绝妙的计划叫好时,大家想到一个问题。很明显,如果只有一张一卡通,那么只有一种方法,即,将其夹入一本书中。当有两张一卡通时,就有了两种选择,即,将两张一卡通夹在一本书里,或者分开夹在不同的书里。当有三张一卡通时,他们就有了5种选择,即:
{{A},{B},{C}} , {{A,B},{C}}, {{B,C},{A}}, {{A,C},{B}} ,{{A,B,C}} 于是,
这个邪恶计划的组织者wf希望了解,如果ACM训练对里有n位帅哥(即有N张一卡通),那么要把这些一卡通夹到书里有多少种不同的方法。
解题思路:
转自:https://blog.csdn.net/qpswwww/article/details/44974693
第一类斯特林数s2[i][j]=将j个互不相同的物品划分成j个非空集合的方案数。
s2[i][j]=s2[i?1][j?1]+js2[i?1][j]
递推式的解释:对于第i个物品有两种情况:
1、前i?1个物品已经划分成了j?1个非空集合,第i个物品单独构成第j个集合。
2、前i?1个物品已经划分成了j个非空集合,第i个物品可以选择放入j个集合当中之一,共j个方案,因此是js2[i?1][j]。
贝尔数bell[i]=i个互不相同的物品,划分成若干个非空集合的方案数。
bell[i]=∑k=0is2[i][k]
这个式子非常显然,即枚举划分成k个非空集合,对s2[i][k]求和即可得到贝尔数。
代码
#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int N=2e3+;
const int MOD=1e3; int ans[N];
int dp[N][N];//dp[i][j]表示将前i个数分j组的方案数 int main(){
dp[][]=;
for(int i=;i<N;i++){
for(int j=;j<=i;j++){
dp[i][j]=(dp[i-][j-]+dp[i-][j]*j)%MOD;
}
}
for(int i=;i<N;i++){
for(int j=;j<=i;j++){
ans[i]=(ans[i]+dp[i][j])%MOD;
}
}
int n;
cin>>n;
while(n--){
int x;
cin>>x;
cout<<ans[x]<<endl;
}
return ;
}