hdu 5282 Senior's String

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5282

     以前一直不明白为什么很多组合数学可以用DP来求。BC上很多都是这类数学问题，也有很好的题解。做过很多这类题以后终于明白了。

所以多做题还是有好处的，嘿嘿。

DP为何可以用来求组合问题，DP无非也是计数，组合也是计数。

首先，你要理解这个组合公式： C(n,k) + C(n,k+1) = C(n+1,k+1)

（下面这段话摘自《算法竞赛入门经典训练指南》 p104）

        这是组合数的递推公式，经常用于预处理，证明如下：n+1 个数里面选 k+1 有两类办法。要么选第1个数，要么不选第一个数。

如果不选，则问题转化为 n 个数里选 k + 1 个数；如果选，则问题转化为 n 个数里选 k 个数。这两类办法是不重复不遗漏的，由加法原理得证。

明白这个性质，了解清楚DP的转移过程，做这个题目就简单多了。

先用dp【i】【j】 表示X的前i个字符中最多有多少在Y的前j个中。这个O(n^2) 处理即可。

然后f【i】【j】 表示 X的前i个字符中最多组合有多少在Y的前j个中，分两种情况。

（1） ： dp[i][j] == dp[i-1][j]

                       f[i][j] += f[i-1][j];    // 表示X的第i个字符不选

（2）：在Y字符串的前j个中找到最后面出现X[i]字符的位置 p

           如果 dp[i-1][p] + 1 == dp[i][j]

                   f[i][j] += f[i-1][p];   // 表示可选X的第i个字符

这两种方案是不重复不遗漏的。

#include<bits/stdc++.h>
using namespace std;
const int maxn = 1010;
const int MOD = 1000000007;
char x[maxn],y[maxn];
int dp[maxn][maxn];
int f[maxn][maxn];
int pre[maxn][27];
int tmp[27];
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%s",x);
        scanf("%s",y);
        memset(dp,0,sizeof(dp));
        memset(pre,0,sizeof(pre));
        memset(tmp,0,sizeof(tmp));
        int n = strlen(x),m = strlen(y);

        for(int i=1;i<=m;i++)
        {
            tmp[ y[i-1] - 'a' ] = i;
            for(int j=0;j<26;j++)
                pre[i][j] = tmp[j];
        }
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=m;j++)
            {
                if( x[i-1] == y[j-1] )
                    dp[i][j] = max(dp[i][j],dp[i-1][j-1]+1);
                else dp[i][j] = max( dp[i-1][j] , dp[i][j-1] );
            }
        }
        memset(f,0,sizeof(f));
        for(int i = 0; i <= n; i++ )
        {
            for(int j = 0; j <= m; j++ )
            {
                if( dp[i][j] == 0 )
                {
                    f[i][j] = 1;
                }
                else {
                    if( dp[i-1][j] == dp[i][j] )
                        f[i][j] = ( f[i][j] + f[i-1][j]) % MOD;
                    if( pre[j][ x[i-1]-'a'] )
                    {
                        int p = pre[j][ x[i-1]-'a'];
                        if( dp[i-1][p-1] + 1 == dp[i][j] )
                        f[i][j] = ( f[i][j] + f[i-1][p-1] ) % MOD;
                    }
                }
            }
        }
       printf("%d\n",f[n][m]);
    }
    return 0;
}