【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）

题解

 自行看码。

代码

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.