【算法笔记】动态规划 线性DP

最长上升子序列（LIS）

一个数列S如果分别是已知数列的单调上升子序列，且是所有符合此条件序列中最长的，则S称为已知序列的最长上升子序列，求S。

拆分问题：求长度为n的数列的最长公共子序列可以转化成求其数列前n-1项的最长公共子序列，再判断是否要加第n项。

状态转移方程：if (a[i] > a[j]) dp[i] = MAX (dp[i], dp[j] + 1);

基础的最长公共子序列代码：

#include<iostream>
using namespace std;
int Max(int x,int y)
{
	if(x>y)
		return x;
	else
		return y;
}
//最长递增序列
int main()
{	
	int a[10] = {2,5,3,7,4,6,8,1,2,9};  //序列a
	int dp[10];  //存放最长递增子序列的长度
	int i,j;
	int result;  //用于存储最大长度
	for(i=0;i<10;i++)  //拆分问题
	{
		dp[i] = 1;  //元素本身算序列的一个子序列，所以初始化为1
		for(j=0;j<i;j++)
		{
			if(a[j]<a[i])
				dp[i] = Max(dp[i],dp[j]+1);  //状态转移方程
		}
		result = Max(result,dp[i]);
	}
	cout<<result<<endl;
	return 0;
}

这里有一个关于最长递增子序列的题目 NYOJ17

单调递增最长子序列

时间限制： 3000 ms  |  内存限制： 65535 KB

难度： 4

描述

求一个字符串的最长递增子序列的长度
如：dabdbf最长递增子序列就是abdf，长度为4

输入

第一行一个整数0<n<20,表示有n个字符串要处理
随后的n行，每行有一个字符串，该字符串的长度不会超过10000

输出

输出字符串的最长递增子序列的长度

样例输入

3
aaa
ababc
abklmncdefg

样例输出

1
3
7

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

char a[10001];
int dp[10001];
int result[21];  //用于存储最大长度

int Max(int x,int y)
{
	if(x>y)
		return x;
	else
		return y;
}

void lis(int n)
{
	result[n] = -1;
	int i,j;
	int num1,num2;
	for(i=0;a[i];i++)  //拆分问题
	{
		dp[i] = 1;  //元素本身算序列的一个子序列，所以初始化为1
		for(j=0;j<i;j++)
		{
			num1 = a[j];
			num2 = a[i];
			if(num1<num2)
				dp[i] = Max(dp[i],dp[j]+1);  //状态转移方程
		}
		result[n] = Max(result[n],dp[i]);
	}
}

//最长递增序列
int main()
{
	int n;
	cin>>n;
	int i,len;
	for(i=0;i<n;i++)
	{
		cin>>a;
		//len = strlen(a);
		lis(i);
	}
	for(i=0;i<n;i++)
		cout<<result[i]<<endl;
	return 0;
}

最长公共子序列（LCS）

给定的两个序列X和Y，当另一个序列Z既是X的子序列又是Y的子序列时，称Z是X和Y的公共子序列。两个序列的公共子序列有很多个，最长公共子序列是子序列中元素个数最多的。

拆分问题：序列X内元素为x1时与序列Y的最长公共子序列，序列X元素为x1和x2时与序列Y的最长公共子序列……序列X元素为x1、x2、…、xi时与序列Y的最长公共子序列。

状态转移方程： 

Eg：X={A,BC,B,D,A,B}  Y={B,D,C,A,B,A}

基础最长公共子序列代码：

#include<iostream>
using namespace std;

int Max(int x,int y)
{
	if(x>=y)
		return x;
	else
		return y;
}

int main()
{
	char x[]={' ','A','B','C','B','D','A','B'};
	char y[]={' ','B','D','C','A','B','A'};
	int dp[8][7];
	int i,j;
	for(i=0;i<7;i++)  //初始化
	{
		dp[0][i]=0;
		dp[i][0]=0;
	}
	dp[7][0] = 0;
	for(i=1;i<8;i++)
	{
		for(j=1;j<7;j++)
		{
			//状态转移方程
			if(x[i] == y[j])
				dp[i][j] = dp[i-1][j-1]+1;
			else
				dp[i][j] = Max(dp[i-1][j],dp[i][j-1]);  //默认如果dp[i-1][j]=dp[i][j-1]取dp[i-1][j]
		}
	}
/*
	//打印表格
	for(i=0;i<8;i++)
	{
		for(j=0;j<7;j++)
			cout<<dp[i][j]<<" ";
		cout<<endl;
	}
*/
	//返回去搜索最长公共子序列的元素
	i = 7;
	j = 6;
	int t = dp[i][j];
	char res[10];
	while(t>0 && i>=0 && j>=0)
	{
		if((dp[i][j] == dp[i-1][j-1]+1) && x[i]==y[j])
		{
			res[t] = x[i];
			i = i-1;
			j = j-1;
			t--;
		}
		else if(dp[i][j] == dp[i-1][j])
			i=i-1;
		else
			j=j-1;
	}
	for(t=1;t<=dp[7][6];t++)
		cout<<res[t]<<" ";
	cout<<endl;
	return 0;
}

这里有一个关于最长公共子序列的题目 NYOJ36

最长公共子序列

时间限制： 3000 ms  |  内存限制： 65535 KB

难度： 3

描述

咱们就不拐弯抹角了，如题，需要你做的就是写一个程序，得出最长公共子序列。
tip：最长公共子序列也称作最长公共子串(不要求连续)，英文缩写为LCS（Longest Common Subsequence）。其定义是，一个序列 S ，如果分别是两个或多个已知序列的子序列，且是所有符合此条件序列中最长的，则 S 称为已知序列的最长公共子序列。

输入

第一行给出一个整数N(0<N<100)表示待测数据组数
接下来每组数据两行，分别为待测的两组字符串。每个字符串长度不大于1000.

输出

每组测试数据输出一个整数，表示最长公共子序列长度。每组结果占一行。

样例输入

2
asdf
adfsd
123abc
abc123abc

样例输出

3
6

#include<iostream>
using namespace std;

char x[1001];
char y[1001];
int dp[1001][1001];
int res[101];

int Max(int x,int y)
{
	if(x>=y)
		return x;
	else
		return y;
}

void lcs(int n)
{
	int i,j;
	int xlen,ylen;
	for(i=0;x[i];i++)  //初始化
	{
		dp[i][0]=0;
	}
	xlen = i;
	for(j=0;y[j];j++)  //初始化
	{
		dp[0][j]=0;
	}
	ylen = j;
	for(i=1;i<=xlen;i++)
	{
		for(j=1;j<=ylen;j++)
		{
			//状态转移方程
			if(x[i-1] == y[j-1])
				dp[i][j] = dp[i-1][j-1]+1;
			else
				dp[i][j] = Max(dp[i-1][j],dp[i][j-1]);  //默认如果dp[i-1][j]=dp[i][j-1]取dp[i-1][j]
		}
	}
	res[n] = dp[xlen][ylen];
}

int main()
{
	int n;
	cin>>n;
	int i;
	for(i=0;i<n;i++)
	{
		cin>>x;
		cin>>y;
		lcs(i);
	}
	for(i=0;i<n;i++)
		cout<<res[i]<<endl;
	return 0;
}