石子合并 区间动态规划

设有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；问题是：找出一种合理的方法，使总的代价最小。输出最小代价）

1.f[L][r]表示从点L到点r的石子合并的最小花费，sum[i]表示前i个沙堆的数量和。
2.将区间[L,r]分成两段[L,k],[k+1,r],这两段分别需要花费f[L][k]和f[k+1][r],再将这两段合并需要sum[L]-sum[r-1]
3.f[L][r]=min(f[L][k]+f[k+1][r]); f[L][r]+=sum[L]-sum[r-1];

区间动态规划：动态规划分为阶段、状态和决策，三者应该从外到里循环。已知区间长为阶段，那么状态就是左端点和右端点，而决策就是不断地枚举合并点k的位置，找寻最优解

#include<bits/stdc++.h>
using namespace std;
const int inf=0x3f3f3f3f;
int a[1005],sum[1005]={0};
int f[1005][1005];

int main()
{
    int n;
    cin>>n;
    for(int i=1;i<=n;i++){
        cin>>a[i];
    }
    //初始化赋值为最大值
    memset(f,inf,sizeof(f));
    for(int i=1;i<=n;i++){
        f[i][i]=0;
        sum[i]=sum[i-1]+a[i];
    }
    //区间dp，首先确定区间长度len,区间长度便是动态规划的阶段。
    for(int len=2;len<=n;len++){     //阶段：区间长度len
        for(int L=1;L<=n-len+1;L++){ //状态左端点L取值范围
            int r=L+len-1;           //区间右端点r
            for(int k=L;k<r;k++){ //决策,连接点k
                f[L][r]=min(f[L][k]+f[k+1][r]+sum[r]-sum[L-1],f[L][r]);
            }
        }
    }
    cout<<f[1][n]<<endl;
    return 0;
}