hdu 5900 QSC and Master（acm/icpc沈阳网络赛，区间DP）

传送门： http://acm.hdu.edu.cn/showproblem.php?pid=5900

题目大意：
给你一个数组，数组的每个元素是一个二元组，用<Key,Value>表示，你每次可以消掉相邻的两个数对，要求他们的key是互质的，你的得分是它们的value，消掉之后后面的会并过来，问你最多可以得到多少分?

题目分析：
首先看到“最多”字样，八成又是dp或者贪心。。。

这道题如果只考虑前i个数这样DP的话，因为每多一个数解法都有可能有很大变化，所以不能简单的用i-1去推i。
所以考虑区间[i,j],因为区间[i,j]要么可以全部约掉，要么从中间选一个最优的点k，使得约掉区间[i,k]和区间[k+1,j]区间得分是最大的。
那么令dp[i][j]代表[i,j]区间内的最优解，那么这个最优解只可能有两种情况，一是把他们全约掉。什么情况可以全约掉呢？很显然，[i+1,j-1]区间全能约掉，且两头的互素，则[i,j]必须全能约掉。

等等，我们的dp[i][j]只记录了最优解呀，没记录能不能全约掉？
别忘了，如果全约掉了，那么得分就是value[i]+…+value[j]嘛~ 要想O(1)的得到这个值，只需做一个前缀和就好啦~~

如果不能全约掉，那就枚举k，使得区间[i,k]和区间[k+1,j]区间得分和最大。
有两种搜法，一种是按长度求，先求所有长为1的区间，再求长为2的……
还可以是加记忆化的dfs，每当算[i,j]时，把[i,j]依赖的区间都先算出来，且算过的就不需要再算，因为[i,j]以外的区间不能影响[i,j]区间内的解。。。。
（亲，如果影响了就不符合动态规划的定义了哦~）

最后一个小tips，答案仍有可能爆int！！！！
时间复杂度应该是$O(n^3)$，N=300所以无压力。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
int t,n;
int k[305],v[305];
ll dp[305][305];
ll sum[305];
bool vis[305][305];
int gcd(int x,int y) {
    return y==0?x:gcd(y,x%y);
}
void dfs(int i,int j) {
    if(vis[i][j])
        return;
    vis[i][j]=true;
    dfs(i+1,j-1);
    if(dp[i+1][j-1] == sum[j-1]-sum[i] && gcd(k[i],k[j])!=1)
        dp[i][j]=dp[i+1][j-1]+v[i]+v[j];
    else {
        for(int k=i;k<j;k++) {
            dfs(i,k);
            dfs(k+1,j);
            dp[i][j]=max(dp[i][j],dp[i][k]+dp[k+1][j]);
        }
    }
}
int main() {
    scanf("%d",&t);
    while (t--) {
        scanf("%d",&n);
        sum[0]=0;
        for(int i=1;i<=n;i++)
            scanf("%d",&k[i]);
        for(int i=1;i<=n;i++) {
            scanf("%d",&v[i]);
            sum[i]=sum[i-1]+v[i];
        }

        memset(dp,0,sizeof(dp));
        memset(vis,0,sizeof(vis));
        for(int i=1;i<=n;i++)
            vis[i][i]=true;

        for(int i=1;i<n;i++) {
            vis[i][i+1]=true;
            if(gcd(k[i],k[i+1])==1)
                dp[i][i+1]=0;
            else
                dp[i][i+1]=v[i]+v[i+1];
        }
        dfs(1,n);
        //debug
        // for(int i=1;i<=n;i++)
        //     for(int j=1;j<=n;j++)
        //         printf("dp[%d][%d]=%I64d,vis=%d\n", i,j,dp[i][j],vis[i][j] );
        printf("%I64d\n", dp[1][n]);
    }
}

接下来随便扯两句。
其实啦~区间DP我们并不陌生，只是不知道这个名字罢了~~

记得大一的时候看的《计算机算法设计与分析》，在DP那章的第一题，求最小的矩阵连乘，那道题就是很典型的区间DP，而且也是在[i,j]内枚举分点k，从[i,k]和[k+1,j]区间里推出[i,j].

顺便，我能吐槽下出题者的英语水平吗？看见好几个病句。。。。。。。。。