POJ1739 基于连通性的状态压缩dp

题目：http://poj.org/problem?id=1739

这种题新手必须仔细钻研大牛们的代码！！！！这一点很重要

大体上与陈丹琪的例题差不多，就是n*m的矩阵求从左下角到右下角的一条满足遍历除障碍外的所有格点的路的方案数。

首先膜拜陈丹琪的论文：http://wenku.baidu.com/view/f49801155f0e7cd18425368e.html

定义dp[l,i,j]为l行以上(包括l行[1,1]到[l,i])中已构成的所有连通块构成连通的方案数(j为当前轮廓线压缩后的状态，用括号表示法)。

1为(,2为),0为无括号，三进制。

十进制数num在三进制中第k位是 num div power(3,k-1) mod 3

类似于二进制取第k位 num and power(2,k-1)>>k-1 = num div power(2,k-1) mod 2 = num>>k-1 mod 2

dp[0,0,0]=1

(注意我的l,i是从0开始到n-1或m-1)

转移就很简单了：引用梁盾大牛的讲述：

1.这个方格上面没有插头，那么也就是意味着路径在这里打了个弯，因此生成一左一右两个插头：

2.这个方格上方或左方有一个插头，那么他可以继续延续，插在下方或者是右方：

3.这个方格上有一右一左两个插头，直接删除合并：

 

4.这个方格上游一左一右两个插头，直接pass，因为用他转移的话路径将会出现不连通的环;

5.这个方格上有两个左插头，同样意味着路径在这里打了个弯，我们需要把这两个左插头删掉，并且把较靠右的左插头所匹配的右插头改成左插头（开始想成把右边第一个右插头删掉，无限WA。。）：

6.这个方格上有两个右插头，同上：

 

 

http://blog.csdn.net/jasonzhu8/archive/2010/08/01/5779518.aspx

有一个特殊情况就是如果i=m 注意不是m-1(m-1时不用特殊处理)

且j div power(3,m)=0 时

就要换行，当前为dp[l,i,j],则转移为：dp[l+1][0][j*3]:=dp[l+1][0][j*3]+dp[l][i][j]

 

详见代码(313MS,神奇的比盾爷的程序快(800MS))

再次注意虽然我说十进制数num在三进制中第k位是 num div power(3,k-1) mod 3 但由于我的行列是从0开始

所以我取位写成了num div power(3,k) mod 3

program poj1739; const th:array[0..9]of longint=(1,3,9,27,81,243,729,2187,6561,19683); var ch:char;bool:array[0..9,0..9]of boolean; dp:array[0..8,0..10,0..21000]of longint; n,m,i,j,l,k,plugleft,plugup,trans,ls,sum:longint; function checkbyte_three(source,place:longint):longint; begin exit(source div th[place] mod 3); end; begin readln(n,m); repeat fillchar(bool,sizeof(bool),false); fillchar(dp,sizeof(dp),0); for i:=0 to n-1 do begin for j:=0 to m-1 do begin read(ch); if ch='#' then bool[i,j]:=true; end; readln; end; dp[0,0,0]:=1; for l:=0 to n do for i:=0 to m+1 do begin if (l=n-1)and(i=m) then break else for j:=0 to th[m+1] do if dp[l,i,j]<>0 then begin if i=m then begin if j div th[m]=0 then dp[l+1][0][j*3]:=dp[l,i,j]+dp[l+1][0][j*3]; continue; end; plugleft:=checkbyte_three(j,i); plugup:=checkbyte_three(j,i+1); if bool[l,i] then if (plugleft=0)and(plugup=0) then begin dp[l,i+1,j]:=dp[l,i,j]+dp[l,i+1,j]; continue; end else continue; if (plugleft=1)and(plugup=2) then continue; if (plugleft=0)and(plugup=0) then begin trans:=j+th[i]+th[i+1]*2; dp[l][i+1][trans]:=dp[l][i+1][trans]+dp[l][i][j]; continue; end; if (plugleft=1)and(plugup=0)or(plugleft=0)and(plugup=1) then begin ls:=j-plugleft*th[i]-plugup*th[i+1]; trans:=ls+th[i]; dp[l][i+1][trans]:=dp[l][i+1][trans]+dp[l][i][j]; trans:=ls+th[i+1]; dp[l][i+1][trans]:=dp[l][i+1][trans]+dp[l][i][j]; continue; end; if (plugleft=2)and(plugup=0)or(plugleft=0)and(plugup=2) then begin ls:=j-plugleft*th[i]-plugup*th[i+1]; trans:=ls+2*th[i]; dp[l][i+1][trans]:=dp[l][i+1][trans]+dp[l][i][j]; trans:=ls+2*th[i+1]; dp[l][i+1][trans]:=dp[l][i+1][trans]+dp[l][i][j]; continue; end; if (plugleft=2)and(plugup=1) then begin trans:=j-2*th[i]-th[i+1]; dp[l][i+1][trans]:=dp[l][i+1][trans]+dp[l][i][j]; continue; end; if (plugleft=1)and(plugup=1) then begin sum:=0; for k:=i+2 to m+1 do begin trans:=checkbyte_three(j,k); if trans=1 then inc(sum) else if trans=2 then if sum=0 then break else dec(sum); end; if k=m+1 then continue; trans:=j-th[k]-th[i]-th[i+1]; dp[l][i+1][trans]:=dp[l][i+1][trans]+dp[l][i][j]; continue; end; if (plugleft=2)and(plugup=2) then begin sum:=0; for k:=i-1 downto -1 do begin trans:=checkbyte_three(j,k); if trans=2 then inc(sum) else if trans=1 then if sum=0 then break else dec(sum); end; if k=-1 then continue; trans:=j-2*th[i]-2*th[i+1]+th[k]; dp[l][i+1][trans]:=dp[l][i+1][trans]+dp[l][i][j]; continue; end; end; end; writeln(dp[n-1][m][1+2*th[m]]); readln(n,m); until n=0; end.