题目链接
快快变大 - 题目 - 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;
}