10723 - Cyborg Genes

好久没做lrj上的动归了，现在感觉每做出一道都是那么的令人兴奋。

这道题搁了好长时间了，第一次记得是暑假开始不久看的，没做出来，现在终于AC了，真是不容易……

动规来解，但应该怎么动？

首先，先求出两个串的最长公共子序列d，然后用两个串的总长度m+n减去d得到的值就是满足题意的最短长度！（这个结论比较直观）

然后，定义状态f(i, j, k)表示第一个串的前i个字符和第二个串的前j个字符满足题意的最短长度为k的串的个数。

状态转移方程：

f(i, j, k) = f(i - 1, j - 1, k - 1) (s[i] == s[j]) 或者 =f(i, - 1 j, k - 1) + f(i, j - 1, k - 1)  (s[i] != s[j])

初始化有点麻烦：f(0, j, j) = f(i, 0, i) = 1;' if(i > k || j > k) f(i, j, k) = 0等

注意题目中答案不能用int存，用unsigned或者long long

还要注意可能有空串，所以要用gets();

程序代码：

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
#define LL long long

LL f[40][40][80];
int d[40][40];
char s[100], t[100];
int n, m, len;
LL dp(int i, int j, int k)
{
    LL & ans = f[i][j][k];
    if(ans != -1) return ans;
    else if(i > k || j > k) return ans = 0;
    else if((i == 0 && j == k) || (j == 0 && i == k)) return ans = 1;
    else if((i == 0 && j != k) || (j == 0 && i != k)) return ans = 0;
    if(s[i - 1] == t[j - 1]) ans = dp(i - 1, j - 1, k - 1);
    else ans = dp(i - 1, j, k - 1) + dp(i, j - 1, k - 1);
    return ans;
}
int main()
{
    //freopen("input.txt", "r", stdin);
    int nCase, con = 1;
    scanf("%d", &nCase);
    getchar();
    while(nCase--){
        gets(s);
        gets(t);
        memset(d, 0, sizeof(d));
        n = strlen(s); m = strlen(t);
        for(int i = 1; i <= n; i++){
            for(int j = 1; j <= m; j++){
                if(s[i - 1] == t[j - 1])
                    d[i][j] = d[i - 1][j - 1] + 1;
                else if(d[i - 1][j] > d[i][j - 1])
                    d[i][j] = d[i - 1][j];
                else d[i][j] = d[i][j - 1];
            }
        }
        len = n + m - d[n][m];
        memset(f, -1, sizeof(f));
        printf("Case #%d: %d %lld\n", con++, len, dp(n, m, len));
    }
    return 0;
}