校内模拟-变形lcs

emmmm a,b2个字符串，a串<=1e6，b串<=1e3，求他们的lcs。

那时间复杂度显然是有一些问题的。。。但是我们会发现lcs一定是<=1e3的，所以用f[i][j]表示b的前i个里面lcs长度是j。

f[i][j]就可以用我们刚才的类似的方式转移了emmm。想法就是考虑第i个字母是否可以作为最长公共子序列的一员。

最后输出就是找到i最大的合法的f[m][i]，输出i就可以了

#include <bits/stdc++.h>
using namespace std;

const int inf=0x3f3f3f3f;
const int maxn=100006;
const int maxm=1006;

int dp[maxm][maxm];
int lena,lenb;
char s1[maxn],s2[maxm];
int nxt[27][maxn];
int wle;

int main(){
	freopen("a.in","r",stdin);
	freopen("a.out","w",stdout);
	scanf("%s",s1+1);
	scanf("%s",s2+1);
	lena=strlen(s1+1);
	lenb=strlen(s2+1);
	memset(dp,inf,sizeof dp);
	memset(nxt,inf,sizeof nxt);
	for (int i=0;i<26;++i){
		for (int j=lena;j>=0;j--){
			nxt[i][j]=nxt[i][j+1];
			if (s1[j]-'a'==i)nxt[i][j]=j;
		}
	}
//	for (int i=0;i<26;++i){
//		printf("\n%c:",(char)('a'+i));
//		for (int j=0;j<=lena;++j){
//			printf("%d ",nxt[i][j]);
//		}
//	}
//	puts("-1-1-");
	for (int i=0;i<=lena;++i)dp[i][0]=0;
	for (int i=1;i<=lenb;++i)
	{
		for (int j=1;j<=i;++j)
		{
//			printf("%d %d\n",i,j);
//			system("PAUSE");
			dp[i][j]=dp[i-1][j];
			if (dp[i-1][j-1]<=lena) dp[i][j]=min(dp[i][j],nxt[s2[i]-'a'][dp[i-1][j-1]+1]);
		}
	}
//	puts("010101");
	for (int i=0;i<=lenb;++i){
		if (dp[lenb][i]==inf)break;
		wle=i;
	}
	printf("%d\n",wle);
	return 0;
}