hdu 4832 chess dp

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