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

如题…

---

【分析】
我们把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;
}