[AHOI2006]基因匹配

https://www.zybuluo.com/ysner/note/1332445

题面

给出两个数列，它们都包含(n)种元素，并且每种元素都恰有(5)个，求它们的最长公共子序列。

• (nleq20000)

解析

一般的(LCS)求法是(O(n^2))的。
设(f[i][j])表示在(s1)中匹配到第(i)个，(s2)中匹配到第(j)个。

fp(i,1,n)
    fp(j,1,n)
    {
      f[i][j]=max(f[i][j-1],f[i-1][j]);
      if(a[i]==b[j]) f[i][j]=max(f[i][j],f[i-1][j-1]+1);
    }

但是(LCS)有时可以转化为(LIS)来做。
假设有两个序列
(s1[6]={a,b,c,a,d,c},s2[7]={c,a,b,e,d,a,b})。
记录(s1)中每个元素在(s2)中出现的位置, 再将位置按降序排列, 则上面的例子可表示为：
(loc(a)={6,2},loc(b)={7,3},loc(c)={1},loc(d)={5})
将(s1)中每个元素的位置按(s1)中元素的顺序排列成一个序列
(s3={6,2,7,3,1,6,2,5,1})
再对(s3)求(LIS)得到的值即为求(LCS)的答案。

感觉这有点像一种等效替代，倒序实际上保证了不能自己向自己转移。

这题转化为(LIS)问题后序列长度为(5*10^5)，用(O(nlogn))的(LCS)算法可以解决。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#define ll long long
#define re register
#define il inline
#define fp(i,a,b) for(re int i=a;i<=b;++i)
#define fq(i,a,b) for(re int i=a;i>=b;--i)
using namespace std;
const int N=5e5+100;
int n,m,pre[N][6],f[N],sta[N],t[N],top,ans;
il ll gi()
{
  re ll x=0,t=1;
  re char ch=getchar();
  while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
  if(ch=='-') t=-1,ch=getchar();
  while(ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar();
  return x*t;
}
il int max(re int x,re int y){return x>y?x:y;}
il void upd(re int x,re int w){for(;x<=n;x+=x&-x) t[x]=max(t[x],w);}
il int que(re int x){re int res=0;for(;x;x-=x&-x) res=max(res,t[x]);return res;}
int main()
{
  m=gi();n=m*5;
  fp(i,1,n)
    {
      re int x=gi();
      pre[x][++pre[x][0]]=i;
    }
  fp(i,1,n)
    {
      re int x=gi();
      fq(j,5,1) sta[++top]=pre[x][j];
    }
  fp(i,1,top)
    f[i]=que(sta[i]-1)+1,upd(sta[i],f[i]),ans=max(ans,f[i]);
  printf("%d
",ans);
  return 0;
}