题目链接: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;
}