容斥原理,书上的办法很巧妙,我这里代码比较难看,注意的是在mod之前要保证里面的数是正的
#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <ctime>
#include <cmath>
#include <vector>
#include <deque>
#include <algorithm>
#include <list>
#include <map>
#include <set>
#include <queue>
#include <stack>
#define REP(i,n) for(int i=0;i<(n);i++)
#define REP1(i,n) for(int i=1;i<=(n);i++)
#define FOR(i,a,b) for (int i=(a);i<=(b);i++)
#define CLR(x,n) memset(x,n,sizeof(x))
#define PN printf("\n")
#define read(x) scanf("%d",&x)
#define read2(x,y) scanf("%d%d",&x,&y)
#define read3(x,y,z) scanf("%d%d%d",&x,&y,&z)
#define write(x) printf("%d",x)
#define write1(x) printf("%d ",x)
#define writeln(x) printf("%d\n",x)
#define write2(x,y) printf("%d %d",x,y)
#define writeln2(x,y) printf("%d %d\n",x,y)
#define write3(x,y,z) printf("%d %d %d",x,y,z)
#define writeln3(x,y,z) printf("%d %d %d\n",x,y,z)
#define delta 1e-5
#define INF 1<<30
#define MOD 1000007
using namespace std;
/*
WA 标程中c[0][0]为1
WA 问题出在计算式,出来可能是负的
AC
*/
int c[500][500];
int main()
{
CLR(c,0);
c[0][0]=1;
REP1(n,499){
c[n][0]=1; c[n][n]=1;
REP1(k,n-1)
c[n][k]=( c[n-1][k-1]+c[n-1][k] ) %MOD;
}
//REP(n,10)
//REP(k,10)
//writeln3(n,k,c[n][k]);
int t;read(t);
REP(i,t){
int m,n,k;read3(m,n,k);
int ans=( c[m*n][k]
-2*c[m*n-m][k]-2*c[m*n-n][k]+4*c[m*n-m-n+1][k]+c[m*n-2*m][k]+c[m*n-2*n][k]
-2*c[m*n-n-2*m+2][k]-2*c[m*n-2*n-m+2][k]+c[m*n-2*n-2*m+4][k]+8*MOD )%MOD;
printf("Case %d: %d\n",i+1,ans);
}
//system("pause");
return 0;
}