四边形不等式（石子合并）

 动态规区间dp做这道题的话应该是n^3，下面的代码优化到了n^2，用四边形不等式优化。

设mid[i][j]是dp[i][j]的最优解的断点，即它左区间的右端点，那么mid[i][j-1]<=mid[i][j]<=mid[i+1][j]，所以在求解dp[i][j]时，枚举k可以只枚举这两个值之间枚举就好，

程序要先枚举区间长度，在枚举左端点，枚举每个区间长度时，他们的k总是只从1到n，只走一遍，所以这就相当于优化了一层，变成了O(n²)的。

比如len长度为3时，dp[1][3]只会枚举mid[1][2]-mid[2][3]，如上。然后dp[2][4]时，枚举mid[2][3]-mid[3][4]。以此类推。

#include<cstdio>
#include<algorithm>
using namespace std;
const int MAXN = 4010;
const int INF = 1e8;
int dp[MAXN][MAXN],a[MAXN],sum[MAXN],mid[MAXN][MAXN],n;

int main()
{
    scanf("%d",&n);
    for (int i=1; i<=n; ++i)
        scanf("%d",&a[i]), sum[i] = sum[i-1]+a[i];
    
    for (int i=1; i<=n; ++i)
        dp[i][i] = 0, mid[i][i] = i;
    
    for (int len=1; len<n; ++len)//求右端点(i+len-1)会减1，先在这减了 
    {
        for (int i=1; i<=n-len; ++i)//左端点i最大值是右端点等于n时，(i+len)=n,i=n-len 
        {
            int j = i+len;//右端点 
            dp[i][j] = INF;
            for (int k=mid[i][j-1]; k<=mid[i+1][j]; ++k) 
            {
                if (dp[i][j]>dp[i][k]+dp[k+1][j])
                {
                    dp[i][j] = dp[i][k]+dp[k+1][j];
                    mid[i][j] = k;
                }
            }
            dp[i][j] += sum[j]-sum[i-1];
        }
    }
    printf("%d",dp[1][n]);
    return 0;
}