最长上升子序列(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
-
第一行一个整数0<n<20,表示有n个字符串要处理
#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
-
第一行给出一个整数N(0<N<100)表示待测数据组数
#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;
}