[模板]求最长公共子序列

如题…

---

【分析】
我们把a序列中的每个元素在b中出现的位置保存起来，再按照降序排列，排列后再代入a的每个对应元素，那就转化为了求这个新的序列的最长上升子序列了。如：a[] = {a, b, c,} b[] = {a,b,c,b,a,d},那么a中的a，b，c在b中出现的位置分别就是{0,4},{1,3},{2}分别按降序排列后代入a序列就是{4,0,2,3,1},之所以要按照降序排列，目的就是为了让每个元素只取到一次。

接下来的问题就是要求最长升序子序列问题了，也就是求LIS

---

【代码】
[代码不是很标准，见谅]

//poj 3260 Common Subsequence#include<iostream>#include<cstdio>#include<cmath>#include<cstring>#include<algorithm>#define ll long long#define M(a) memset(a,0,sizeof a)#define fo(i,j,k) for(i=j;i<=k;i++)using namespace std;const int mxn=1e6;int tp,n,m;int b[1005][2005],c[10005],d[10005];char s1[1005],s2[1005];inline int LIS(){    int i,j,len=0;    fo(i,1,tp)    {        if(c[i]>d[len]) d[++len]=c[i];        else if(c[i]<d[1]) d[1]=c[i];        d[lower_bound(d+1,d+len+1,c[i])-d]=c[i];    }    return len;}int main(){    int i,j;    scanf("%d%d",&n,&m);    while(scanf("%s%s",s1+1,s2+1)==2)    {        M(b);M(d);tp=0;        int len1=strlen(s1+1),len2=strlen(s2+1);        fo(i,1,len2)          b[s2[i]][++b[s2[i]][0]]=i;        fo(i,1,len1)          for(j=b[s1[i]][0];j;j--)            c[++tp]=b[s1[i]][j];        printf("%d\n",LIS());    }    return 0;}