题意:
给出a,b,c三个字符串,问ab两个字符串最长不包含c的公共子序列。
解题思路:
基本的思想还是求最长共子序列的dp。
但是需要扩展一维,第三维k表示dp[i][j]最长公共子序列的后缀和c的前缀匹配到的位置的下一个位置,每次a[i]==b[j]的时候, 我们去找一下a[i],加上a[i]的子序列和c匹配到的位置,这里利用下kmp就可以,从k开始跑。
然后用一个结构体记录状态,这样比较好想。
代码:
#include <bits/stdc++.h>
using namespace std;
const int inf=0x3f3f3f3f;
struct p
{
int x, y, z;
}pre[105][105][105], e;
char a[105];
char b[105];
char c[105];
int dp[105][105][105];
int nex[105];
void get_next(const char P[], int nex[])
{
int q, k;
int m=strlen(P);
nex[0]=0;
for(q=1, k=0; q<m; ++q)
{
while(k>0 && P[q]!=P[k])
{
k=nex[k-1];
}
if(P[q]==P[k])
{
k++;
}
nex[q]=k;
}
return;
}
void dfs(p t)
{
int x, y, z;
x=t.x, y=t.y, z=t.z;
if(x==0 && y==0 && z==0)return;
// printf("%d %d %d %d\n", x, y, z, dp[x][y][z]);
p q=pre[x][y][z];
int i=q.x, j=q.y, k=q.z;
dfs(q);
if(dp[i][j][k]+1==dp[x][y][z])printf("%c", a[i]);
}
int main()
{
int i, j, k;
scanf("%s%s%s", a, b, c);
get_next(c, nex);
memset(dp, 0, sizeof dp);
int n=strlen(a), m=strlen(b), q=strlen(c);
for(i=0; i<n; i++)
{
for(j=0; j<m; j++)dp[i][j][0]=0;
}
for(i=0; i<=n; i++)
{
for(j=0; j<=m; j++)
{
for(k=0; k<q; k++)
{
if(dp[i+1][j][k]<dp[i][j][k])
{
dp[i+1][j][k]=dp[i][j][k];
pre[i+1][j][k].x=i;
pre[i+1][j][k].y=j;
pre[i+1][j][k].z=k;
}
if(dp[i][j+1][k]<dp[i][j][k])
{
dp[i][j+1][k]=dp[i][j][k];
pre[i][j+1][k].x=i;
pre[i][j+1][k].y=j;
pre[i][j+1][k].z=k;
}
if( a[i]==b[j])
{
int w=k;
while(w>0 && a[i]!=c[w])w=nex[w-1];
if(a[i]==c[w])w++;
else w=0;
if(dp[i+1][j+1][w]<dp[i][j][k]+1)
{
dp[i+1][j+1][w]=dp[i][j][k]+1;
pre[i+1][j+1][w].x=i;
pre[i+1][j+1][w].y=j;
pre[i+1][j+1][w].z=k;
// printf("%d %d %d %d %c %d\n", i+1, j+1, w, k, a[i], dp[i+1][j+1][w]);
}
}
}
}
}
int ma=0;
for(k=0; k<q; k++)
{
if(dp[n][m][k]>ma)
{
ma=dp[n][m][k];
e.x=n;
e.y=m;
e.z=k;
}
}
// printf("e%d %d %d %d\n", e.x, e.y, e.z, dp[e.x][e.y][e.z]);
if(ma==0)return 0*printf("0\n");
dfs(e);
}