SDUT 动态规划 最长公共子序列

最长公共子序列

Problem Description
从一个给定的串中删去（不一定连续地删去）0个或0个以上的字符，剩下地字符按原来顺序组成的串。例如：
“ ”，“a”，“xb”，“aaa”，“bbb”，“xabb”，“xaaabbb”都是串“xaaabbb”的子序列。（例子中的串不包含引号。）
编程求N个非空串的最长公共子序列的长度。限制：2<=N<=100；N个串中的字符只会是数字0，1，…，9或小写英文字母a，b，…，z；每个串非空且最多含100个字符；N个串的长度的乘积不会超过30000。

Input
文件第1行是一个整数T，表示测试数据的个数（1<=T<=10）。接下来有T组测试数据。各组测试数据的第1行是一个整数Ni，表示第i组数据中串的个数。各组测试数据的第2到N+1行中，每行一个串，串中不会有空格，但行首和行末可能有空格，这些空格当然不算作串的一部分。

Output
输出T行，每行一个数，第i行的数表示第i组测试数据中Ni个非空串的最长公共子序列的长度。
Sample Input
1
3
ab
bc
cd

Sample Output
0

解题思路
　 题目中有一句话：N个串的长度的乘积不会超过30000，这个地方告诉了我们数据的规模并且还强调了N个串长度的乘积，这里就要提一下，在求两字符序列的最长公共字符子序列这个问题中，我们为了方便理解，会画一个m * n的表格（m代表字符串x的长度，n代表字符串y的长度，x，y为任意两字符序列）来观察公共字符个数，两字符序列的最长公共字符子序列问题还可以看成二维公共子序列问题，因为解题的时候我们会建立一个二维数组。
　　该题可以看成一个多维字符求最长公共子序列，但是题目说：N个子符串的乘积不会超过30000。这里我们就可以把这道多维问题转化为一个一维问题，参照求两字符序列的最长公共字符子序列问题把本该建立的多维表格给转化为一维数组。每一个index都对应着唯一的情况，即：多维字符串的长度为x1,x2,x3……xi时是否有公共字符。
　　具体思路如下：
　　建立一个二维字符数组保存字符，建立一个一维数组len来保存每一个字符数组的长度，一维数组L来保存多维字符数组长度的动态变化，用于求index（这部分请结合下面的代码利于理解），我们从每一个字符数组的末尾开始匹配，若都相等则把每一个多维数组长度L减一（所以我称它为动态长度数组），所以最长长度ret = DP（L） + 1；若只要有一个不同，则开始枚举每一种情况，详细请看代码：

#include <bits/stdc++.h>
using namespace std;
char a[110][110];
int dp[30010];
int len[110];
int n;
int DP(int *x)
{
    int i,j,index,ret;
    for(i=1;i<=n;i++)
        if(x[i]==0)
            return 0;
    for(index=x[n]-1,i=n-1;i>=1;i--)
        index=index*len[i]+x[i]-1;
    if(dp[index]>=0)
        return dp[index];
    for(i=2;i<=n;i++)
        if(a[1][x[1]-1]!=a[i][x[i]-1])
        break;
    if(i>n)
    {
        for(j=1;j<=n;j++)
            x[j]--;
        ret=DP(x)+1;
        for(j=1;j<=n;j++)
            x[j]++;
    }
    else
    {
        ret=0;
        int t;
        for(j=1;j<=n;j++)
        {
            x[j]--;
            t=DP(x);
            ret=max(ret,t);
            x[j]++;
        }
    }
    dp[index]=ret;
    return ret;
}
int main()
{
    int t;
    int L[110];
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
        {
            scanf("%s",a[i]);
            len[i]=L[i]=strlen(a[i]);
        }
        memset(dp,-1,sizeof(dp));
        printf("%d\n",DP(L));
    }
    return 0;
}