jzoj 3059. 【NOIP2012模拟10.26】雕塑

Description

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

解析

由于急着打第一题，这题只简单地打了个暴力就溜了。

正解

看到数据范围n<=20,m<=10,方案种数<=2^63-1。

由于方案总数较大，若采用单纯的搜索算法，很难做到不超时。

其实观察一下数据规模中的n和m，m<=10,这是一条重要的信息，m比n*n少得多，那么雕塑“安置在禁区内的方案数”就必然比“安置在禁区外的”方案数少得多，

"安置在禁区内的雕塑仍可以用dfs，但是必须搞清一个问题：安置在禁区内的方案总数与求不在禁区的方案数又有什么关系呢，这就恰好就是容斥原理中淘汰原则的十分巧妙的运用。

再来进行进一步的分析：用rk表示把k个雕塑放在了n*n的方格上，并且这k个雕塑都处在禁区放置位置上的方案数，根据容斥原理的推论n个雕塑都安置在非禁区内的方案数等于

n!-r1*(n-1)!+r2*(n-2)-r3*(n-3)!...(-1)^k*rk*(n-k)!...+(-1)^n*rn

（摘自jzoj题解）

 

CODE（本蒟蒻用的是容斥）

var
        h,l,pow,r:array[0..1000]of int64;
        a:array[0..100,0..5]of int64;
        i,j,n,m:longint;
        ans:int64;
procedure dfs(x,y,p:longint);
var
        i:longint;
begin
        if y=0 then begin
                inc(r[p]);
                exit;
        end
        else begin
                for i:=x+1 to m do begin
                        if (h[a[i,1]]=0)and(l[a[i,2]]=0) then begin
		                h[a[i,1]]:=1;
                                l[a[i,2]]:=1;
		                dfs(i,y-1,p);
                                h[a[i,1]]:=0;
                                l[a[i,2]]:=0;
                        end;
                end;
        end;
end;
begin
        readln(n,m);
        for i:=1 to m do read(a[i,1],a[i,2]);
        for i:=1 to m do dfs(0,i,i);
	ans:=1;
	pow[0]:=1;pow[1]:=1;
        for i:=2 to n do begin
                pow[i]:=pow[i-1]*i;
                ans:=ans*i;
        end;
	for i:=1 to m do
	        if odd(i) then
                        ans:=ans-r[i]*pow[n-i]
	        else
                        ans:=ans+r[i]*pow[n-i];
        writeln(ans);
end.