区间dp

区间dp

石子合并

题目描述：有n堆石子排成一排，每堆石子有一定的数量。将n堆石子并成为一堆，每次只能合并相邻的两堆石子，合并的花费为这两堆石子的总数。经过n-1次合并后成为一堆，求总的最小花费。
输入：第一行是整数n，表示有n堆石子。第二行有n个数，分别表示这n堆石子的数目。
输出：总的最小花费。
输入样例：
3
2 4 5
输出样例：
17
提示：样例的计算过程是：第一次合并2+4=6；第二次合并6+5=11；总花费6+11=17。

for(int i=1; i<=n; i++) 
    dp[i][i] = 0;
for(int len=2; len<=n; len++)       //len：区间[i,j]的长度，从小区间扩展到大区间
	for(int i=1; i<=n-len+1; i++){  // 区间起点i
   		int j = i + len - 1;        // 区间终点j，i<j<=n
        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] + w[i][j]);
	}

四边形不等式优化

只需一个简单的优化操作，就能把上面代码的复杂度变为O(n2) O(n^2)O(n
2
)。这个操作就是把循环i≤k<j i ≤ k < ji≤k<j改为：
    s[i][j−1]≤k≤s[i+1][j] s[i][j-1] ≤ k ≤ s[i+1][j]s[i][j−1]≤k≤s[i+1][j]
  其中s[i][j] s[i][j]s[i][j]记录从i到j的最优分割点。在计算dp[i][j] dp[i][j]dp[i][j]的最小值时得到区间[i,j] [i, j][i,j]的分割点k kk，记录在s[i][j] s[i][j]s[i][j]中，用于下一次循环。

for(i = 1;i <= n;i++){
    dp[i][i] = 0;                        
    s[i][i] = i;                        //s[][]的初始值
}
for(int len = 2; len <= n; len++)       
	for(int i = 1; i <= n-len+1; i++){   
   		int j = i + len - 1;            
   		for(k = s[i][j - 1]; k <= s[i + 1][j]; k++){   //缩小循环范围
            if(dp[i][j] > dp[i][k] + dp[k + 1][j] + w[i][j]){  //是否更优
   			   dp[i][j] = dp[i][k] + dp[k + 1][j] + w[i][j];
               s[i][j] = k;                 //更新最佳分割点
            }
   		}
	}

字符串修改

题目描述：给定两个长度相等的字符串A、B，由小写字母组成。一次操作，允许把A中的一个连续子串（区间）都转换为某个字符（就像用刷子刷成一样的字符）。要把A转换为B，问最低操作数是多少?
输入：第一行是字符串A，第二行是字符串B。两个字符串的长度不大于100。
输出：一个表示答案的整数。
输入样例：
zzzzzfzzzzz
abcdefedcba
输出样例：
6
提示：第1次把zzzzzfzzzzz转换为aaaaaaaaaaa，第2次转为abbbbbbbbba，第3次转为abccccccccba…

算法思路：
（1）如果原串与结果的字符都不同时，等同于在空白字符串基础上修改得到结果。
所以首先考虑空白字符串，dp[i][j]代表从i到j区间内需要修改的最小步数。
1.如果B[i]==B[j]，此时dp [i][j] = dp [i+1][j]或dp [i][j] = dp [i][j-1]。
2.如果B[i]!=B[j]，此时类似石子合并
dp [i][j] = min( dp [i][k] + dp [k+1][j] ) (i <= k < j )
3.递推边界：修改单个字符需要一次操作。
dp [i][i] = 1

（2）当原串有与结果相同的字符时，可以考虑保留。
所以要在原来的dp数组上进行修改。
1.如果A[i] ==B [i]，此时有dp [0][i] = dp [0][i-1]。(i=0时dp[0][0]=0)
2.如果A[i] != B[i],此时有
dp[0][i]=min( dp [0][k] + dp [k+1][j] )(i <= k < j)

#include<iostream>
#include<string>
#include<algorithm>
using namespace std;

string A, B;
int dp[100][100];

int main()
{
	cin >> A >> B;
	if (A.length() != B.length())
		return 0;

	int n = A.length();

	for (int i = 0; i < n; i++)
	{
		dp[i][i] = 1;
	}

	for (int len = 2; len <= n; len++)
	{
		for (int i = 0; i <= n - len; i++)
		{
			int j = i + len - 1;
			
			dp[i][j] = 100000;
			if (B[i] == B[j])
			{
				dp[i][j] = dp[i][j - 1];
			}
			else
			{
				for (int k = i; k < j; k++)
				{
					dp[i][j] = min(dp[i][j], dp[i][k] + dp[k + 1][j]);
				}
			}

		}
	}

	for (int i = 0; i < n; i++)
	{
		if (!i&&A[i] == B[i])
		{
			dp[0][i] = 0;
		}
		if (i&&A[i] == B[i])
		{
			dp[0][i] = dp[0][i - 1];
		}
		else
		{
			for (int k = 0; k < i; k++)
			{
				dp[0][i] = min(dp[0][i], dp[0][k] + dp[k + 1][i]);
			}
		}
	}

	cout << dp[0][n - 1] << endl;

	system("pause");
	return 0;
}