dp1[i][j]表示只走x轴走j步到i位置有多少总走法,dp2同,dp方程就很好写
wa了无数发,发现MOD写在INF上了
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<map>
using namespace std;
#define MOD 9999991
const int INF=;
const double eps=1e-;
typedef long long ll;
#define cl(a) memset(a,0,sizeof(a))
#define ts printf("*****\n");
const int MAXN=;
int n,m,tt;
int x,y,k;
int c[][];
int sum1[],sum2[];
int dp1[MAXN][],dp2[MAXN][];
void add(int &a,int b)
{
a=(a+b)%MOD;
}
void fun()
{
cl(dp1),cl(dp2);
c[][]=c[][]=c[][]=;
for(int i=;i<=;i++)
{
c[i][]=;
for(int j=;j<i;j++)
{
c[i][j]=c[i-][j]+c[i-][j-];
c[i][j]%=MOD;
}
c[i][i]=;
}
dp1[x][]=;
for(int i=;i<=k;i++)
{
for(int j=;j<=n;j++)
{
if(j+<=n) add(dp1[j][i],dp1[j+][i-]);
if(j+<=n) add(dp1[j][i],dp1[j+][i-]);
if(j->=) add(dp1[j][i],dp1[j-][i-]);
if(j->=) add(dp1[j][i],dp1[j-][i-]);
}
}
dp2[y][]=;
for(int i=;i<=k;i++)
{
for(int j=;j<=m;j++)
{
if(j+<=m) add(dp2[j][i],dp2[j+][i-]);
if(j+<=m) add(dp2[j][i],dp2[j+][i-]);
if(j->=) add(dp2[j][i],dp2[j-][i-]);
if(j->=) add(dp2[j][i],dp2[j-][i-]);
}
}
cl(sum1),cl(sum2);
for(int i=;i<=k;i++)
{
for(int j=;j<=n;j++)
{
add(sum1[i],dp1[j][i]);
}
}
for(int i=;i<=k;i++)
{
for(int j=;j<=m;j++)
{
add(sum2[i],dp2[j][i]);
}
}
}
int main()
{
int i,j;
#ifndef ONLINE_JUDGE
freopen("1.in","r",stdin);
#endif
int ca=;
scanf("%d",&tt);
while(tt--)
{
scanf("%d%d%d%d%d",&n,&m,&k,&x,&y);
fun();
printf("Case #%d:\n",ca++);
ll sum=;
for(i=;i<=k;i++)
{
sum += (long long)c[k][i]*sum1[i]%MOD*sum2[k-i]%MOD;
sum %= MOD;
}
printf("%d\n",(int)sum);
}
}