2018ACM-ICPC EC-Final 现场赛I题 Misunderstanding...Missing 倒着DP

（有任何问题欢迎留言或私聊 && 欢迎交流讨论哦

Catalog

Problem:

 很多acm群里有题面PDF了，我就不赘述了。简单说一下，你有n次操作，每次操作有3种选择，1.造成A+ai点伤害；2.永久给D增加bi；3.永久给A增加ci。(每次操作前执行：A+=D)。问最后最多造成多少伤害。

Solution：

 这么明显的dp却没有想到，真实太菜了。当时想到要从后往前推，但就是卡着不知道怎么转移最优，如何平衡3种操作。赛后看到qls的状态定义，终于理解qls的意思了。菜是原罪啊。

状态表示：
$dp[i][j][k]$表示从后开始执行到第$i$步，执行了$j$次操作1，选择操作1的下标和为$k$所造成的最高伤害。

初始化：

memset(dp, -1, sizeof(dp));//必须初始化为-1，因为你要确保状态的正确性，也就是j次操作1的下标和真的是k才行。
dp[n][1][n] = cw[n].a;//因为最后一次肯定是要打伤害的

状态转移方程：
$$dp[i][j][k]=max(dp[i+1][j][k]+max(c\ast j,b\ast (k-i\ast j)),dp[i+1][j-1][k-i]+a);$$

• 第一部分是第$i$次选择操作2或3，第二部分是第$i$次选择操作1。
• 如果这一次选择操作3，显然后面的$j$次伤害需要增加$c\ast j$。
• 同理：选择操作2，后面每次操作需要增加$b\ast (x-i)$点伤害，$x$是选择操作1的编号，求个和之后就是$b\ast (k-i\ast j)$了。
• 到这里，这题就解决啦。
  
  

AC_Code：

第一个代码爆内存了，其他两个都可以ac

//这份代码常数应该挺大，注意优化一下
#include<bits/stdc++.h>

using namespace std;
typedef long long LL;

const int MXN = 1e5 + 5;

int n, m;
LL dp[101][101][5052];
struct lp{
    int a, b, c;
}cw[105];

int main(int argc, char const *argv[]) {
    int tim; scanf("%d", &tim);
    while(tim --) {
        scanf("%d", &n);
        for(int i = 1; i <= n; ++i) {
            scanf("%d%d%d", &cw[i].a, &cw[i].b, &cw[i].c);
        }
        memset(dp, -1, sizeof(dp));
        dp[n][1][n] = cw[n].a;
        int tmp = n*(n+1)/2;
        for(int i = n-1; i >= 1; --i) {
            for(LL j = n-i+1,a=cw[i].a,b=cw[i].b,c=cw[i].c; j >= 1; --j) {
                for(int k = tmp; k >= n; --k) {
                    if(dp[i+1][j][k] != -1) dp[i][j][k] = dp[i+1][j][k]+max(c*j, b*(k-i*j));
                    if(k >= i+n &&j >= 2&& dp[i+1][j-1][k-i] != -1)
                        dp[i][j][k] = max(dp[i][j][k], dp[i+1][j-1][k-i]+a);
                }
            }
        }
        LL ans = 0;
        for(int i = 1; i <= n; ++i) {
            for(int j = 1; j <= tmp; ++j) {
                ans = max(ans, dp[1][i][j]);
            }
        }
        printf("%lld\n", ans);
    }
    return 0;
}

滚动数组：

#include<bits/stdc++.h>
#define A cw[i].a
#define B cw[i].b
#define C cw[i].c
using namespace std;
typedef long long LL;
const int INF = 0x3f3f3f3f;
const int MXN = 1e5 + 7;

int n, m;
struct lp{
    LL a, b, c;
}cw[105];
LL dp[5051][101];

int main() {
    int tim; scanf("%d", &tim);
    while(tim --) {
        scanf("%d", &n);
        m = n*(n+1)/2;
        for(int i = 0; i <= m; ++i) {
            for(int j = 0; j <= n; ++j) {
                dp[i][j] = -1;
            }
        }
        for(int i = 1; i <= n; ++i) scanf("%lld%lld%lld", &A, &B, &C);
        dp[n][1] = cw[n].a;
        //dp[0][0] = 0;
        for(int i = n-1; i >= 1; --i) {
            for(int k = m; k >= n; --k) {
                for(int j = n-i+1; j >= 1; --j) {
                    if(dp[k][j] != -1) dp[k][j] += max(j*C,(k-i*j)*B);
                    if(k-i>=n&&j>=2&&dp[k-i][j-1] != -1) {
                        dp[k][j] = max(dp[k][j], dp[k-i][j-1] + A);
                    }
                }
            }
        }
        LL ans = 0;
        for(int i = 1; i <= m; ++i) {
            for(int j = 1; j <= n; ++j) {
                ans = max(ans, dp[i][j]);
            }
        }
        printf("%lld\n", ans);
    }    
    return 0;
}

如果你不会滚动数组，也可以这样优化空间：

#include<bits/stdc++.h>
#define A cw[i].a
#define B cw[i].b
#define C cw[i].c
using namespace std;
typedef long long LL;
const int INF = 0x3f3f3f3f;
const int MXN = 1e5 + 7;

int n, m;
struct lp{
    LL a, b, c;
}cw[105];
LL dp[5051][101][2];

int main() {
    int tim; scanf("%d", &tim);
    while(tim --) {
        scanf("%d", &n);
        m = n*(n+1)/2;
        for(int i = 0; i <= m; ++i) {
            for(int j = 0; j <= n; ++j) {
                dp[i][j][0] = dp[i][j][1] = -1;
            }
        }
        for(int i = 1; i <= n; ++i) scanf("%lld%lld%lld", &A, &B, &C);
        dp[n][1][0] = cw[n].a;
        //dp[0][0][0] = 0;
        int p = 1, q = 0;
        for(int i = n-1; i >= 1; --i) {
            for(int k = m; k >= n; --k) {
                for(int j = n-i+1; j >= 1; --j) {
                    if(k-i>=n&&j>=2&&dp[k-i][j-1][q] != -1) dp[k][j][p] = dp[k-i][j-1][q] + A;
                    if(dp[k][j][q] != -1) {
                        dp[k][j][p] = max(dp[k][j][p], dp[k][j][q]+max(j*C,(k-i*j)*B));
                    }
                }
            }
            p = !p; q = !q;
        }
        LL ans = 0;
        for(int i = 1; i <= m; ++i) {
            for(int j = 1; j <= n; ++j) {
                ans = max(ans, dp[i][j][q]);
            }
        }
        printf("%lld\n", ans);
    }    
    return 0;
}

Problem Description：