题目描述:
给定n与1~n的两个全排列a,b,请计算a与b的最大公共子序列长度
输入
第一行:整数n
第二行:序列a
第三行:序列b
输出
一行,最大公共子序列长度
样例输入
5
1 2 3 4 5
3 2 1 4 5
样例输出 3
数据范围:1≤n≤100000
由于数据范围巨大,传统的
O
(
n
2
)
O(n^2)
O(n2)算法肯定不行。但LCS问题不像LIS那样有优化方法,所以不能按照LCS的方法做。
由于给的是1~n的全排列,所以序列中所有元素都是相同的,不同的只是排列顺序罢了。而LCS是从前往后比较的,所以如果一个序列中的从左往右的某几个元素在另一个序列中出现的位置是单调递增的,那么这些元素就构成了LCS的一部分。
因此,可以用一个数组loc[]标记一个数组中的元素在另一个数组中出现的位置,这个数组的LIS就是原来两个序列的LCS。由于LIS存在
O
(
n
l
o
g
n
)
O(nlogn)
O(nlogn)的做法,所以能够通过。
#include <cstdio>
#include <algorithm>
#define INF 0x3f3f3f3f
const int M=100002;
using namespace std;
int dp[M]={-1},loc[M],a[M]; //loc[a[i]]标记a[i]在b[]中出现的位置
int main(){
int n;
scanf("%d",&n);
for(int i=0;i<n;i++) scanf("%d",&a[i]);
for(int i=0;i<n;i++){ //b数组没有必要设立
int b;
scanf("%d",&b);
loc[b]=i;
}
fill(dp+1,dp+n+1,INF);
int res=0;
for(int i=0; i<n; i++) {
if(loc[a[i]]>dp[res]) {
dp[++res]=loc[a[i]];
}
else {
int pos=lower_bound(dp,dp+res,loc[a[i]])-dp;
dp[pos]=loc[a[i]];
}
}
printf("%d",res);
return 0;
}
启示:如果此题不是全排列,但能保证两个序列的元素相同,只是位置不同,那么可以先离散化再按照上面的方法处理(当然也可以用map)。