LCS:Longest Common Subsequence,最长公共子序列
话不多说,上题:
P1439【模板】最长公共子序列
传送门
首先,我看到这个洛谷的题,我就想到了vjudgeOJ上的一道题目Common Subsequence,感觉题目很像,我就去看了看原题,结果发现都是求最长公共子序列,就兴冲冲的用相同的算法解题,交了一发,原地自闭。
第一发:
#include<bits/stdc++.h>
using namespace std;
int dp[10005][10005];
int a[10005],b[10005];
int main()
{
int n;
cin>>n;
for(int i=0;i<n;++i) cin>>a[i];
for(int i=0;i<n;++i) cin>>b[i];
memset(dp,0,sizeof(dp));
for(int i=1; i<=n; ++i)
{
for(int j=1; j<=n; ++j)
{
if(a[i-1]==b[j-1])
dp[i][j]=dp[i-1][j-1]+1;
else
dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
}
}
cout<<dp[n][n]<<endl;
return 0;
}
经过分析,我才发现,这道题MLE是因为,没看清题意数组开小了,但是我当我把数组开到题目规定的数组大小
【数据规模】
对于50%的数据,n≤1000
对于100%的数据,n≤100000
运行程序发现出错,二维数组是不能开到100000。
而且时间复杂度为T(n^2),直接超时了。
所以应该采取一种T(nlogn)的方法,经过无数次的失败,于是就有了第二发。
第二发:
#include <bits/stdc++.h>
using namespace std;
int a[100005];
int p[100005];
int low[100005];
int bf(int r,int x)
{
int l=1,mid;
while(l<=r)
{
mid=l+r>>1;
if(p[low[mid]]>p[x])
r=mid-1;
else
l=mid+1;
}
return l;
}
int main()
{
int n;
cin>>n;
for(int i=1;i<=n;++i) cin>>a[i];
for(int i=1;i<=n;++i)
{
int x;
cin>>x;
p[x]=i;
}
int ans=0;
for(int i=1; i<=n; ++i)
{
if(p[a[i]]>p[low[ans]])
low[++ans]=a[i];
else
low[bf(ans,a[i])]=a[i];
}
cout<<ans<<endl;
return 0;
}
这种方法思路核心跟LIS:Longest Increasing Subsequence(最长递增子序列)类似,即二分+插入。
就是把第二个序列离散化,其作用在于用第一个序列快速确定前后位置关系,即满足第二个元素在第二个序列离散化后位置大于第一个,那么就是两个序列的共同元素,并且位置关系也是正确的子序列位置关系。
至于二分+插入的作用是:如果下一个 元素不满足这关系,就给它在low中安排一个恰当的位置,但值得注意的是ans不会变化。