啦啦队 (容斥原理)

【题目描述】在一个m行n列的矩形格里放k个相同的石子，文友多少种方法？每个格子最多放体格石子，所有石子都要用完，并且第一行，最后一行，第一列，最后一列都得有石子。
【输入格式】
输入第一行为数据数组T（T<=50），每组数据包含整数m，n，k(2<= n,m <=20,k<=500)。
【输出格式】
对于每组数据，输出方案总数除以1000 007的余数。

#include<catdio>
#include<cstring>
using namespace std;

const int MOD = 1000007;
const int MAXN = 500;
int main()
{
	memset(c,0,sizeof(c));
	for(int i=0;i<MAXN;i++)
	{
		c[i][0]=1;c[i][i]=1;//有定义可知每个数的边界值
		for(int j=0; j<i ;j++)
		{
			c[i][j]=(c[i-1][j-1]+c[i-1][j])%MOD;//由递推公式来用数组c来存储c(m,n)的值，行坐标是数据下底而上底是要取得个数
		}
	}
	int T;
	scanf("%d",&T);
	while(T--)
	{
		int n,m,k;//输入行数，列数，要去的石子的个数
		scanf("%d",&n,&m,&k);
		for(int i=0;i<16;i++)//用从 1~16 的二进制表示了4中情况的组合情况，可以把这当做一个模板。
		{
			int c=n;
			int l=m;//l和c分别表示行数列数可选择的个数 
			int b=0;//表示集合数
			int ans=0;
			if(i&1) c--,b++;//第一列没有石头，可以石头的列数有r减一行。
			if(i&2) c--,b++
			if(i&4) l--,b++;
			if(i&8) l--,b+++;
			if(b&1)
				ans=(ans+MOD-c[l*c][k])%MOD;//保证差不为负数
			else
				ans=(ans+c[l*c][k])%MOD;						
		}
		printf("Case %d:\n",ans);
	}
	return 0;
}

分析：设满足“第一行没有石子“的方案集为A,最后一行没有石子的方案为B，第一列没有石子方案集为C，最后一列没有石子的方案集为D，全集为S，则所求答案就是“在S中但不在集合A，B，C，D任何一个集合中·“的元素个数，可以理解为容斥原理求解。