动态规划（二）典型题目总结

最长上升子序列

f[i]记录以a[i]为结尾的最长上升子序列长度，初始条件是f[i] = 1；

递推式：

如果（a[i] > a[j] &&& f[j] + 1 > f[i])则更新f[i] = f[j] + 1;

更新过程

此题的递推式由小推大（最优子结构），由过去推现在（无后效性）
代码：

#include<bits/stdc++.h>
using namespace std;
const int N = 100;
int a[N];
int f[N];
int main(){
	int n;
	int ans = 1;
	cin >> n;
	for(int i = 1; i <= n; i++){
		cin >> a[i];
		f[i] = 1;
	}
	for(int i = 2; i <= n; i++){
		for(int j = 1; j < i; j++){
			if(a[i] > a[j]){
				f[i] = max(f[j] + 1,f[i]);
				ans = max(ans,f[i]);
			}
		}
	}
	cout << ans;
}

最长公共子序列（LCS)

题目描述：

给定两个字符串，输出其最长公共子序列的长度。
输入：ADABEC和DBDCA
输出：3 解释：最长公共子序列是DBC,其长度为3.（子序列可以是不连续的）

用f[i][j]记录序列a[1…i]和b[1…j]的最长公共子序列长度。
然后考虑末尾元素a[i]与b[j]是否在公共子序列中：

递归版本

#include<bits/stdc++.h>
#include<iostream>
using namespace std;
const int MAXN = 1000 + 10;
char str1[MAXN],str2[MAXN];
int dp[MAXN][MAXN];
int p[MAXN][MAXN]; 
void LCS_print(int i,int j){
	if(i == 0|| j == 0) return;
	if(str1[i - 1] == str2[j -1]){
		LCS_print(i - 1,j - 1); //递归 
		cout << str1[i - 1]; //回溯 
	}else if(dp[i][j - 1] > dp[i - 1][j]){
		LCS_print(i,j - 1);
	}else{
		LCS_print(i - 1,j);
	}
}
int main(){
	int n,m;
	gets(str1);
	gets(str2);
	n = strlen(str1);
	m = strlen(str2);
	for(int i = 1; i <= n; i++){
		for(int j = 1; j <= m; j++){
			if(str1[i - 1] == str2[j - 1]){
				dp[i][j] = dp[i - 1][j - 1] + 1;
				p[i][j] = 1;
			}else{
				dp[i][j] = max(dp[i - 1][j],dp[i][j - 1]);
			}
		}
	} 
	cout << dp[n][m] << endl;
	LCS_print(n,m); 
} 

非递归数组版

#include<bits/stdc++.h>
#include<iostream>
using namespace std;
const int MAXN = 1000 + 10;
char str1[MAXN],str2[MAXN];
int dp[MAXN][MAXN];
int p[MAXN][MAXN]; 
char s[MAXN];
/*
void LCS_print(int i,int j){
	if(i == 0|| j == 0) return;
	if(str1[i - 1] == str2[j -1]){
		LCS_print(i - 1,j - 1); //递归 
		cout << str1[i - 1]; //回溯 
	}else if(dp[i][j - 1] > dp[i - 1][j]){
		LCS_print(i,j - 1);
	}else{
		LCS_print(i - 1,j); 
	}
} 
*/
int main(){
	int n,m;
	gets(str1);
	gets(str2);
	n = strlen(str1);
	m = strlen(str2);
	for(int i = 1; i <= n; i++){
		for(int j = 1; j <= m; j++){
			if(str1[i - 1] == str2[j - 1]){
				dp[i][j] = dp[i - 1][j - 1] + 1;
				p[i][j] = 1;
			}else if(dp[i - 1][j] > dp[i][j - 1]){
				dp[i][j] = dp[i - 1][j];
				p[i][j] = 2;
			}else{
				dp[i][j] = dp[i][j - 1];
				p[i][j] = 3;
			}
		}
	} 
	cout << dp[n][m] << endl;
	int i = n,j = m,k = dp[n][m];
	while(i > 0 && j > 0){
		if(p[i][j] == 1){
			s[k--] = str1[i - 1];
			i--;j--;
		}else if(p[i][j] == 2){
			i--;
		}else{
			j--;
		}
	}
	for(int i = 1; i <= dp[n][m]; i++){
		cout << s[i];
	}
//	LCS_print(n,m); 
} 

最长公共字串

#include<bits/stdc++.h>
using namespace std;
const int MAXN = 1010;
char str1[MAXN];
char str2[MAXN];
int f[MAXN][MAXN];
int n,m,xindex,yindex;
char s[MAXN];
int maxn = 0;
int main(){
	gets(str1);
	gets(str2);
	n = strlen(str1);
	m = strlen(str2);
	for(int i = 1; i <= n; i++){
		for(int j = 1; j <= m; j++){
			if(str1[i -1] == str2[j -1]){
				f[i][j] = f[i - 1][j - 1] + 1;
			}else{
				f[i][j] = 0;
			}
			if(f[i][j] > maxn){
				maxn = f[i][j];
				 xindex = i;
				 yindex = j;
			}
			
		}
	}
	cout << maxn << endl;int mm = maxn;
	for(int i = xindex; i > 0; i--){
		s[maxn--] = str1[i - 1];
	}
	for(int i = 1; i <= mm; i++){
		cout << s[i];
	}
}

编辑距离

字符串编辑距离就是从一个字符修改到另一个字符时，编辑单个字符（比如修改、插入、删除）所需要的最少次数。
从字符串"NOTV"修改为字符串"LOVER"需4次单字符编辑操作：
LOVT(N改为L)
LOV(删除T)
LOVE(插入E）
LOVER(插入R)
因此编辑距离为4

问题：

输入两个字符串a、b，输入从字符串a修改到字符串b时的编辑距离。

输入：

a = “NOTV”,b = “LOVER”

输出：

4