题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4832
思路:感觉这道题是道好题,考的知识点挺多,都不难,但要组合起来,还是很有难度的
总结下知识点
(1)问题分解, 把一个复杂的问题分解成若干个独立的简单的问题, 并加以求解, 这是紫书P237的原话。以后对于涉及到行和列的问题, 就要联想到行和列是否无关。
对于这道题,从a点开始走k步,并没有要求最后走到的点,所以行和列 可以随便走。
(2)设dp[i][j]为第i步走到第j行时的方案数, 那么dp[i]j = sum( dp[i-1[j+d] ) d = 1, -1, 2, -2;
(3) 有k步,假设有a步走了行,b步走了列,问总共有多少种走法,就是c(k,a)种走法。 所以要先处理下杨辉三角
加上看代码后,应该就可以明白了
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;
const int mod = 9999991;
const int maxn = 1010;
int dir[] = {1, -1, 2, -2};//行列移动方向
int n, m, k, x, y;
int c[maxn][maxn];
int dp1[maxn][maxn], dp2[maxn][maxn];
int sum1[maxn], sum2[maxn];
void init() //计算组合数,杨辉三角
{
c[0][0] = 1;
for(int i=1; i<1010; i++)
{
c[i][0] = c[i][i] = 1;
for(int j=1; j<i; j++)
{
c[i][j] = (c[i-1][j] + c[i-1][j-1]) % mod;
}
}
}
bool judge(int x, bool flag)
{
if(flag)//行
{
if(x<1 || x>n)
return false;
return true;
}
else//列
{
if(x<1 || x>m)
return false;
return true;
}
}
void row()
{
dp1[0][x] = 1;//第0步在第x行的方案数是1
for(int i=1; i<=k; i++)
{
for(int j=1; j<=n; j++)
{
for(int d=0; d<4; d++)
{
int xx = j + dir[d];//上一步所在的行
if(judge(xx, true))
{
dp1[i][j] = (dp1[i][j] + dp1[i-1][xx]) % mod;
}
}
}
}
for(int i=0; i<=k; i++)
{
for(int j=1; j<=n; j++)
sum1[i] = (sum1[i] + dp1[i][j]) % mod;//sum1[i]表示第i步可以走到的行的方案数
}
}
void column()
{
dp2[0][y] = 1;
for(int i=1; i<=k; i++)
{
for(int j=1; j<=m; j++)
{
for(int d=0; d<4; d++)
{
int yy = j + dir[d];
if(judge(yy, false))
{
dp2[i][j] = (dp2[i][j] + dp2[i-1][yy]) % mod;
}
}
}
}
for(int i=0; i<=k; i++)
{
for(int j=1; j<=m; j++)
sum2[i] = (sum2[i] + dp2[i][j]) % mod;
}
}
int main()
{
init();
int t, cas = 1;
scanf("%d", &t);
while(t--)
{
scanf("%d%d%d%d%d", &n, &m, &k, &x, &y);
memset(dp1, 0, sizeof(dp1));
memset(dp2, 0, sizeof(dp2));
memset(sum1, 0, sizeof(sum1));
memset(sum2, 0, sizeof(sum2));
row();
column();
ll ans = 0;
for(int i=0; i<=k; i++)
ans = (ans + (ll)c[k][i]*sum1[i]%mod*sum2[k-i]%mod) % mod;//注意这里可能会溢出
printf("Case #%d:\n", cas++);
printf("%I64d\n", ans);
}
return 0;
}