石子合并(四边形不等式优化)

石子合并

题目链接
设有N堆石子排成一排，其编号为1，2，3，…，N。
每堆石子有一定的质量，可以用一个整数来描述，现在要将这N堆石子合并成为一堆。
每次只能合并相邻的两堆，合并的代价为这两堆石子的质量之和，合并后与这两堆石子相邻的石子将和新堆相邻，合并时由于选择的顺序不同，合并的总代价也不相同。
例如有4堆石子分别为 1 3 5 2， 我们可以先合并1、2堆，代价为4，得到4 5 2， 又合并 1，2堆，代价为9，得到9 2 ，再合并得到11，总代价为4+9+11=24；
如果第二步是先合并2，3堆，则代价为7，得到4 7，最后一次合并代价为11，总代价为4+7+11=22。
问题是：找出一种合理的方法，使总的代价最小，输出最小代价。

输入格式

第一行一个数N表示石子的堆数N。
第二行N个数，表示每堆石子的质量(均不超过1000)。

输出格式

输出一个整数，表示最小代价。

数据范围

1≤N≤300

输入样例：

4
1 3 5 2

输出样例：

22

算法分析

朴素做法的区间dp 时间复杂度为O(n³),当n大于1000的时候就很容易超时
所以我们使用四边形不等式优化,能用到这个优化,要满足,求和需要取最优中间结点
而且具有单调性就是小区间如果包含于大区间,大区间一定比小区间大
还有就是如果如果对于i≤i+1<j≤j+1，有w(i,j)+w(i+1,j+1)≤w(i+1,j)+w(i,j+1)，我们称函数w满足四边形不等式。
w为区间代价（可以形象理解为两个交错区间的w的和不超过小区间与大区间的w的和）
这样我们的时间复杂度降到了O(n²)
具体的四边形不等式优化参考

代码实现

#include<iostream>
#include<algorithm>
using namespace std;
int n;
const int maxn=1e3+5;
int dp[maxn][maxn];
int p[maxn][maxn];      //p[i][j]表示i到j最优解的分界点位置
int s[maxn];
int main()
{
    cin>>n;
    for(int i=1;i<=n;i++)
    {
        cin>>s[i];
        s[i]=s[i]+s[i-1];
        p[i][i]=i;
    }
    for(int len=2;len<=n;len++)
    {
        for(int i=1;i+len-1<=n;i++)
        {
            int j=i+len-1;
            dp[i][j]=1e9;
            for(int k=p[i][j-1];k<=p[i+1][j];k++)
            {
                int dj=dp[i][k]+dp[k+1][j]+s[j]-s[i-1];
                if(dp[i][j]>=dj)
                {
                    dp[i][j]=dj;
                    p[i][j]=k;
                }
            }
        }
    }
    cout<<dp[1][n]<<endl;
    return 0;
}