最长公共子序列（稀疏序列）nlogn解法

首先这种做法只能针对稀疏序列，

比如这种情况：

abc abacabc。

会输出5

，，，，就比较尴尬，

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<cmath>
 #include<algorithm>
 #define lli long long int
 using namespace std;
 const int MAXN=;
 const int maxn=0x3f;
 void read(int &n)
 {
     char c='+';int x=;bool flag=;
     while(c<''||c>''){c=getchar();if(c=='-')flag=;}
     while(c>=''&&c<='')
     x=(x<<)+(x<<)+c-,c=getchar();
     flag==?n=-x:n=x;
 }
 int a[MAXN];
 int b[MAXN];
 int n;
 int c[MAXN];// 预处理后的数组
 int tot=;
 int dp[MAXN];
 char aa[MAXN];
 char bb[MAXN];
 void work()
 {
         memset(dp,0x3f,sizeof(dp));
         tot=;
         memset(c,,sizeof(c));
         memset(a,,sizeof(a));
         memset(b,,sizeof(b));
         for(int i=;i<strlen(aa);i++)
             a[i+]=aa[i]-;

         for(int i=;i<strlen(bb);i++)
             b[i+]=bb[i]-;

         for(int i=;i<=strlen(bb);i++)
         {
             for(int j=strlen(aa);j>=;j--)
                 if(b[i]==a[j])
                     c[++tot]=j;
         }
         //for(int i=1;i<=tot;i++)
         //    printf("%d ",c[i]);
         for(int i=;i<=tot;i++)
         {
             int p=upper_bound(dp+,dp+tot+,c[i])-dp;
             if(dp[p]==c[i])
                 continue;
             dp[p]=c[i];
         }
 }
 int  main()
 {
     //ios::sync_with_stdio(false);
     while(~scanf("%s %s",aa,bb))
     {
         work();
         for(int i=;i<=tot+;i++)
         {
             if(dp[i]>=maxn)
             {
                 printf("%d\n",i-);
                 break;
             }
         }
     }

     return ;
 }