【bzoj4275】[ONTAK2015]Badania naukowe

Description

给定三个数字串A,B,C，请找到一个A,B的最长公共子序列，满足C是该子序列的子串。
Input

第一行包含一个正整数n(1<=n<=3000)，表示A串的长度。
第二行包含n个正整数，其中第i个数表示Ai。
第三行包含一个正整数m(1<=m<=3000)，表示B串的长度。
第四行包含m个正整数，其中第i个数表示Bi。
第五行包含一个整数k(0<=k<=3000)，表示C串的长度。
第六行包含k个正整数，其中第i个数表示Ci。
Output

输出一个整数，即满足条件的最长公共子序列的长度，如果无解输出-1。特别的，如果k为0且无解，请输出0。
Sample Input

7

1 2 2 3 1 1 2

6

1 2 1 3 1 2

2

3 2
Sample Output

4
HINT

找到的最长个公共子序列为(1,2,3,2)。

题解
记f[i][j]为A,B从前往后做最长公共子序列
g[i][j]为A,B从后往前做最长公共子序列
x[i]表示A从i位置开始最近到x[i]位置出现C子序列
y[i]表示B从i位置开始最近到y[i]位置出现C子序列
ans=max(f[i-1][j-1]+g[x[i]+1][y[i]+1]+|C|)

代码

#include<bits/stdc++.h>
typedef long long ll;
const ll mod=1000000007LL;
const int inf=1000000000;
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if (ch=='-')f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n,m,C,a[3005],b[3005],c[3005],A[3005],B[3005],ans=-1;
int f[3005][3005],g[3005][3005];
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();
    C=read();for (int i=1;i<=C;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-1][j],f[i][j-1]);
    if (!C) return printf("%d",f[n][m]),0;
    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+1][j],g[i][j+1]);
    for(int i=1;i<=n;i++)
        for(int j=i,p=1;j<=n;j++)
        {
            if(a[j]==c[p])p++;
            if(p>C){A[i]=j;break;}
        }
    for(int i=1;i<=m;i++)
        for(int j=i,p=1;j<=m;j++)
        {
            if(b[j]==c[p])p++;
            if(p>C){B[i]=j;break;}
        }
    for (int i=1;i<=n;i++)
        if (A[i]) for (int j=1;j<=m;j++)
            if (B[j]) ans=max(f[i-1][j-1]+g[A[i]+1][B[j]+1]+C,ans);
    cout<<ans;
    return 0;
}