具体说来就是分成1组,分成2组,分成3组……分成n组的数量全部加起来
设m[i]表示分成i组的有多少人 dp[i]表示当n=i的时候的分法
则当多出一个i时,这个i可以和别人一组,也可以自己一组
自己一组的分法有dp[i-1]个,和别人一起的分法根据分成k组有m[k]个分发,则就有km[k]个分发法(对于m[k]中的每个分法,第i个可以到k组里的任意一个)
所以dp[i]=dp[i-1]+sigma(km[k])(k:1->i-1)
然后就维护m数组,只要照上面的意思维护就可以了
m[k]=i*temp[k]+temp[k-1] (k:1->i)
temp[k]=m[k] (k:1->i);
#include <iostream>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<math.h>
#include<vector>
#define ENTER printf("\n");
//快速打印一个数组
#define ll long long
#define NUM 100
#define MOD 1000
using namespace std;
//构造函数 析构函数怎么打印出地址
int n;
int m[2000 + 1];
int temp[2000 + 1];
int dp[2000 + 1];
int t;
int main()
{
dp[0] = 0; dp[1] = 1;
m[0] = 0; m[1] = 1;
temp[0] = 0; temp[1] = 1;
for (int i = 2; i <= 2000; i++)
{
for (int j = 1; j <= i - 1; j++)
{
dp[i] += (j * m[j]) % MOD;
}
dp[i] = (dp[i] + dp[i - 1]) % MOD;
for (int j = 1; j <= i; j++)
{
m[j] = ((j * temp[j]) % MOD + temp[j - 1] % MOD) % MOD;
}
for (int j = 1; j <= i; j++)
{
temp[j] = m[j];
}
}
scanf("%d", &t);
while (t--)
{
scanf("%d", &n);
printf("%d\n", dp[n]);
}
return 0;
}