最长公共子序列和其变形-CSDN博客

LCS是Longest Common Subsequence的缩写，即最长公共子序列。一个序列，如果是两个或多个已知序列的子序列，且是所有子序列中最长的，则为最长公共子序列。（摘自百度百科）

一.DP通法

　　给定串1和串2，求其LCS

　　譬如给定2个序列：

1 2 3 4 5 3 4 2 5

求其最长的公共子序列。

我们可以用

　　如果不相同，即无法更新公共元素，考虑继承：

$dp[ i ] [ j ] = max(dp[ i-1 ][ j ] , dp[ i ][ j-1 ]);$

如果当前的串

#include<iostream>
#include<cstdio>    
#include<cctype>
#include<cmath>
#define R register int
using namespace std;
const int N=5010;
int n,m,dp[N][N];
char c1[N],c2[N];

int main()
{
    scanf("%d%d",&n,&m);
    for(R i=1;i<=n;i++) scanf("%c",&c1[i]);
    for(R i=1;i<=m;i++) scanf("%c",&c2[i]);;
    for(R i=1;i<=n;i++)
        for(R j=1;j<=m;j++)
        {
            dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
            if(c1[i]==c2[j]) dp[i][j]=max(dp[i][j],dp[i-1][j-1]+1);
        }
    printf("%d\n",dp[n][m]);
    return 0;
}

二.可是总有一些独特的：洛谷 P1439 【模板】最长公共子序列

对于这个题而言，朴素算法是 $n^2的，会被10^5卡死，所以我们可以考虑nlogn的做法：$

因为两个序列都是

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cctype>
 4 #include<cmath>
 5 #define R register int
 6 using namespace std;
 7 const int N=100010;
 8 int f[N],a[N],b[N],c[N];
 9 
10 inline int g()
11 {
12     int ret=0,fix=1; char ch;
13     while(!isdigit(ch=getchar())) fix=ch=='-'?-1:fix;
14     do ret=(ret<<3)+(ret<<1)+ch-'0'; while(isdigit(ch=getchar()));
15     return ret*fix;
16 }
17 
18 int main()
19 {
20     R n=g();
21     for(R i=1;i<=n;i++) a[i]=g(),f[a[i]]=i;
22     for(R i=1;i<=n;i++) b[i]=g(),b[i]=f[b[i]];
23     R len=0;
24     for(R i=1;i<=n;i++)
25         if(b[i]>c[len]) c[++len]=b[i];
26         else
27         {
28             R pos=lower_bound(c+1,c+len+1,b[i])-c;
29             c[pos]=b[i];    
30         }
31     printf("%d\n",len);
32     return 0;
33 }