多校第八场 Hdu 4681 String

看了解题报告才明白，但是解题报告用的是LCS算法，这个算法我不熟悉，所以先暴力一下左子序列和右子序列。。

题意：给出三个序列A,B,C,求一个最长序列D使得：D是A、B的子序列，C是D的子串，输出D的长度

思路：用数组pA[num1][0]，pA[num1][1]分别记录C序列在A序列的起始位置和最近的终点位置

           用数组pB[num1][0]，pB[num1][1]分别记录C序列在B序列的起始位置和最近的终点位置

           Left[i][j]表示A中i位置下B中j位置下最长的公共左子序列

           Right[i][j]表示A中i位置下B中j位置下最长的公共右子序列

对于任意A，B 的起始位置，求出A，B 在起始位置之前的最长公共
子序列，用同样的方法求出与起始位置相应的结束位置以后的最长公
共子序列。两者之和求最大值再加上C的长度即为所求

先进行预处理

代码：

#include <iostream>
#include <stdio.h>
#include <cstring>
#include <cmath>
#include <vector>
#include <algorithm>
#include <map>
#include <queue>
#include <stack>
#define lson l,mid,num<<1
#define rson mid+1,r,num<<1|1
using namespace std;
char a[1005],b[1005],c[1005];
int Left[1005][1005],Right[1005][1005];
int pA[1005][2];
int pB[1005][2];
int main()
{
    int t,la,lb,lc,g=1;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%s",a+1);//要从下标1开始是因为为保证Left[i-1][j-1]合理
        scanf("%s",b+1);
        scanf("%s",c+1);
        la=strlen(a+1);
        lb=strlen(b+1);
        lc=strlen(c+1);
        memset(Left,0,sizeof(Left));
        memset(Right,0,sizeof(Right));
        for(int i=1; i<=la; i++)//预处理
        {
            for(int j=1; j<=lb; j++)
            {
                int cost=0;
                if(a[i]==b[j])
                    cost=1;
                Left[i][j]=Left[i-1][j-1]+cost;
                Left[i][j]=max(Left[i][j],max(Left[i][j-1],Left[i-1][j]));
            }
        }
        for(int i=la; i>=1; i--)//预处理
        {
            for(int j=lb; j>=1; j--)
            {
                int cost=0;
                if(a[i]==b[j])cost=1;
                Right[i][j]=Right[i+1][j+1]+cost;
                Right[i][j]=max(Right[i][j],max(Right[i+1][j],Right[i][j+1]));
            }
        }
        int num1=0,num2=0;
        for(int i=1; i<=la; i++)
        {
            if(a[i]==c[1])
            {
                int j=2,ii=i+1;
                while(j<=lc)
                {
                    while(ii<=la && a[ii]!=c[j])ii++;
                    if(ii>la)break;
                    j++;
                    ii++;
                }
                if(j>lc)
                {
                    pA[num1][0]=i;
                    pA[num1][1]=ii-1;
                    num1++;
                }

            }
        }
        for(int i=1; i<=lb; i++)
        {
            if(b[i]==c[1])
            {
                int j=2,ii=i+1;
                while(j<=lc)
                {
                    while(ii<=lb && b[ii]!=c[j])ii++;
                    if(ii>lb)break;
                    j++;
                    ii++;
                }
                if(j>lc)
                {
                    pB[num2][0]=i;
                    pB[num2][1]=ii-1;
                    num2++;
                }
            }
        }
        int ans=0;
        for(int i=0; i<num1; i++)
        {
            for(int j=0; j<num2; j++)
            {
            ans=max(ans,Left[pA[i][0]-1][pB[j][0]-1]+Right[pA[i][1]+1][pB[j][1]+1]);
            }
        }
        printf("Case #%d: %d\n",g++,ans+lc);
    }
    return 0;
}