【模板】最长公共子序列

主要是个人讲解

传送门（洛谷）

这道题两种方法

一 ： 动态规划（但是此题规模过大过不了）：

这是一个状态！！！！！（只需找到状态转移）

我们考虑最后一个匹配方案
第一：
我们可以先得出，如果最后一个是匹配的，那么他们前一个也一定是匹配的，得出转移方程：$dp[i][j]=dp[i-1][j-1]+1$
第二（最后一个不匹配）：
如果是A串最后一个不匹配，那么我们考虑A串的前一个和B串匹配，就是dp[i-1][j]
如果是B串最后一个不匹配，那么我们考虑B串的前一个和A串匹配，就是dp[i][j-1]
得出转移方程：$dp[i][j]=max(dp[i-1][j],dp[i][j-1])$

Code

#include<bits/stdc++.h>
#define ll long long
#define in __int128
#define rep(i,a,b) for(register int i=(a);i<=(b);i++)
#define don(i,a,b) for(register int i=(a);i>=(b);i--)
using namespace std;
const int maxn=1e6+10;
const int maxm=1e3+10;
int n;
int a[100000+10],b[100000+10];
int dp[2010][2010];
 
void lcs(int i,int j)
{
    for(i=1; i<=n; i++) {
        for(j=1; j<=n; j++) {
            if(a[i-1] == b[j-1])
                dp[i][j] = dp[i-1][j-1] + 1;
            else dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
        }
    }
}
 
int main()
{
	freopen("input.txt","r",stdin);
    scanf("%d",&n);
    rep(i,1,n) scanf("%d",&a[i]);
    rep(i,1,n) scanf("%d",&b[i]);
    lcs(1,1);
    printf("%d\n", dp[n][n]);
    return 0;
}

二：转化为最长上升子序列

我们把一个串当作比对串，把第一个串处理一下，相当于离散化，就是用1,2,3…n
去处理A串，例如：A串:5,3,2,1,4 B串:3,5,4,2,1
5,3,2,1,4用1,2,3,4,5处理后对应关系就为5-1,3-2,2-3,1-4,4-5
那么B串成了2,1,5,3,4
因为A串处理后的顺序是1,2,3…n 上升的，那么我们就找B串的最长上升子序列（最好手工模拟一下）
可以这样看

红色连线部分就恰好匹配

Code

#include<bits/stdc++.h>
#define ll long long
#define in __int128
#define rep(i,a,b) for(register int i=(a);i<=(b);i++)
#define don(i,a,b) for(register int i=(a);i>=(b);i--)
using namespace std;
const int maxn=1e6+10;
const int maxm=1e3+10;
int n;
int a[maxn],b[maxn],san[maxn],len[maxn];
int d[maxn];

template <class t> inline void read(t &x)
{
    x=0;int f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)){x=10*x+ch-'0';ch=getchar();}
    x*=f;
}

template <class t> inline void write(t x)
{
    if(x<0){putchar('-');x=~x+1;}
    if(x>9) write(x/10);
    putchar(x%10+48);
}

/*--------------------------------------------------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------------------------------------------------*/

void readdata()
{
	read(n);
	rep(i,1,n) {
		read(a[i]);
		san[a[i]]=i;
	}
	rep(i,1,n) read(b[i]);
}

void work()
{
	/*int lenth=0;
	for(int i=1;i<=n;i++)
    {
        if(san[b[i]]>len[lenth])
        {
        	lenth++;
            len[lenth]=san[b[i]];
            continue;
        }
        int k=lower_bound(len,len+lenth,san[b[i]])-len;
        len[b[i]]=san[k];
    }*/
//----------------------------------------------------------------------
    int l=-1;//找最长上升子序列的过程（不懂上升子序列自己看）
	rep(i,1,n) {
		int j=lower_bound(d,d+l,san[b[i]])-d;
		d[i]=j+1;
		if(l<d[i]) l=d[i];
		d[j]=san[b[i]];
	}
	write(l);
    //printf("%d\n",lenth);
}

int main()
{
	freopen("input.txt","r",stdin);
	readdata();
	work();
	return 0;
}