大体题意:
给你两个字符串,求LCS的长度,一个不超过100W,一个不超过1000
思路:
正常的LCS的dp肯定不能进行,这样会超时。
令dp[i][j]表示当前枚举的LCS的长度为第i位,字符为第二个字符串的的第j 个字符。dp[i][j] 为第一个字符串的的前哪个位置。
那么假设dp[i][j] 是已知的话,那么考虑转移第j+1个字符,假设对LCS没有贡献的话,只能转移 dp[i][j+1] = min(dp[i][j]);
否则有贡献的话,既能转移 DP[i][j+1] 也能转移 dp[i+1][j+1]
不要漏掉状态,否则会WA掉。
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#define Siz(x) (int)x.size()
#define Max(a,b) ((a)>(b)?(a):(b))
using namespace std;
const int inf = 0x3f3f3f3f;
char s1[1000000+7];
char s2[1007];
int T;
int dp[1007][1007];
vector<int>g[2][26];
vector<int>::iterator it;
char t[1000000+7];
int main(){
scanf("%d",&T);
while(T--){
scanf("%s%s",s1+1,s2+1);
int len1 = strlen(s1+1);
int len2 = strlen(s2+1);
if (len1 < len2){
swap(len1,len2);
strcpy(t+1,s1+1);
strcpy(s1+1,s2+1);
strcpy(s2+1,t+1);
}
for (int i = 0; i < 2; ++i){
for (int j = 0; j < 26; ++j) g[i][j].clear();
}
for (int i = 1; i <= len1; ++i){
int id=s1[i]-'a';
g[0][id].push_back(i);
}
memset(dp,inf,sizeof dp);
for (int i = 1; i <= len2; ++i){
int id=s2[i]-'a';
g[1][id].push_back(i);
}
bool ok = 0;
for (int j = 1; j <= len2; ++j){
int id = s2[j]-'a';
if (Siz(g[0][id]))dp[1][j] = g[0][id].front(),ok = 1;
}
int ans = ok;
// printf("%d\n",ans);
for (int i = 1; i < 1007; ++i){
for (int j = 1; j < len2; ++j){
if (dp[i][j] != inf){
int id = s2[j+1] - 'a';
if (Siz(g[0][id]))it = lower_bound(g[0][id].begin(),g[0][id].end(),dp[i][j]+1);
// printf("%c %d %d %d\n",id+'a',dp[i][j]+1,i,j);
if (!Siz(g[0][id]) || it == g[0][id].end() ){
dp[i][j+1] = min(dp[i][j],dp[i][j+1]);
ans = Max(ans,i);
}
else {
// printf("haha\n");
dp[i+1][j+1] = min(*it,dp[i+1][j+1]);
dp[i][j+1] = min(dp[i][j+1],dp[i][j]); /// 状态转移不要漏掉。。
ans = Max(ans,i+1);
}
}
}
}
printf("%d\n",ans);
}
return 0;
}
/**
111
abcdefghijklmnopqrstuvwxyz
bbddee
mmmi
immm
3
ijh
hjaij
2
**/