区间dp

1022 石子归并 V2

N堆石子摆成一个环。现要将石子有次序地合并成一堆。规定每次只能选相邻的2堆石子合并成新的一堆，并将新的一堆石子数记为该次合并的代价。计算将N堆石子合并成一堆的最小代价。

 

例如： 1 2 3 4，有不少合并方法

1 2 3 4 => 3 3 4(3) => 6 4(9) => 10(19)

1 2 3 4 => 1 5 4(5) => 1 9(14) => 10(24)

1 2 3 4 => 1 2 7(7) => 3 7(10) => 10(20)

 

括号里面为总代价可以看出，第一种方法的代价最低，现在给出n堆石子的数量，计算最小合并代价。

输入

第1行：N（2 <= N <= 1000)
第2 - N + 1：N堆石子的数量（1 <= A[i] <= 10000)

输出

输出最小合并代价

输入样例

4
1
2
3
4

输出样例

19

 

区间dp，四边形不等式优化 O(n^2)

#include<bits/stdc++.h>
using namespace std;
const int maxn = 2009;
long long dp[maxn][maxn];
long long s[maxn][maxn];
long long a[maxn];
long long sum[maxn];
int main()
{
    long long i,j,k,n,m,ans;
    scanf("%lld",&n);
    memset(dp,0x3f,sizeof(dp));
    for(i = 1; i <= 2*n; i++)//开两倍数组模拟环
    {
        if(i<=n)
            scanf("%lld",&a[i]);
        else
            a[i] = a[i-n];
        dp[i][i] = 0;
        s[i][i] = i;
        sum[i] = sum[i-1]+a[i];
    }


    for(i = 2; i<=n; i++)//枚举区间长度
    {
        for(j = 1; j+i<=2*n; j++)//枚举起点
        {
            int ends = j+i-1;
            for(k = s[j][ends-1]; k<= s[j+1][ends]; k++)//枚举中间点，四边形不等式优化
            {
                if(dp[j][k]+dp[k+1][ends]+sum[ends]-sum[j-1]<dp[j][ends])
                {
                    dp[j][ends] = dp[j][k]+dp[k+1][ends]+sum[ends]-sum[j-1];
                    s[j][ends] = k;
                }
            }
        }
    }
    ans = 0x3f3f3f3f;
    for(i = 1; i<=n; i++)
    {
        ans = min(ans,dp[i][i+n-1]);
    }
    printf("%lld\n",ans);
    return 0;
}