HDU 4681 string 求最长公共子序列的简单DP+暴力枚举

先预处理，用求最长公共子序列的DP顺着处理一遍，再逆着处理一遍。

再预处理串a和b中包含串c的子序列，当然，为了使这子序列尽可能短，会以c 串的第一个字符开始 ，c 串的最后一个字符结束

将这些起始位置先记录下来，然后枚举这些位置，最大的值输出，看一下代码，你就会顿悟了····哈哈。

贴代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1005
using namespace std;
char a[2][N],b[N];//a[0]为串a,a[1]为串b,b为串c
int dp1[N][N],dp2[N][N];
int len[3],num[2];//len[2]为串c的长度
struct node
{
    int s,e;
} p[2][N];
void init()//dp顺序，逆序求一次最长公共子序列
{
    for(int i=0; i<=len[1]+1; ++i)
    {
        dp1[0][i] = 0;
        dp2[len[0]+1][i]=0;
    }
    for(int i=0; i<=len[0]+1; ++i)
    {
        dp1[i][0] = 0;
        dp2[i][len[1]+1] =0;
    }
    for(int i=1; i<=len[0]; ++i)
    {
        for(int j=1; j<=len[1]; ++j)
        {
            if(a[0][i] == a[1][j] ) dp1[i][j] = dp1[i-1][j-1]+1;
            else dp1[i][j] = max(dp1[i-1][j],dp1[i][j-1]);
        }
    }
    for(int i=len[0]; i>=1; --i)
    {
        for(int j=len[1]; j>=1; --j)
        {
            if(a[0][i]==a[1][j]) dp2[i][j]=dp2[i+1][j+1]+1;
            else dp2[i][j] = max(dp2[i+1][j],dp2[i][j+1]);
        }
    }
}
//从串中找一个包含c串的以c[0]为开始，以c[len-1]为结束的一段，记录下起始位置
int pipei(int ind,int start)
{
    int i,j=1;
    for(i=start+1; i<=len[ind]; ++i)
    {
        if(a[ind][i] == b[j]) ++j;
        if(j==len[2]) break;
    }
    if(j <len[2]) return -1;
    else return i;
}
//预处理出a串和b串中所有如上所述的子串的起始位置
void yumeiju(int ind)
{
    num[ind] = 0;
    for(int i=1; i<=len[ind]; ++i)
    {
        if(a[ind][i] == b[0])
        {
            int e = pipei(ind,i);
            if(e == -1) continue;
            p[ind][num[ind]].s=i;
            p[ind][num[ind]].e = e;
            ++num[ind];
        }
    }
}
//枚举所有的可能，求解
int meiju()
{
    init();
    yumeiju(0);
    yumeiju(1);
    int ans=0;
    for(int i=0; i<num[0]; ++i)
    {
        for(int j=0; j<num[1]; ++j)
        {
            if(dp1[p[0][i].s-1][p[1][j].s-1] + dp2[p[0][i].e+1][p[1][j].e+1] > ans)
                ans = dp1[p[0][i].s-1][p[1][j].s-1] + dp2[p[0][i].e+1][p[1][j].e+1] ;
        }
    }
    return ans +len[2];
}
int main()
{
//    freopen("in.txt","r",stdin);
    int t;
    scanf("%d",&t);
    a[0][0]='#';
    a[1][0] ='#';
    for(int ser =1; ser<=t; ++ser)
    {
        scanf("%s%s%s",a[0]+1,a[1]+1,b);
        len[0]= strlen(a[0])-1;
        len[1] = strlen(a[1])-1;
        len[2] = strlen(b);
        printf("Case #%d: %d
",ser,meiju());
    }
    return 0;
}