2019.1.5【NOIP提高组】模拟B组 JZOJ 3059 雕塑

Question

【问题描述】

Wcyz为了迎接百年校庆，美化校园，请了校友笨笨将n座雕塑，准备安置在校园内，整个校园可以抽象成一个n*n的大网格，每个1*1网格最多只能安置一座雕塑，但是某些1*1的网格上恰好是一个食堂或湖泊，这些网格是不能安置雕塑的，每个雕塑的造型相同，这样同一种安置方案中交换排列都算一种。任意雕塑在同一行或同一列是不合法的方案

学校想知道有多少种安置方案，笨笨想从中选择最好的一种方案，笨笨想请你告诉他方

案种数。

Input

【输入格式】

第一行，两个整数n，m（n<=20,m<=10），用空格隔开，n表示n*n的大网格，m表示不能安置雕塑的位置


第二行至m+1行，每行两个数x，y，用空格分开，表示坐标（x,y）的1*1 的网格上不能安置雕塑。

Output

【输出格式】

一个数，方案种数（方案种数<=2^63-1）

【Sample Input】

6 7
1 1
2 1
2 2
3 3
3 4
4 3
4 4

【Sample Output】

184

【Data Constraint】

Hint

$n<=20,m<=10$

---

思路

方案数=总共的方案数-放在障碍物处的方案数

于是乎我们就可以枚举障碍物处的填补状况了

$However$,这样会多算，因为有些部分是公共控制的，所以，要用荣斥，偶数个时+，奇数个时减即可

---

代码

#include<cstdio>
#include<cstring>
using namespace std;int x[10],y[10],n,m;
long long ans,P[21],k;
bool ok,vx[21],vy[21];
int count(int x)//算出当前状态下选的雕塑的个数
{
	int ans=0;
	for(;x;x-=x&-x) ans++;
	return ans;
}
signed main()
{
	P[0]=1;
	for(register int i=1;i<21;i++) P[i]=P[i-1]*i;//预处理阶乘
	scanf("%d%d",&n,&m);
	for(register int i=0;i<m;i++) scanf("%d%d",x+i,y+i);
	ans=P[n];
	for(long long i=1;i<(1<<m);i++)
	{
		ok=true;k=n;
		memset(vx,0,sizeof(vx));
		memset(vy,0,sizeof(vy));
		for(register int j=0;j<m;j++)
		if((i>>j)&1)
		{
			if(vx[x[j]]||vy[y[j]]) {ok=false;break;}//同一个位置不能重复占用，占用了说明该状态不合法
			k--;
			vx[x[j]]=true;vy[y[j]]=true;
		}
		if(ok) if(count(i)&1) ans-=P[k];else ans+=P[k];//偶数个+，奇数个减
	}
	printf("%lld",ans);
}