CF GYM 1020561 I. Misunderstood … Missing解题报告（逆序动态规划）

题目链接：http://codeforces.com/gym/102056/problem/I
题目描述：

给定$D$，$A$两个初始值为$0$的属性，接下来进行$n$轮， 每轮有三个选项，你可以选择其一。三个选项为：

1. 对敌人造成$A+ai$的伤害；
2. 使自己的$A$增加$ci$；
3. 使自己的$D$增加$bi$。

每回合开始时$A$会增加$D$，问最多能造成多少伤害。

题解：

因为有三个选项对应了三种决择， 而考虑到贪心的不正确性， 我们可以考虑$dp$， 首先对于最后一回合时我们肯定选择攻击不会选择使$A$或者$D$增加， 但是$dp$又有一个问题：如果我们正着$dp$当我们选择$D$或者$A$增加时会对后面回合造成影响， 不满足$dp$的无后效性， 因为我们想到了最后一次必定攻击，我们考虑一下能不能倒着$dp$。
首先考虑如果空间够用， 我们尝试把攻击的位置进行状态压缩， $dp[i][j]$表示到$i$轮攻击状态为$j$能打出的最大伤害。如果本轮攻击， 那么当前就由$dp[i+1][k]+a[i]$（k为不包含$i$轮攻击的状态）转移过来（这里不增加$A$是因为将$A$的增加放在了后两种情况考虑）；如果选择$A$增加，那么之前处理过的伤害会增加，则可以由$dp[i+1][j]+k\ast c[i]$ （$k$为$j$与$i$的攻击次数之差）转移而来；如果我们这轮让$D$增加， 那么我们会在对应攻击的地方伤害会有增加我们设后续攻击的位置为$p1，p\mathrm{2......}pk$那么考虑伤害的增量：$|i-p1|\ast bi,|i-p2|\ast bi......|i-pk|\ast bi$我们把他们加起来后发现增量为:$|i\ast k\ast bi-(p1+p2+......+pk)\ast bi|$其中$k$为攻击次数， 那么我们发现选择$D$增加，伤害的增量与$k$(攻击次数），$bi$，$(p1+p2+......+pk)$有关，所以我们可以这样定义状态$dp[i][j][k]$表示到$i$轮，对应攻击次数$j$，攻击位置的和为$k$能打出的最大伤害，此时转移就变成了：$dp[i][j][k]=dp[i+1][j-i][k-i]+a[i]$（攻击）$dp[i][j][k]=dp[i+1][j][k]+c[i]\ast j$ (选择增加A）$dp[i][j][k]=dp[i+1][j][k]+|i\ast j\ast bi-k\ast bi|$（选择增加D）。
最后考虑空间问题由于只与$i$与$i+1$有关，所以滚动数组滚一滚。
代码：

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>

const int maxn = 110;

typedef long long ll;

using namespace std;

int n, t;
ll ans;
ll a[maxn], b[maxn], c[maxn], dp[2][maxn][(1+maxn)*maxn/2];

ll read(){
	ll x = 0, f = 1;
	char ch;
	do{
		ch = getchar();
		if (ch == '-') f = -1;
	}while(ch < '0' || ch > '9');
	do{
		x = x * 10 + ch - 48;
		ch = getchar();
	}while(ch >= '0' && ch <= '9');
	return x * f;
}
int main(){
	t = read();
	while(t --){
		n = read();
		for (int i = 1;i <= n;i ++) a[i] = read(), b[i] = read(), c[i]= read();
		int pos = 1, laspos = 0;
		memset(dp, -1, sizeof(dp));
		dp[0][1][n] = a[n];
		ans = a[n];
		for (int i = n-1;i >= 1;i --){
			memset(dp[pos], -1, sizeof(dp[pos]));
			for (int j = 1;j <= n-i+1;j ++){
				for (int k = n;k <= (n + i) * (n - i + 1) / 2;k ++){
					if (dp[laspos][j-1][k-i] != -1) dp[pos][j][k] = max(dp[pos][j][k], dp[laspos][j-1][k-i] + a[i]);
					if (dp[laspos][j][k] != -1) 
						dp[pos][j][k] = max(dp[pos][j][k], dp[laspos][j][k] + (k - j * i) * b[i]),
						dp[pos][j][k] = max(dp[pos][j][k], dp[laspos][j][k] + j * c[i]);
					ans = max(dp[pos][j][k], ans);
				}
			}
			laspos = pos;
			pos ^= 1;
		}
		cout << ans << endl;
	}
	return 0;
}