所谓区间dp,顾名思义就是在一段区间上的动态规划。它既要满足dp问题的最优子结构和无后效性外,还应该符合在区间上操作的特点。我的理解是往往会对区间进行合并操作。抑或是单个元素(可看成一个小区间)跨区间进行操作。例如括号匹配问题,石子合并问题(通过多次的相邻合并,最后实质上会产生跨区间的合并,如果你把其中的石子看作参考系的话就很容易感觉出来),还有在整数中插入运算符号的问题(利用运算符的优先级以及交换律可看出)
这样以来,如果我们要得知一个大区间的情况,由于它必定是由从多个长度不一的小区间转移而来(转移情况未知),我们可以通过求得多个小区间的情况,从而合并信息,得到大区间。
对于一个长度为n的区间,确定它的子区间需要首尾两个指针,显然子区间数量级为n2,那区间dp的复杂度也就为n2
for (int len = 1; len < n; len++) { //操作区间的长度
for (int i = 0, j = len; j <= n; i++, j++) { //始末
//检查是否匹配(非必须)
for (int s = i; s < j; s++) {
//update
}
}
}
很多个石子合并成一堆,
最后一定会进行到还有两堆石子
然后把他们合并成一堆(这一步的花费是这一堆石子的总大小),再一个,最后这两堆石子一定(1~k)个和(k+1~n)的石子合并而来,所以我们就有了一个最后一步的公式:
dp[1][n] = min(dp[1][n] , dp[1][k] + dp[k+1][n] + sum[1][n])
然后,我们再看,dp[1][k]和dp[k+1][n]是不是也是通过上面的那个公式算出来的,而且要使得dp[1][n]最小那么,dp[1][k]和dp[k+1][n]一定也是如此算出来的最小值(也就是经过合并成这一堆的最小花费)
因此,我们要算出个大区间的石子合成最小值,是不是就要算出某两个小区间的最小值,如此下去,我们是不是就需要从区间长度的从小到大这样dp下去(因为我们只知道长度为1的区间值啊),最后算出这个dp[1][n]的最小值,所以我们就有了下面的代码:(有的同学可能有疑问,为什么要算出中间一些区间的dp呢,你看哈,我们的dp[1][k]是不是由某个dp[1][s]和dp[s+1][k]算出来的,那么我们这个dp[s+1][k]是不是就是中间的某个区间了)
#include<stdio.h>
#include<iostream>
#include<math.h>
#include<algorithm>
#define N 210
const int INF = 1e9;
using namespace std;
int dp[N][N];
int sum[N];
int a[N];
int main() {
int n;
while (~scanf("%d", &n)) {
for (int i = 1; i <= n; i++) {
scanf("%d", &a[i]);
sum[i] = sum[i - 1] + a[i];
}
for (int len = 1; len < n; len++) { //操作区间的长度
for (int i = 1, j = len + 1; j <= n; i++, j++) { //始末
//检查是否匹配(非必须)
dp[i][j] = INF;
for (int k = i; k < j; k++) {
dp[i][j] = min(dp[i][j], dp[i][k] + dp[k + 1][j] + sum[j] - sum[i - 1]);
}
}
}
printf("%d\n", dp[1][n]);
}
return 0;
}
好了,接下来我们来看看另一个经典的区间dp问题,括号匹配问题,这里我们就不讲题目的大意了,直接来讲题目的思路吧
我们要求出一个区间最大的括号匹配长度,同样,这个区间一定可以分成两个区间,但是,这里的分法就有了两种哦,没错,就是模板中检查是否匹配的那种啦!我们来举一个例子:(()(())),我们来拆这个字符串,因为是匹配问题,我们是不是有两种拆法,一种是类似于上面石子并归拆成两个区间,另一种就是最左边和最右边的匹配,然后加上中间的,下面就来看代码理解下吧!
#include<stdio.h>
#include<string.h>
#include<iostream>
#include<math.h>
#include<algorithm>
#define N 110
const int INF = 1e9;
using namespace std;
char str[N];
int dp[N][N];
bool ck(int i, int j) {
if ((str[i] == '(' && str[j] == ')') || (str[i] == '[' && str[j] == ']')) {
return true;
} else {
return false;
}
}
int main(int argc, const char * argv[]) {
while (~scanf("%s", str)) {
if (str[0] == 'e') break;
int len;
len = strlen(str);
memset(dp, 0, sizeof(dp));
for (int l = 1; l < len; l++) { //len = j - i 为当前区间长度
for (int i = 0, j = l; j < len; i++, j++) { // i++, j++
if (ck(i, j)) { // 匹配
dp[i][j] = dp[i + 1][j - 1] + 2;
}
// 讨论区间合并情况,求最大值
for (int k = i; k < j; k++) {
dp[i][j] = max(dp[i][j], dp[i][k] + dp[k + 1][j]);
}
}
}
printf("%d\n", dp[0][len - 1]);
}
return 0;
}