最长公共子序列（稀疏序列）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=4001;
const int maxn=0x3f;
void read(int &n)
{
    char c='+';int x=0;bool flag=0;
    while(c<'0'||c>'9'){c=getchar();if(c=='-')flag=1;}
    while(c>='0'&&c<='9')
    x=(x<<1)+(x<<3)+c-48,c=getchar();
    flag==1?n=-x:n=x;
}
int a[MAXN];
int b[MAXN];
int n;
int c[MAXN];// 预处理后的数组
int tot=0;
int dp[MAXN]; 
char aa[MAXN];
char bb[MAXN];
void work()
{
        memset(dp,0x3f,sizeof(dp));
        tot=0;
        memset(c,0,sizeof(c));
        memset(a,0,sizeof(a));
        memset(b,0,sizeof(b));
        for(int i=0;i<strlen(aa);i++)
            a[i+1]=aa[i]-48;
    
        for(int i=0;i<strlen(bb);i++)
            b[i+1]=bb[i]-48;

        for(int i=1;i<=strlen(bb);i++)
        {
            for(int j=strlen(aa);j>=1;j--)
                if(b[i]==a[j])
                    c[++tot]=j;
        }
        //for(int i=1;i<=tot;i++)
        //    printf("%d ",c[i]);
        for(int i=1;i<=tot;i++)
        {
            int p=upper_bound(dp+1,dp+tot+1,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=1;i<=tot+1;i++)
        {
            if(dp[i]>=maxn)
            {
                printf("%d
",i-1);
                break;
            }
        }
    }
    
    return 0;
}