UVA - 10723 Cyborg Genes DP：最长公共子序列扩展

传送门

题意：输入两个串，找一个最短的串使两个串都是它的子序列（可以不连续），求最短的长度和解的个数。

观察可知最短长度就是两串长度之和减去最长公共子序列的长度。一开始还想求出这个长度来之后暴个数。。真是太愚蠢了= =，分析一下，当s1[i]==s2[j]的时候，个数一定和c[i-1][j-1]的个数一样，就是加上当前的字母；如果不相等，那么分三种情况，d[i-1][j]>d[i][j-1]，那么c[i][j] = c[i-1][j]，即在末尾加上当前字母，d[i-1][j]<d[i][j-1]，那么c[i][j] = c[i][j-1]，同在末尾加字母，c[i][j] = c[i-1][j]+c[i][j-1]，两种情况都可以。

然后就蠢了，忘了还有空串的情况，wa了好几次。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int maxn = 55;
char s1[maxn],s2[maxn];
int d[maxn][maxn];
int T,len1,len2,ans,cases = 0;
long long c[maxn][maxn];
int main(){
    //freopen("in.txt","r",stdin);
    scanf("%d%*c",&T);
    while(T--){
        gets(s1+1);
        gets(s2+1);
        len1= strlen(s1+1);
        len2 = strlen(s2+1);
        memset(d,0,sizeof d);
        memset(c,0,sizeof c);
        c[0][0] = 1;
        for(int i=1;i<=len1;i++)
            c[i][0] = 1;
        for(int i=1;i<=len2;i++)
            c[0][i] = 1;
        for(int i=1;i<=len1;i++)
            for(int j=1;j<=len2;j++){
                if(s1[i]!=s2[j]){
                    d[i][j] = max(d[i-1][j],d[i][j-1]);
                    if(d[i-1][j]>d[i][j-1]) c[i][j] = c[i-1][j];
                    else if(d[i-1][j]<d[i][j-1]) c[i][j] = c[i][j-1];
                    else c[i][j] = c[i-1][j]+c[i][j-1];
                }
                if(s1[i]==s2[j]) {
                    d[i][j] = max(d[i][j],d[i-1][j-1]+1);
                    c[i][j] = c[i-1][j-1];
                }
            }
        ans = len1+len2-d[len1][len2];
        printf("Case #%d: %d %lld\n",++cases,ans,c[len1][len2]);
    }
    return 0;
}