题目传送门
题意: 给你两个长度为n(n<=1e5)的排列,求两个排列的最长公共子序列。
朴素LCS(最长公共子序列): 对于时间复杂度为O(n^2)的朴素lcs,我们用二维数组f[i][j]表示第一个序列的前i个和第二个序列的前j个,他们的最长公共子序列的长度,那么状态转移方程:a[i]==b[j]:dp[i][j]=dp[i-1][j-1]+1 else:dp[i][j]=max(dp[i-1][j],dp[i][j-1]) 。
但是这题很明显,时间和空间复杂度都不符合要求,但是这题是对于两个长度相同的排列,我们只需要在把第二个排列的每个数转化为它在另一个排列中的位置,再求第二个排列的最长上升子序列(O(nlogn)可求)就行。(自己想想,挺好懂的 )
朴素方法:
int a[N],b[N],dp[N][N];
signed main()
{
int n;
cin>>n;
for(int i=1;i<=n;i++)
cin>>a[i];
for(int i=1;i<=n;i++)
cin>>b[i];
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
if(a[i]==b[j])
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;
}
本题AC代码:
#include<bits/stdc++.h>
#define endl '\n'
#define null NULL
#define ls p<<1
#define rs p<<1|1
#define fi first
#define se second
#define mp make_pair
#define pb push_back
#define ll long long
#define int long long
#define vi vector<int>
#define mii map<int,int>
#define pii pair<int,int>
#define ull unsigned long long
#define all(x) x.begin(),x.end()
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define ct cerr<<"Time elapsed:"<<1.0*clock()/CLOCKS_PER_SEC<<"s.\n";
char *fs,*ft,buf[1<<20];
#define gc() (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<20,stdin),fs==ft))?0:*fs++;
inline int read(){int x=0,f=1;char ch=gc();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=gc();}
return x*f;}
using namespace std;
const int N=1e5+5;
const int inf=0x7fffffff;
const int mod=998244353;
const double eps=1e-6;
const double PI=acos(-1);
int a[N],b[N],c[N],d[N];
signed main()
{
int n;
cin>>n;
for(int i=1;i<=n;i++)
{
cin>>a[i];
c[a[i]]=i;
}
for(int i=1;i<=n;i++)
{
cin>>b[i];
b[i]=c[b[i]];
}
int len=1;
d[1]=b[1];
for(int i=1;i<=n;i++)
{
if(b[i]>d[len])
d[++len]=b[i];
else
d[lower_bound(d+1,d+len+1,b[i])-d]=b[i];
}
cout<<len<<endl;
}