【区间dp】#558. 快快变大

题目链接

快快变大 - 题目 - Daimayuan Online Judge

题意

可以对数列进行的操作是合并 a[x]与a[x+1] 使之变成 a[x] * a[x+1] %mod 同时获得积分 (a[x]

-a[x+1])^2 直到该数列只剩下最后一个数组 求解 通过有限次的该操作 可以获得的最大积分是多少

思路

一道区间dp的题目  关键是表示出他的动态方程

我们用 dp[x][y] 表示合并从x到y 获得的最大积分数

那么动态方程为：

dp[i][j]=max(dp[i][j],dp[i][k]+dp[k+1][j]+cost)

cost指代把dp[i][k] 和 dp[k+1][j] 的两个区间合并起来的花费

则cost= (x-y)^2

其中 x=a[i]*..*a[k]         y=a[k+1]*...*a[j]

所以 我们需要用一个二维数组sum[x][y]来存储从 a[x] 乘到 a[y] 的结果

至此，整道题的分析就结束了

代码

#include<map>
#include<cmath>
#include<stack>
#include<queue>
#include<cstdio>
#include<vector>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#define Endl "\n"
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
const int maxn = 500;
const int INF = 0x3f3f3f;
const int mod =1000003 ;
const double pi = acos(-1); 
inline ll read(){
	ll x = 0, f = 1; char ch; ch = getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
ll n;
ll a[maxn];
ll dp[maxn][maxn];
ll sum[305][305];
ll x;
void solve() {
	n=read();
	memset(dp,0,sizeof(dp));
	for(int i=1;i<=n;i++)
	{
		a[i]=read();
	}
	for(int i=1;i<=n;i++)
	{
		sum[i][i]=a[i];
		for(int j=i-1;j>0;j--)
			sum[j][i]=(sum[j+1][i] * a[j] )%mod;	
	}
	for(int len=2;len<=n;len++)
		for(int i=1;i<=n;i++)
		{
			int j=i+len-1;
			if(j>n) break;
			for(int k=i;k<j;k++)
				dp[i][j]=max(dp[i][j],dp[i][k]+dp[k+1][j]+(sum[i][k]-sum[k+1][j])*(sum[i][k]-sum[k+1][j]));
		}
	cout<<dp[1][n]<<Endl;
	return ; 
} 

int main() {
//	for(int T = read(); T; T--)
		solve();
	return 0;
}