codevs1050

题目地址：http://codevs.cn/problem/1050/

分析：

最开始想直接用状压做，发现怎么都想不出来，就和当年的多行多米诺骨牌（这道题至少最后还是把普通状压做法看懂了）。

直到听到 @tsyao 神牛说这个是轮廓线状压。

大白书（LRJ算法竞赛）第六章第一节好像就是吧。

 

主体思路：

一。做法思路比较简单，写起来蛋疼死我了= =，用四进制来保存状态（连通性问题，第六章第三个例题算是这个题的提高）。直接位运算比较方便，但是容易出错，注意运算符的优先级。

二。把整个棋盘分成5*N个阶段，对于每个阶段枚举上一个阶段的所有情况，做轮廓线DP。注意当前点在行首的特判，因为这个时候它的左边是没有点的。还有第一行第一个点的情况，我的做法是开始虚构一行，将其除00000外所有状态的值都赋值为inf，00000 赋值为0。

三。做完DP后，答案就是最后一个阶段里值最小的一个。

 

一些小细节：

一。用四进制来判断连通性的正确性证明：

首先用四进制可以用位运算，这个是选择它的好处之一。

而这道题比较简单，一个阶段只有5个格子，对于该阶段每一个没有联通的黑色（假设白色为0的话，黑色为正整数），其值不同，已经连接好的值是一样的，我们令最初出现的黑色为1，与它联通的都为1，其他的为2，3，4，5，依此类推。

(1)但是由于一个阶段只有5个格子，而且如果五个格子里面没有1却有其他数字的黑色，说明不合法，所以黑色最大为5。

(2)一个阶段中5个格子最多只有3种黑色：如果5个格子均在一行，最多只有  “黑 白 黑 白 黑” 这种情况，如果不在一排，也只有   “line1： 黑 白 黑 白 line2：黑”  这种情况。所以 0 1 2 3就可以把所有情况表示出来，因此使用四进制可以解决。

二。在每次到达一个新阶段并且枚举上一个阶段后：

(1)将旧阶段左移2位，判断：1-> 其上方是否为黑色，是的话直接赋值为上方黑色的值。 2-> 其上方如果不是黑色，左边如果有格子且为黑色，则赋值为该黑色的值。 3-> 否则赋值为3。

(2)处理该阶段，将相邻的黑色的值赋值为他们中值较小的黑色的值。

(3)上面的处理完之后检查一下该阶段，如果其中只有1和3的黑色，将所有3赋值为2。

(4)处理完一个阶段后，要将四进制中第六位移出(即二进制中第11 12位)，并且进行判断，如果移出的不为1，且在现有状态中没有其他数值和它相同的格子，那么新状态不合法，要赋值inf。

(5)第一个1之前的所有格子不处理，最后一个1之后的格子都不处理(正确性很好证明，想想就明白了)。

 

注：

(1)主要就是这些问题了，如果还有小问题，可以自己模拟检查错误，手写一两组数据就行了。这里感谢wr_1737的手写数据帮助我AC= =。

(2)我DEBUG的方法很原始，PRINTF大法= =。打出来的16W行数据我看了两三遍，DEBUG接近两个小时才过= =，细节错误一个一个找出来的。

(3)下面发一个代码，里面的DEBUG过程我注释掉了，可以删了再看，还有一些DEBUG过程中加的思路注释，可以看看。代码结合上面的分析看起来会容易点。

代码：http://paste.ubuntu.com/7299284/

/*我是一个从来不发BLOG骗访问量的真诚的人*/

/*

为了那16W行的DEBUG数据，觉得有帮助的同学点一个赞吧！

(

求验证的同学可以把DEBUG注释全打开，然后输入以下输入：

5
10101
01010
10101
01010
10101

)

代码：

const
    maxn=102;
type
    node=record
      x,y:longint;
    end;
    aa=array[0..6]of longint;
var
    f:array[0..maxn,0..1204]of longint;
    a:array[0..maxn]of longint;
    fa:array[1..4]of longint;
    flag:array[1..4]of boolean;
    n:longint;

procedure init;
var
    i,j:longint;
    s:char;
begin
    readln(n);
    for i:=1 to n do
      begin
        for j:=1 to 5 do
          begin
            read(s);
            a[i]:=a[i]<<1+ord(s)-ord('0');
          end;
        readln;
      end;
    while n>0 do
      begin
        if a[n]>0 then break;
        dec(n);
      end;
    if n=0 then
    begin
      write(0);
      halt;
    end;
end;

procedure change(var a:aa;b,c:longint);
var
    i:longint;
begin
    for i:=1 to 5 do
      if a[i]=b then
      begin
        a[i]:=c;
        if (a[i-1]<>0) and (a[i-1]<10) then change(a,a[i-1],c);
        if (a[i+1]<>0) and (a[i+1]<10) then change(a,a[i+1],c);
      end;
end;

procedure get(var a:aa);
var
    i,c:longint;
begin
    c:=10;
    for i:=1 to 5 do
      if (a[i]<>0) and (a[i]<10) then
      begin
        inc(c);
        change(a,a[i],c);
      end;
    for i:=1 to 5 do
      if a[i]>0 then dec(a[i],10);
end;

function bit(x:longint):longint;
begin
    if x=0 then exit(0);
    exit(bit(x-(x and -x))+1);
end;

var
    q:array[0..maxn*1024]of node;

procedure work;
var
    head,tail,i,j,k,ans,save:longint;
    s,t:aa;
    flag:boolean;
begin
    fillchar(f,sizeof(f),1);
    t[0]:=0;
    t[6]:=0;
    ans:=500;
    f[0,0]:=0;
    q[1].x:=0;
    q[1].y:=0;
    head:=1;
    tail:=1;
    while head<=tail do
      begin
        save:=q[head].y;
        for i:=1 to 5 do
          begin
            s[i]:=q[head].y and 3;
            q[head].y:=q[head].y>>2;
          end;
        q[head].y:=save;
        if q[head].x=n then
        begin
          flag:=true;
          for i:=1 to 5 do
            if s[i]>1 then flag:=false;
          if flag then
          if ans>f[q[head].x,q[head].y] then ans:=f[q[head].x,q[head].y];
          inc(head);
          continue;
        end;
        for i:=0 to 31 do
          if i and a[q[head].x+1]=0 then
          begin
            for j:=1 to 5 do
              t[j]:=(((a[q[head].x+1]+i)>>(j-1))and 1)*(j+3);
            k:=0;
            for j:=1 to 5 do
              if (s[j]>0) and (t[j]>0) then
              begin
                t[j]:=s[j];
                k:=k or (1<<s[j]);
              end;
            flag:=true;
            for j:=1 to 5 do
              if (s[j]>0) and (k and (1<<s[j])=0) then flag:=false;
            if flag=false then continue;
            get(t);
            k:=0;
            for j:=5 downto 1 do
              k:=k<<2+t[j];
            if f[q[head].x+1,k]>500 then
            begin
              inc(tail);
              q[tail].x:=q[head].x+1;
              q[tail].y:=k;
            end;
            if f[q[head].x+1,k]>f[q[head].x,q[head].y]+bit(i) then f[q[head].x+1,k]:=f[q[head].x,q[head].y]+bit(i);
          end;
        inc(head);
      end;
    write(ans);
end;

begin
    init;
    work;
end.