石子合并（四边形优化第一题）

满足两个性质

 

1、区间包含的单调性：如果对于 i≤i'<j≤j'，有 w(i',j)≤w(i,j')，那么说明w具有区间包含的单调性。（可以形象理解为如果小区间包含于大区间中，那么小区间的w值不超过大区间的w值）

2、四边形不等式：如果对于 i≤i'<j≤j'，有 w(i,j)+w(i',j')≤w(i',j)+w(i,j')，我们称函数w满足四边形不等式。（可以形象理解为两个交错区间的w的和不超过小区间与大区间的w的和）

 

下面给出两个定理：

　　　　1、如果上述的 w 函数同时满足区间包含单调性和四边形不等式性质，那么函数 m 也满足四边形不等式性质

　　　　　  我们再定义 s(i,j) 表示 m(i,j) 取得最优值时对应的下标（即 i≤k≤j 时，k 处的 w 值最大，则 s(i,j)=k）。此时有如下定理

　　　　2、假如 m(i,j) 满足四边形不等式，那么 s(i,j) 单调，即 s(i,j)≤s(i,j+1)≤s(i+1,j+1)。

然后存一个s[i][j]表示这个区间最优值的下标，然后以后每次就k就不用枚举了

 

for(k=s[i][j-1];k<=s[i+1][j];k++){
                if(dp2[i][j]<dp2[i][k]+dp2[k+1][j]){
                    dp2[i][j]=dp2[i][k]+dp2[k+1][j];
                    s[i][j]=k;
                }
            }

这样即可

 

全部代码

 

#include<cstdio>
#include<cmath>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<vector>
#include<set>
#include<queue>
#include<iostream>
using namespace std;

int dp1[105][105];
int dp2[105][105];
int sum[105];
int qq[105];
int s[105][105];
int main(){
    freopen("in.txt","r",stdin);
    int i,j,k,f1,f2,f3,f4,t1,t2,t3,t4,n,m;
    cin >> n;
    for(i=1;i<=n;i++){
        cin >> sum[i];
        sum[i]+=sum[i-1];
    }
    int len1;
    for(i=0;i<=n+1;i++)s[i][i]=i;
    for(len1=1;len1<n;len1++){
        for(i=1;i<n;i++){
            j=i+len1;
            if(j>n)break;
            dp1[i][j]=1e9+7;
            for(k=s[i][j-1];k<=s[i+1][j];k++){
                if(dp1[i][j]>dp1[i][k]+dp1[k+1][j]){
                    dp1[i][j]=dp1[i][k]+dp1[k+1][j];
                    s[i][j]=k;
                }
            }
            dp1[i][j]+=sum[j]-sum[i-1];
        }
    }
    for(i=0;i<=n+1;i++)s[i][i]=i;
    for(len1=1;len1<n;len1++){
        for(i=1;i<n;i++){
            j=i+len1;
            if(j>n)break;
            //for(k=i;k<j;k++)
            //dp[i][k],dp[k+1][j]
            for(k=s[i][j-1];k<=s[i+1][j];k++){
                if(dp2[i][j]<dp2[i][k]+dp2[k+1][j]){
                    dp2[i][j]=dp2[i][k]+dp2[k+1][j];
                    s[i][j]=k;
                }
            }
            dp2[i][j]+=sum[j]-sum[i-1];
        }
    }
    cout << dp1[1][n] << endl;
    cout << dp2[1][n] << endl;
    return 0;
}