4275: [ONTAK2015]Badania naukowe DP

来源：互联网发布：淘宝网李宁运动服编辑：程序博客网时间：2024/06/16 16:37

考虑把C串放到A,B串中，如果匹配成功，则从两个端点向两边做LCS，就是答案了。
详细的说一下：
用fi,j表示A[1..i],B[1..j]的LCS，用gi,j表示A[i..n],B[j..m]的LCS。
用di表示A串中从i开始匹配成功的最后端点。
用ei表示B串中从i开始匹配成功的最后端点。
ans=max(f[i−1][j−1]+f[d[i]+1][e[i]+1]+lenc)。

#include<iostream>#include<cstdio>using namespace std;int n,m,l,ans=-1;int a[3005],b[3005],c[3005],d[3005],e[3005];int f[3005][3005],g[3005][3005];inline int read(){    int a=0,f=1; char c=getchar();    while (c<'0'||c>'9') {if (c=='-') f=-1; c=getchar();}    while (c>='0'&&c<='9') {a=a*10+c-'0'; c=getchar();}    return a*f;}int main(){    n=read();    for (int i=1;i<=n;i++) a[i]=read();    m=read();    for (int i=1;i<=m;i++) b[i]=read();    l=read();    for (int i=1;i<=l;i++) c[i]=read();    for (int i=1;i<=n;i++)          for (int j=1;j<=m;j++)            if (a[i]==b[j]) f[i][j]=f[i-1][j-1]+1;            else f[i][j]=max(f[i][j-1],f[i-1][j]);    for (int i=n;i;i--)        for (int j=m;j;j--)            if (a[i]==b[j]) g[i][j]=g[i+1][j+1]+1;            else g[i][j]=max(g[i][j+1],g[i+1][j]);    if (!l) {cout << f[n][m] << endl; return 0;}    for (int i=1;i<=n;i++)        for (int j=i,now=1;j<=n;j++)        {            if (a[j]==c[now]) now++;            if (now>l) {d[i]=j; break;}        }    for (int i=1;i<=m;i++)        for (int j=i,now=1;j<=m;j++)        {            if (b[j]==c[now]) now++;            if (now>l) {e[i]=j; break;}        }    for (int i=1;i<=n;i++)         for (int j=1;j<=m;j++)            if (d[i]&&e[j]) ans=max(ans,f[i-1][j-1]+g[d[i]+1][e[j]+1]+l);    cout << ans << endl;    return 0;}

0 0