区间DP入门题 51Nod 1021 石子合并

题目：http://www.51nod.com/Challenge/Problem.html#!#problemId=1021

题意：给出n堆石子，要把这n堆石子合并成1堆，每次只能选相邻的两个合并成一堆，每次合并的代价是这两堆石子的总数之和，问最少代价是多少。

解析：一个区间的最小合并代价是由这个区间内的子区间合并的最小值得到，并且要加上他们之前所花费的代价sum[j] - sum[i-1]。

代码：

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <string>
#include <vector>
#include <set>
using namespace std;

const int inf = 0x3f3f3f3f;
#define LL          long long
#define ll          LL
#define zero(a)     memset(a,0,sizeof(a))
#define minf(a)     memset(a,inf,sizeof(a))
#define rep(i,a,b)  for(int (i)=(a); (i)<=(b); ++(i))
#define rd(a)       scanf("%d",&(a))
#define rld(a)       scanf("%lld",&(a))
#define pf(a)       printf("%d\n",(a))

const int N = 105;
int dp[N][N];
int a[N];
int sum[N]={0};

int main()
{
    int n;
    rd(n);
    rep(i, 1, n)    rd(a[i]), sum[i] = sum[i-1]+a[i];
    rep(len, 2, n){
        rep(i, 1, n) {
            int j = i + len - 1;
            if (j > n) break;
            dp[i][j] = inf;
            rep(k, i, j - 1) {
                dp[i][j] = min(dp[i][j], dp[i][k] + dp[k + 1][j]);
            }
            dp[i][j] += sum[j] - sum[i-1];
        }
    }
    pf(dp[1][n]);
    return 0;
}

---

区间dp通用模板：

for(i = 1;i <= n;i++)
    dp[i][i] = 初始值；//一般为0
for(len = 2;len <= n;len++){//len选择区间长度
    for(i = 1;i <= n;i++){//枚举起点
        j = i + len - 1;//合并终点
        if(j > n)break;//不可越界
        for(k = i;k < j;k++)//枚举分割点，寻找最优分割
            dp[i][j] = max(dp[i][j], dp[i][k] + dp[k + 1][j] + w[i][j]);//状态转移
    }
}