简介:
n*m的矩阵中放置k个石子,要求第一行,最后一行,第一列,最后一列必须有石子
分析:
要是没有特殊的限制,在n*m的矩阵中放置k个石子,会求吧:
C(n*m,k)
实际上我们确实可以把这个问题转化一下:
设A为第一行没有石子的方案
B为第一列没有石子的方案
C为最后一行没有石子的方案
D为最后一列没有石子的方案
S为全集
最后答案就是在”S中但是不在ABCD任何一个集合“中的方案数
这样我们就可以用容斥原理求解了
在容斥的时候,我们要减去一些集合的∪,又要加上一些集合的∪,我们可以发现,
只要是奇数个集合的∪一定是减去,偶数个集合的∪一定是加上
那我们怎么枚举集合的所有组合状态呢:
最耿直的:四重循环
但是,这道题我们要换一种高mei级shi精zhao妙shi的方法
我们利用状压dp的思路,用二进制上的每一位代表一个集合
这样只要一层循环即可
tip
状态的枚举要从0开始
C[0][0]要初始化为1
//这里写代码片
#include<cstdio>
#include<cstring>
#include<iostream>
#define ll long long
using namespace std;
const int mod=1000007;
int n,m,k;
int C[503][503];
void cl()
{
int i,j;
memset(C,0,sizeof(C));
C[0][0]=1; C[1][1]=1; C[1][0]=1;
for (i=2;i<=500;i++)
{
C[i][0]=1; C[i][i]=1;
for (j=1;j<i;j++)
C[i][j]=C[i-1][j]+C[i-1][j-1],C[i][j]%=mod;
}
}
int main()
{
int T;
cl();
scanf("%d",&T);
for (int cas=1;cas<=T;cas++)
{
scanf("%d%d%d",&n,&m,&k);
int ans=0;
for (int i=0;i<=(1<<4)-1;i++) //枚举集合的组合状态
{
int sum=0,r=n,c=m;
if (i&1) sum++,r--; //A
if (i&2) sum++,c--; //B
if (i&4) sum++,r--; //C
if (i&8) sum++,c--; //D
if (sum&1) //奇数
ans=(ans+mod-C[r*c][k])%mod;
else //偶数
ans=(ans+C[r*c][k])%mod;
}
printf("Case %d: %d\n",cas,ans);
}
return 0;
}