题目链接:https://vjudge.net/problem/HDU-2512
一卡通大冒险
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 2572 Accepted Submission(s): 1741
Problem Description
因为长期钻研算法, 无暇顾及个人问题,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张一卡通),那么要把这些一卡通夹到书里有多少种不同的方法。
{{A},{B},{C}} , {{A,B},{C}}, {{B,C},{A}}, {{A,C},{B}} ,{{A,B,C}} 于是,
这个邪恶计划的组织者wf希望了解,如果ACM训练对里有n位帅哥(即有N张一卡通),那么要把这些一卡通夹到书里有多少种不同的方法。
Input
包含多组数据,第一行为n,表示接下来有n组数据。以下每行一个数x,表示共有x张一卡通。(1≤x≤2000).
Output
对每组数据,输出一行:不同的方法数,因为这个数可能非常大,我们只需要它除以1000的余数。
Sample Input
4
1
2
3
100
1
2
3
100
Sample Output
1
2
5
751
2
5
751
Author
BUAA Campus 2007
Source
题解:
单纯的第二类斯特林数。由于没有要求组数,因此对 S[n][k] 求和, 其中1<=k<=n 。
代码一:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <cmath>
#include <queue>
#include <stack>
#include <map>
#include <string>
#include <set>
using namespace std;
typedef long long LL;
const int INF = 2e9;
const LL LNF = 9e18;
const int MOD = 1e3;
const int MAXN = 2e3+; LL S[MAXN][MAXN], f[MAXN]; void init()
{
// memset(S, 0, sizeof(S));
for(int i = ; i<MAXN; i++)
{
S[i][] = ; S[i][i] = ;
for(int j = ; j<i; j++)
S[i][j] = ((j*S[i-][j])%MOD + S[i-][j-])%MOD;
} memset(f, , sizeof(f));
for(int i = ; i<MAXN; i++)
for(int j = ; j<=i; j++)
f[i] = (f[i] + S[i][j])%MOD;
} int main()
{
int T, n;
scanf("%d", &T);
init();
while(T--)
{
scanf("%d", &n);
printf("%d\n", f[n]);
}
}
代码二:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <cmath>
#include <queue>
#include <stack>
#include <map>
#include <string>
#include <set>
using namespace std;
typedef long long LL;
const int INF = 2e9;
const LL LNF = 9e18;
const int MOD = 1e3;
const int MAXN = 2e3+; LL S[][MAXN], f[MAXN]; void init()
{
int cur = ;
memset(f, , sizeof(f));
for(int i = ; i<MAXN; i++)
{
S[cur][] = ; S[cur][i] = ;
for(int j = ; j<i; j++)
S[cur][j] = ((j*S[!cur][j])%MOD + S[!cur][j-])%MOD; for(int j = ; j<=i; j++)
f[i] = (f[i] + S[cur][j])%MOD; cur ^= ;
}
} int main()
{
int T, n;
scanf("%d", &T);
init();
while(T--)
{
scanf("%d", &n);
printf("%d\n", f[n]);
}
}