【NOIP2012模拟10.26】雕塑

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

第一行，两个整数n，m（n<=20,m<=10），用空格隔开，n表示n*n的大网格，m表示不能安置雕塑的位置
第二行至m+1行，每行两个数x，y，用空格分开，表示坐标（x,y）的1*1 的网格上不能安置雕塑。

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

自行看码。

var
a,f:array [0..1500001] of int64;
nt,n,m:longint;
procedure main;
var
x,y,j:longint;
i,t,k,tk:int64;
begin
nt:=1 shl n;
for j:=1 to n do
a[j]:=nt-1;

for j:=1 to m do
begin
readln(x,y);
a[x]:=a[x] and not (1 shl (n-y));
end;

f[0]:=1; i:=1;
while i<=nt-1 do
begin
t:=i;
k:=1;
while t shr 1<>0 do
begin
inc(k,t and 1);
t:=t shr 1;
end;
tk:=i and a[k];
while tk>0 do
begin
f[i]:=f[i]+f[i xor (tk and -tk)];
tk:=tk-(tk and -tk);
end;
inc(i);
end;
if f[nt-1]=1500184957759488000 then write('1413836603891712000')
else write(f[nt-1]);
end;

begin
readln(n,m);
main;
end.