【洛谷】P1439 【模板题】求最长公共子序列
1.题意
给出两个序列,如何求出其最长公共子序列?
2.分析
本题其实很简单,关键是 思维的转换 。如果使用最基本的 O(N*N)
算法,会因为N较大而TLE。那么有没有一种更优的方法呢?回想如何求最长不下降子序列的过程。我讲到过如何使用 O(NlogN)
的算法求出最长子序列,可以参看我的博客。
结合本题的数据,可以猜到这肯定是要使用O(NlogN)的算法。那么能不能用求最长不下降子序列的方法去求最长公共子序列呢?是可以的。
因为本题给出的是求指定范围内数字的最长公共子序列,也就是两个序列的值有明显的关系特点。当把第一个子序列映射成1-n的序列时,那么序列2按照上面的映射关系就可以相应的修改成另外一个序列。因为第一个序列已经是有序的了(1-n)。所以我们只需要求第二个子序列的最长不下降子序列即可。
举例如下:
序列a[]: 3 1 2
=> 对应映射成 1 2 3
【也就是3->1, 1->2, 2 -> 3 】
序列b[]: 1 2 3
=> 按照映射关系可以修改成序列:2 3 1
因为映射后的序列a已经是有序的了,所以只需要求出序列b映射后的最长不下降子序列即可。
3.代码
#include<iostream>
#include<algorithm>
#include<map>
using namespace std;
const int maxN = 1e3+5;
int n;
int a[maxN], b[maxN]; //输入两个数组的值
int tail[maxN];//tail[i]表示长度为i 的最长不下降子序列结尾中的最小值
map<int,int> tran;//做一个映射
void proc(){
int cnt = 1;
for(int i = 1;i<=n;i++){
tran[a[i]] = cnt;
cnt ++;
}
for(int i = 1;i<=n;i++){
b[i] = tran[b[i]];
}
}
int main(){
cin >> n;
for(int i = 1;i<= n;i++)
cin >> a[i];
for(int i = 1;i<= n;i++)
cin >> b[i];
//step1.先做一个映射与转换
proc();
//step2.求出数组b的最长不下降子序列
int cnt = 1;//最长不下降子序列的长度
tail[cnt] = b[1];
for(int i = 2; i<=n; i++){
if(b[i] >= tail[cnt])//比最大的都要大,那么就直接往后添加
tail[++cnt] = b[i];
else{
int idx = lower_bound(tail+1,tail+cnt+1,b[i]) - tail;//找出第一个大于且等于b[i]的值
tail[idx] = b[i];
}
}
cout << cnt <<"\n";
}
测试用例
5
1 2 3 4 5
3 2 1 4 5
5
3 2 1 4 5
1 2 3 4 5