题目链接
http://acm.hdu.edu.cn/showproblem.php?pid=5900
题意:给出序列Ai.key和Ai.value,若当前相邻的两个数Ai.key和Ai+1.key的最大公约数大于1,则可以把这两个数消去,同时消去Ai.value和Ai+1.value,每次消去得到的分数为Ai和Ai+1的value值,问最大可能得分。
注意:当Ai和Ai+1被消去后,Ai-1和Ai+2成为了新的相邻的数。若符合条件则可以继续消去。
思路:很明显是区间DP,但是我比赛中如何也A不了。原因有两个:1.没有注意到位运算的低优先级。2.没有把所有情况考虑清楚。
设f[i][j]为从第i个数到第j个数所能得到的最大得分,j-i+1从2开始依次枚举,最后枚举到f[1][n](考虑DP的无后效性)
分类讨论即可。先考虑不合并的情况,则f[i][j]被f[i][k]+f[k + 1][j]依次更新。
再考虑合并的情况f[i][j]可以被f[i][k] + a[k + 1] + a[k + 2] + f[k +3][j]更新。
特殊情况讨论:当a[i + 1]到a[j - 1]中所有的数都可以被消去,则a[i]和a[j]作为相邻的数也可以被消去。这种情况的讨论非常重要!!!
代码送上
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <functional>
using namespace std;
#define REP(i,n) for(int i(0); i < (n); ++i)
#define rep(i,a,b) for(int i(a); i <= (b); ++i)
#define dec(i,a,b) for(int i(a); i >= (b); --i)
#define for_edge(i,x) for(int i = H[x]; i; i = X[i])
const int N = 100000 + 10;
const int M = 10000 + 10;
const int Q = 1000 + 10;
const int A = 30 + 1;
typedef long long LL;
LL f[Q][Q];
bool c[Q][Q];
bool ret[Q][Q];
int T;
int n;
LL a[Q], b[Q];
LL gcd(LL a, LL b){ return b == 0 ? a : gcd(b, a % b); }
int main(){
scanf("%d", &T);
while (T--){
scanf("%d", &n);
memset(a, 0, sizeof a);
memset(b, 0, sizeof b);
memset(c, false, sizeof c);
memset(ret, false, sizeof ret);
memset(f, 0, sizeof f);
rep(i, 1, n) scanf("%lld", a + i);
rep(i, 1, n) scanf("%lld", b + i);
rep(i, 1, n - 1) rep(j, i + 1, n) if (gcd(a[i], a[j]) != 1) c[i][j] = true;
rep(i, 1, n - 1) if (gcd(a[i], a[i + 1]) > 1) ret[i][i + 1] = true;
for (int len = 4; len <= n; len += 2){
rep(i, 1, n - len + 1){
int j = i + len - 1;
bool flag = false;
if (c[i][i + 1] && ret[i + 2][j]) flag = true;
if (c[j - 1][j] && ret[i][j - 2]) flag = true;
if (c[i][j] && ret[i + 1][j - 1]) flag = true;
for (int k = i + 2; k <= j - 3; k += 2) if (ret[k][k + 1] && ret[i][k - 1] && ret[k + 2][j]){ flag = true; break;}
if (flag) ret[i][j] = true;
}
}
rep(i, 1, n - 1) if (c[i][i + 1]) f[i][i + 1] = b[i] + b[i + 1];
rep(i, 1, n - 2) f[i][i + 2] = max(f[i][i + 1], f[i + 1][i + 2]);
rep(len, 4, n){
rep(i, 1, n - len + 1){
int j = i + len - 1;
rep(k, i, j) f[i][j] = max(f[i][j], f[i][k] + f[k + 1][j]);
if (len % 2 == 0){ if (ret[i + 1][j - 1] && c[i][j]) f[i][j] = max(f[i][j], b[i] + b[j] + f[i + 1][j - 1]); }
f[i][j] = max(f[i][j], f[i + 1][j]);
f[i][j] = max(f[i][j], f[i][j - 1]);
if (c[i][i + 1]) f[i][j] = max(f[i][i + 1] + f[i + 2][j], f[i][j]);
if (c[j - 1][j]) f[i][j] = max(f[j - 1][j] + f[i][j - 2], f[i][j]);
rep(k, i + 1, j - 2) if (c[k][k + 1]) f[i][j] = max(f[i][j], f[i][k - 1] + f[k][k + 1] + f[k + 2][j]);
}
}
printf("%lld\n", f[1][n]);
}
return 0;
}