hdu 4832 dp ***

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<map>
using namespace std;
#define MOD 9999991
const int INF=9999991;
const double eps=1e-5;
typedef long long ll;
#define cl(a) memset(a,0,sizeof(a))
#define ts printf("*****\n");
const int MAXN=2050;
int n,m,tt;
int x,y,k;
int c[1050][1050];
int sum1[1050],sum2[1050];
int dp1[MAXN][1050],dp2[MAXN][1050];
void add(int &a,int b)
{
a=(a+b)%MOD;
}
void fun()
{
cl(dp1),cl(dp2);
c[0][0]=c[1][0]=c[1][1]=1;
for(int i=2;i<=1010;i++)
{
c[i][0]=1;
for(int j=1;j<i;j++)
{
c[i][j]=c[i-1][j]+c[i-1][j-1];
c[i][j]%=MOD;
}
c[i][i]=1;
}
dp1[x][0]=1;
for(int i=1;i<=k;i++)
{
for(int j=1;j<=n;j++)
{
if(j+1<=n)  add(dp1[j][i],dp1[j+1][i-1]);
if(j+2<=n)  add(dp1[j][i],dp1[j+2][i-1]);
if(j-1>=1)  add(dp1[j][i],dp1[j-1][i-1]);
if(j-2>=1)  add(dp1[j][i],dp1[j-2][i-1]);
}
}
dp2[y][0]=1;
for(int i=1;i<=k;i++)
{
for(int j=1;j<=m;j++)
{
if(j+1<=m)  add(dp2[j][i],dp2[j+1][i-1]);
if(j+2<=m)  add(dp2[j][i],dp2[j+2][i-1]);
if(j-1>=1)  add(dp2[j][i],dp2[j-1][i-1]);
if(j-2>=1)  add(dp2[j][i],dp2[j-2][i-1]);
}
}
cl(sum1),cl(sum2);
for(int i=0;i<=k;i++)
{
for(int j=1;j<=n;j++)
{
add(sum1[i],dp1[j][i]);
}
}
for(int i=0;i<=k;i++)
{
for(int j=1;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=1;
scanf("%d",&tt);
while(tt--)
{
scanf("%d%d%d%d%d",&n,&m,&k,&x,&y);
fun();
printf("Case #%d:\n",ca++);
ll sum=0;
for(i=0;i<=k;i++)
{
sum += (long long)c[k][i]*sum1[i]%MOD*sum2[k-i]%MOD;
sum %= MOD;
}
printf("%d\n",(int)sum);
}
}