bzoj1087: [SCOI2005]互不侵犯King (codevs2451) 状压dp

唔...今天学了状压就练练手...

点我看题

这题的话，我感觉算是入门题了QAQ...

然而我还是想了好久...

大致自己推出了方程，但是一直挂，调了很久选择了题解 坚持不懈的努力的调代码。

然后发现题解的方程和我一毛一样呀QAQ

然后就又开始无限调调调QAQ最后发现自己统计答案没放循环里...气哇QAQ看都之后哇的一声哭出来。

好吧进入正题吧。

设 f[i,j,k] 表示 放到第 i 行，第 i 行状态为 j 且总共放了 k 个的方案数。

然后方程就是  f[i,j,k+num[j]]+=f[i-1,x,k]

就是枚举一个 x 表示 第 i-1 层的状态，然后判其是否合法。

刚开始的时候懒的预处理，后来发现预处理代码更短就改了QAQ是的我是懒兔纸。

所以我萌预处理一个 num[i] 表示 状态为 i 中的 1 个数有多少个,也就是放了多少个国王。

顺便对于 不合法的状态 i 的num[i]=-1

然后对于枚举的每一个 x 直接暴力判 j 和 x 这两个状态合起来后是否合法。

暴力判的原因是我懒得预处理QAQ

然后 就枚举 k 表示 前 i-1 行已经放了 k 个国王。

辣么k+num[j] 就是第 i 行放完后的个数。

显然 k 的范围是 num[x]<=k<=m-num[j]

然后就愉快的水过去了，虽然wa了一次，因为longlong。

总的来说今天的状态压缩学得还不错=v=

 var
     i,j,k,x:longint;
     n,m:longint;
     f:array[..,..,..]of int64;
     ans:int64;
     num:array[..]of longint;
 function isnot(i,x:longint):boolean;
 begin
   exit(( << (i-))and x>);
 end;
 function ok(x:longint):longint;
 var i:longint;
     num:longint;
 begin
   num:=;
   for i:= to n do
   begin
     if isnot(i,x) and((i>)and isnot(i-,x)) then exit(-);
     if isnot(i,x) then inc(num);
   end;
   exit(num);
 end;
 function check(x,y:longint):boolean;
 var i:longint;
 begin
   for i:= to n do
   begin
     if isnot(i,x) then
       if ((i>)and(isnot(i-,y)))or(isnot(i,y))or((i<n)and(isnot(i+,y))) then exit(false);
   end;
   exit(true);
 end;
 begin
   read(n,m);
   f[,,]:=;
   for j:= to ( << n)- do
   num[j]:=ok(j);
   for i:= to n do
   for j:= to ( << n)- do
   begin
     if num[j]>= then
     begin
       for x:= to( << n)- do
       if (num[x]>=)and(check(j,x)) then
       begin
         for k:=num[x] to m do
         if k+num[j]<=m then
           inc(f[i,j,k+num[j]],f[i-,x,k]);
       end;
     end;
     if i=n then inc(ans,f[i,j,m]);
   end;
   writeln(ans);
 end.

bzoj1087