关键点是行和列是独立的。
移动K步可以视为x方向移动i步(y向不动),y方向移动K-i步(x向不动)。分别预处理出x和y方向上移动i和K-i步时的种数,然后使用排列组合枚举i。相乘即可。
注意,计算组合数时可能使用了除法这将使模运算出错。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<cmath>
#include<map>
#include<string>
#include<set>
#include<queue>
#include<algorithm>
using namespace std;
typedef long long LL;
const int maxn=1010;
const LL mod=9999991;
LL C[maxn][maxn];
LL dp1[maxn][maxn],dp2[maxn][maxn];
LL dp[maxn][maxn];
const int mv[4]= {1,2,-1,-2};
void init()
{
C[0][0]=1;
for(int i=1; i<=1000; ++i)
{
C[i][0]=C[i-1][0];
for(int j=1; j<=i; ++j)
{
C[i][j]=C[i-1][j-1]+C[i-1][j];
C[i][j]%=mod;
}
}
}
int main()
{
init();
int T,kase=0;
scanf("%d",&T);
while(T--)
{
int n,m,K,x,y;
scanf("%d%d%d%d%d",&n,&m,&K,&x,&y);
memset(dp1,0,sizeof(dp1));
memset(dp2,0,sizeof(dp2));
memset(dp,0,sizeof(dp));
dp1[0][0]=dp1[0][x]=1;
for(int i=1; i<=K; ++i)
{
for(int j=1; j<=n; ++j)
{
for(int k=0; k<4; ++k)
{
int s=j+mv[k];
if(1<=s&&s<=n)
dp1[i][j]+=dp1[i-1][s];
dp1[i][j]%=mod;
}
dp1[i][0]+=dp1[i][j];
dp1[i][0]%=mod;
}
}
dp2[0][0]=dp2[0][y]=1;
for(int i=1; i<=K; ++i)
{
for(int j=1; j<=m; ++j)
{
for(int k=0; k<4; ++k)
{
int s=j+mv[k];
if(1<=s&&s<=m)
dp2[i][j]+=dp2[i-1][s];
dp2[i][j]%=mod;
}
dp2[i][0]+=dp2[i][j];
dp2[i][0]%=mod;
}
}
LL ans=0;
for(int i=0; i<=K; ++i)
{
int j=K-i;
ans+=C[K][i]*dp1[i][0]%mod*dp2[j][0]%mod;
ans%=mod;
}
printf("Case #%d:\n%I64d\n",++kase,ans);
}
return 0;
}