hdu 5000 dp+猜

题意：

给n（2000）个T [ i ]，1 <= T [ i ]，并且所有T[ i ]的和小于等于2000。

每个T代表着一个能力，可以取值为[ 0 , T[i] ]。

每个人有n个能力，如果这个能力的值全部都大于等于另一个人的能力，则另一个人就不能生存了。

现在给出所有能力的最大值T[i]，问最多能留下多少人共存。

解析：

想要共存，就要把这些人都凑开，一项能力大，一项能力小这样子。

然后开始我一直wa，因为猜想的地方错了，开始的猜想是，所有的组合每个人的能力值加起来为所有能力最大的那个这样子最大。

今天看了下题解，就改了一个地方就过了，所有的组合，每个人的能力值加起来为所有T[i]的和 sum / 2，这个时候才是最多的。

然后就是dp的地方了，dp[ i ] [ j ]表示，前 i 个数，凑成 j 所用的方法数。

状态转移方程：

dp[i][j] = (dp[i][j] + dp[ i - 1 ] [ j - k ]) % mod

其中，k为 0 到 T[i]。

代码：

#pragma comment(linker, "/STACK:1677721600")
#include <map>
#include <set>
#include <cmath>
#include <queue>
#include <stack>
#include <vector>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <climits>
#include <cassert>
#include <iostream>
#include <algorithm>
#define pb push_back
#define mp make_pair
#define LL long long
#define lson lo,mi,rt<<1
#define rson mi+1,hi,rt<<1|1
#define Min(a,b) ((a)<(b)?(a):(b))
#define Max(a,b) ((a)>(b)?(a):(b))
#define mem(a,b) memset(a,b,sizeof(a))
#define FIN freopen("in.txt", "r", stdin)
#define FOUT freopen("out.txt", "w", stdout)
#define rep(i,a,b) for(int i=(a); i<=(b); i++)
#define dec(i,a,b) for(int i=(a); i>=(b); i--)

using namespace std;
const int mod = 1e9 + 7;
const double eps = 1e-4;
const double ee = exp(1.0);
const int inf = 0x3f3f3f3f;
const int maxn = 2000 + 10;
const double pi = acos(-1.0);
const LL iinf = 0x3f3f3f3f3f3f3f3f;

LL readT()
{
    char c;
    int flg = 0;
    LL ret = 0;
    while(c = getchar(), (c < '0' || c > '9') && c != '-');
    if(c == '-') flg = 1;
    else ret = c ^ 48;
    while( c = getchar(), c >= '0' && c <= '9') ret = ret * 10 + (c ^ 48);
    return flg ? - ret : ret;
}

int a[maxn];
LL dp[maxn][maxn];

int main()
{
#ifdef LOCAL
    FIN;
#endif // LOCAL
    int ncase;
    scanf("%d", &ncase);
    while (ncase--)
    {
        int n;
        scanf("%d", &n);
        int sum = 0;
        rep(i, 1, n)
        {
            scanf("%d", &a[i]);
            sum += a[i];
        }
        sum /= 2;
        mem(dp, 0);
        dp[0][0] = 1;
        rep(i, 1, n)
        {
            rep(j, 0, sum)
            {
                dec(k, a[i], 0)
                {
                    if (j - k < 0)
                        continue;
                    dp[i][j] = (dp[i][j] + dp[i - 1][j - k]) % mod;
                }
            }
        }
        printf("%I64d\n", dp[n][sum] % mod);
    }
    return 0;
}