2017年浙工大程序设计迎新赛决赛—网络同步赛 G-取数游戏2（区间dp）

G-取数游戏2

时间限制：C/C++ 1秒，其他语言2秒
空间限制：C/C++ 131072K，其他语言262144K
64bit IO Format: %lld

题目描述

给定两个长度为n的整数列A和B，每次你可以从A数列的左端或右端取走一个数。假设第i次取走的数为ax，则第i次取走的数的价值vi=bi⋅ax，现在希望你求出∑vi的最大值。

输入描述:

第一行一个数T，表示有T组数据。
对于每组数据，第一行一个整数n，
接下来两行分别给出A数列与B数列。

输出描述:

每一组数据输出一行，最大的∑vi。

示例1

输入

2
2
1 1000
2 1
5
1 3 5 2 4
1 2 3 4 5

输出

2001
52

说明

对于第二个样例，
第一次从左边取走a1，v1=a1⋅b1=1,
第二次从左边取走a2，v2=a2⋅b2=6,
第三次从右边取走a5，v3=a5⋅b3=12,
第四次从右边取走a4，v4=a4⋅b4=8,
第五次取走剩下的a3，v5=a3⋅b5=25。
总价值∑vi=1+6+12+8+25=52

备注:

T≤10
1≤n≤103
1≤ai,bi≤103

问题分析

蒟蒻的第一道区间dp，过了好久，还是写一下题解吧。
区间dp就如其名，找状态自然要从区间里找。题里要求每次从a数列左面或者右面取一个数，但是会乘上bi，这就导致不同的取法会造成不同的结果。所以如何确定最优子结构呢？我们可以把问题分解为求当前子区间可获得的最大价值。假设当前子区间可获得的最大值确定之后，我们再去确定其它区间可获得的最大值。最后dp[1][n]就是答案。每次选择都有两个选择，取左还是取右，都去比较一下就好，用dp[j][k]保存当前区间可获得的最大价值。边界条件为区间起点等于终点时，dp[j][k] = a[i]*b[n]（j==k，1<=i<=n）。
dp[j][k]=max(dp[j+1][k]+a[j]*b[n-i],dp[j][k-1]+a[k]*b[n-i])|1<=i<=n,1<=j<=n-i

好，接下来上AC code（^_^）

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int N = 1e3+2;
int a[N],b[N],dp[N][N];

int main()
{
    // freopen("in.txt","r",stdin);
    int t,n;
    cin>>t;
    while(t--)
    {
        cin>>n;
        memset(dp,0,sizeof dp);
        for(int i = 1; i <= n; i++)
            cin>>a[i];
        for(int j = 1; j <= n; j++)
            cin>>b[j];
        for(int i=1;i<=n;i++) //边界
            dp[i][i]=a[i]*b[n];
        for(int i = 1; i <= n; i++) //枚举区间长度
        {
            for(int j = 1; j+i <= n; j++)  //枚举区间起点
            {
                int k = i+j;  //区间终点
                //当前区间可获得的最大价值
                dp[j][k]=max(dp[j+1][k]+a[j]*b[n-i],dp[j][k-1]+a[k]*b[n-i]);//取左面还是右面
            }
        }
        cout<<dp[1][n]<<endl;
    }
    return 0;
}

其实，这道题还可以拿记忆化搜索来做。（这样或许更好理解一些）

int solve(int j,int k,int x)
{
    if(j==k)//边界
    return a[j]*b[x];
    if(dp[j][k]!=0)
    return dp[j][k];
    //记忆化(当前区间取左还是取右)
    return dp[j][k] = max(solve(j+1,k,x+1)+a[j]*b[x],solve(j,k-1,x+1)+a[k]*b[x]);
}

最后输出solve(1,n,1)即可。