CodeForces 128 C.Games with Rectangle（组合数学）

Description

一个$n\times m$的矩阵，每次可以取出一个子矩阵代替原来的矩阵（其边界不能与原矩阵重合），问有多少种进行$k$次操作的方案，两种方案视为不用当且仅当在某一步取的矩阵不一样

Input

三个整数$n,m,k(1\le n,m,k\le 1000)$

Output

输出进行$k$次操作的方案数，结果模$10^9+7$

Sample Input

3 3 1

Sample Output

1

Solution

问题转化为在这个矩阵里画$k$个越来越小的矩阵且保证矩阵不相交，其实就是从非边界的$n-1$的行分界线选出$2k$个作为这$k$个矩阵的上下边界，从$m-1$个非边界的列分界线选出$2k$个作为这$k$个矩阵的左右边界，故方案数为$C_{n-1}^{2k}\cdot C_{m-1}^{2k}$

Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
const int INF=0x3f3f3f3f,maxn=1005;
#define mod 1000000007
int fact[maxn],inv[maxn];
void init(int n=1000)
{
    fact[0]=1;
    for(int i=1;i<=n;i++)fact[i]=(ll)i*fact[i-1]%mod;
    inv[1]=1;
    for(int i=2;i<=n;i++)inv[i]=mod-(ll)(mod/i)*inv[mod%i]%mod;
    inv[0]=1;
    for(int i=1;i<=n;i++)inv[i]=(ll)inv[i-1]*inv[i]%mod;
}
int C(int n,int m)
{
    if(m>n)return 0;
    return (ll)fact[n]*inv[m]%mod*inv[n-m]%mod;
}
int main()
{
    init();
    int n,m,k;
    while(~scanf("%d%d%d",&n,&m,&k))printf("%d\n",(ll)C(n-1,2*k)*C(m-1,2*k)%mod);
    return 0;
}