代码:石子合并

石子合并
链接：https://ac.nowcoder.com/acm/problem/51170
来源：牛客网

题目描述:
设有N堆沙子排成一排，其编号为1，2，3，…… ，N（N≤300）。每堆沙子有一定的数量，可以用一个整数来描述，现在要将这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)。

输出描述:
合并的最小代价

区间DP

#include <bits/stdc++.h>
using namespace std;
const int N = 10000001;
int value[1001];
int dp[1001][1001];
int main()
{
    int n;
    cin >> n;
    for (int i = 1; i <= n; i++)
        for (int j = 1; j <= n; j++)
            if (i != j)
                dp[i][j] = N;
    int temp;
    value[0] = 0;
    for (int i = 1; i <= n; ++i)
    {
        cin >> temp;
        //前缀和处理
        value[i] = value[i - 1] + temp;
        dp[i][i] = 0;
    }
    //区间dp
    for (int l = 1; l <= n; ++l)
    //枚举长度
        for (int i = 1; i <= n; ++i)
        //枚举起点
        {
            int j = i + l ;
            if (j > n)
                break;
            for (int k = i; k < j; ++k)
                dp[i][j] = min(dp[i][j], dp[i][k] + dp[k + 1][j] + value[j] - value[i - 1]);
                //状态转移方程
        }
    cout << dp[1][n] << endl;
}

四边形不等式优化

#include <bits/stdc++.h>
using namespace std;
const int N = 10000001;
int value[1001];
int dp[1001][1001];
int minn[1001][1001];
int main()
{
    int n;
    cin >> n;
    for (int i = 1; i <= n; i++)
        for (int j = 1; j <= n; j++)
            if (i != j)
                dp[i][j] = N;
    int temp;
    value[0] = 0;
    for (int i = 1; i <= n; ++i)
    {
        cin >> temp;
        //前缀和处理
        value[i] = value[i - 1] + temp;
        dp[i][i] = 0;
        minn[i][i] = i;
    }
    //区间dp
    for (int l = 1; l <= n; ++l)
        //枚举长度
        for (int i = 1; i <= n; ++i)
        //枚举起点
        {
            int j = i + l;
            if (j > n)
                break;
            for (int k = minn[i][j - 1]; k <= minn[i + 1][j]; ++k)
            //四边形不等式优化
                if (dp[i][k] + dp[k + 1][j] + value[j] - value[i - 1] < dp[i][j])
                {
                    dp[i][j] = dp[i][k] + dp[k + 1][j] + value[j] - value[i - 1];
                    minn[i][j] = k;
                }
        }
    cout << dp[1][n] << endl;
}