HDU 4681 String 最长公共子序列

题目链接

http://acm.hdu.edu.cn/showproblem.php?pid=4681

题意：

给你a,b,c三个串，构造一个d串使得d是a，b的子序列，并且c是d的连续子串。求d最大的长度。

题解：

枚举a，b串开始匹配c的位置，（结束的位置可以贪心出来），然后前后都用最长公共子序列来跑就可以了。

O(n^2)预处理，O(n^2)枚举。

#pragma comment(linker, "/STACK:102400000,102400000")
#include<cstdio>
#include<cstring>
#include<vector>
#include<algorithm>
#include<queue>
using namespace std;

typedef long long LL;
const int maxn = ;

char sa[maxn], sb[maxn], sc[maxn];
int dp1[maxn][maxn], dp2[maxn][maxn];

int pa[maxn], pb[maxn];
void get_p(char* ss, int* pp) {
    memset(pp, -, sizeof(int)*maxn);

    int lss = strlen(ss),lsc=strlen(sc);
    for (int i = ; i < lss; i++) {
        int k = ,j;
        for (j = i; j < lss; j++) {
            if (k == lsc) break;
            if (sc[k] == ss[j]) {
                k++;
            }
        }
        if (k == lsc){
            pp[i] = j - ;
        }
    }
}

void init() {
    memset(dp1, , sizeof(dp1));
    memset(dp2, , sizeof(dp2));
}

int main() {
    int tc,kase=;
    scanf("%d", &tc);
    while (tc--) {
        init();
        scanf("%s%s%s", sa, sb, sc);
        get_p(sa, pa), get_p(sb, pb);
        int la = strlen(sa), lb = strlen(sb),lc=strlen(sc);
        for (int i = ; i < la; i++) {
            for (int j = ; j < lb; j++) {
                if (sa[i] == sb[j]) {
                    dp1[i][j] = max(dp1[i][j], i>&&j>?dp1[i - ][j - ]+:);
                }
                dp1[i][j] = max(dp1[i][j], i>?dp1[i - ][j]:);
                dp1[i][j] = max(dp1[i][j], j >  ? dp1[i][j - ] : );
            }
        }
        for (int i = la-; i >= ; i--) {
            for (int j = lb - ; j >= ; j--) {
                if (sa[i] == sb[j]) {
                    dp2[i][j] = max(dp2[i][j], dp2[i + ][j + ] + );
                }
                dp2[i][j] = max(dp2[i][j], dp2[i + ][j]);
                dp2[i][j] = max(dp2[i][j], dp2[i][j + ]);
            }
        }
        int ans = ;
        for (int i = ; i < la - ; i++) {
            for (int j = ; j < lb - ; j++) {
                if (pa[i] != - && pb[j] != -) {

                    int lef = i >  && j >  ? dp1[i - ][j - ] : ;
                    int rig = pa[i] < la -  && pb[j] < lb -  ? dp2[pa[i] + ][pb[j] + ] : ;
                    ans = max(ans, lef + lc + rig);
                }
            }
        }
        printf("Case #%d: %d\n",++kase, ans);
    }
    return ;
}