问题描述:
东东有两个序列A和B。
他想要知道序列A的LIS和序列AB的LCS的长度。
注意,LIS为严格递增的,即a1<a2<…<ak(ai<=1,000,000,000)。
Input:
第一行两个数n,m(1<=n<=5,000,1<=m<=5,000)
第二行n个数,表示序列A
第三行m个数,表示序列B
Output:
输出一行数据ans1和ans2,分别代表序列A的LIS和序列AB的LCS的长度
Simple Input
5 5
1 3 2 5 4
2 4 3 1 5
Simple Output
3 2
解题思路:
第一个子问题:
这里用到动态规划的思想,求某个位置作为区间右端点的最长上升子序列只需往前找所有值(Aj<Ai)比它小的f(j)的最大值,再加一就是f(i),即对于j<i且Aj<Ai,f[i]=max(f[j])+1。
最长上升子序列
⚫ 状态:定义 fi 表示以 Ai 为结尾的最长上升序列的方程。
⚫ 初始化:f1 = 1
⚫ 转移过程
⚫ 输出答案:max{f[i], i=1…n}
⚫ 时间复杂度:O(n^2)
第二个子问题:
最长公共子序列
⚫ 设计状态:假设 f[i][j] 为 A1, A2, …, Ai 和 B1, B2, …, Bj 的 LCS 长度
⚫ 初始化:初始 f[1][0] = f[0][1] = f[0][0] = 0
⚫ 转移方程:当 Ai == Bj 时,f[i][j] = f[i-1][j-1] + 1
⚫ 否则 f[i][j] = max(f[i-1][j], f[i][j-1])
⚫ 输出答案:f[n][m]
⚫ 时间复杂度:O(nm)
实验代码:
#include<iostream>
using namespace std;
int n,m,len1,len2,ansA,ansAB;
const int sz=5050;
long long int a[sz],b[sz],arrayA[sz],arrayB[sz];
int _inf=-(10e8+7);
int f[5001][5001];
void getAns_A()
{
a[0]=_inf;
for(int i=1;i<=5001;i++)
{
int j,maxx=0;
for(j=i-1;j>=0;j--)
{
if(a[j]<a[i])
{
if(arrayA[j]>maxx)
maxx=arrayA[j];
}
}
arrayA[i]=maxx+1;
}
for(int i=0;i<=5001;i++)
if(arrayA[i]>ansA)
ansA=arrayA[i];
}
void getAns_AB()
{
for(int i=0;i<=n;i++)
f[0][i]=0;
for(int i=0;i<=m;i++)
f[i][0]=0;
for(int i=1;i<=m;i++)
{
for(int j=1;j<=n;j++)
{
if(b[i]==a[j])
f[i][j]=f[i-1][j-1]+1;
else
f[i][j]=max(f[i-1][j], f[i][j-1]);
}
}
ansAB=f[m][n];
}
int main(void)
{
cin>>n>>m;
for(int i=1;i<=n;i++)
scanf("%d",&a[i]);
for(int i=1;i<=m;i++)
scanf("%d",&b[i]);
getAns_A();
cout<<ansA<<" ";
getAns_AB();
cout<<ansAB<<endl;
return 0;
}