【POJ1739、travel、画圈圈】连通性状态压缩DP

 

其实这所谓连通性状态压缩dp并不麻烦，有模板可循

 

算法流程

1.       预处理出所有有效状态，这个有效状态的定义是所有括号是配对的，即括号序列合法，并给每个状态一个唯一对应的哈希值（3进制），我用0表示无插头，1表示左插头，-1表示右插头

2.       预处理出所有转移，形如g[i,j,k]，表示轮廓线的状态为i，转折点为j，第k号转移到哪个状态（编号，不记哈希值）。一般是以下几种转移（但对于每个特定的格子只有两种转移）

0 0 –> 0 0 or 1 -1 空格不管或在空格新建两个插头

0 1 –> 0 1 or 1 0 将一边的插头延续，2个方向均可

0 -1 -> 0 -1 or -1 0将一边的插头延续，2个方向均可

1 0 -> 1 0 or 0 1将一边的插头延续，2个方向均可

-1 0 -> -1 0 or 0 -1将一边的插头延续，2个方向均可

1 -1 -> 0 0 将一个连通块合并，“收圈”

-1 1 -> 0 0 合拢两条挤到一起的“线”

1 1 … -1 -1 -> 0 0 … 1 -1合拢两条挤到一起的“线”

1 1 … -1 -1 -> 1 -1 … 0 0合拢两条挤到一起的“线”

需要注意的是，找某状态对应的编号应用hash，不能裸找，裸找能慢十几倍

注：图片来自盾哥博客

3.       搞完了预处理，dp就容易多了，一般是根据题目看某些特殊的格子只能有特定的转移，讨论一下即可。这里的空间消耗一般会比较大，所以一般得写滚动数组

 

有了这些预处理，dp就方便多了，直接写就可以了

 

【POJ1739】

 

直接用连通性状态压缩dp求哈密尔顿路径有点麻烦，所以可以看成是求从(n,1)->(n,m)再从(n,m)->(n+1,m)->(n+1,1)->(n,1)的一条哈密尔顿环即可

具体的话相当于直接dp，ans=f[n+1,1,第二个有左括号，第m+1有右括号的状态]

program syj; const maxs=1000; maxn=9; type state=array[1..maxn]of longint; var n,m,i,j,k,l,r,u,t,old,now,z:longint; ok:boolean; a,b:state; cc:char; s:array[0..maxs]of state; h:array[0..20000]of longint; g,gg:array[0..maxs,0..maxn,0..1]of longint; f:array[0..1,0..maxs]of longint; function hash:longint; var i:longint; begin hash:=0; for i:=1 to m+1 do hash:=hash*3+a[i]+1; end; procedure dfs(i,k:longint); begin if i>m+1 then begin if k=0 then begin inc(t); s[t]:=a; h[hash]:=t; end; exit; end; a[i]:=0; dfs(i+1,k); a[i]:=1; dfs(i+1,k+1); if k>0 then begin a[i]:=-1; dfs(i+1,k-1); end; end; procedure find(kk:longint); begin if j=m then begin if a[m+1]<>0 then exit; move(a[1],a[2],m*4); a[1]:=0; end; g[i,j,kk]:=h[hash]; gg[i,j,kk]:=z; end; begin assign(input,'pku1739.in');reset(input); assign(output,'pku1739.out');rewrite(output); readln(n,m); while n>0 do begin fillchar(g,sizeof(g),0); fillchar(h,sizeof(h),0); t:=0; dfs(1,0); for i:=1 to t do begin b:=s[i]; for j:=1 to m do begin if (b[j]=0)and(b[j+1]=0) then begin a:=b; z:=1; find(0); z:=0; a[j]:=1;a[j+1]:=-1; find(1); end else if b[j]=0 then begin a:=b; find(0); a[j]:=a[j+1];a[j+1]:=0; find(1); end else if b[j+1]=0 then begin a:=b; find(0); a[j+1]:=a[j];a[j]:=0; find(1); end else if (b[j]=-1)and(b[j+1]=1) then begin a:=b; a[j]:=0;a[j+1]:=0; find(0); end else if (b[j]=1)and(b[j+1]=1) then begin a:=b; a[j]:=0;a[j+1]:=0; k:=1; for l:=j+2 to m+1 do begin inc(k,a[l]); if k=0 then break; end; a[l]:=1; find(0); end else if (b[j]=-1)and(b[j+1]=-1) then begin a:=b; a[j]:=0;a[j+1]:=0; k:=-1; for l:=j-1 downto 1 do begin inc(k,a[l]); if k=0 then break; end; a[l]:=-1; find(0); end; end; end; fillchar(f,sizeof(f),0); f[0,1]:=1; old:=0;now:=1; for i:=1 to n do begin for j:=1 to m do begin read(cc); fillchar(f[now],sizeof(f[now]),0); for k:=1 to t do for r:=0 to 1 do if g[k,j,r]>0 then if not((cc='#')and(gg[k,j,r]=0)or(cc='.')and(gg[k,j,r]=1)) then inc(f[now,g[k,j,r]],f[old,k]); old:=now;now:=1-now; end; readln; end; fillchar(a,sizeof(a),0); a[2]:=1;a[m+1]:=-1; for k:=1 to t do begin ok:=true; for l:=1 to m+1 do if s[k,l]<>a[l] then begin ok:=false;break; end; if ok then break; end; writeln(f[old,k]); readln(n,m); end; close(input);close(output); end.  

 

travel

题意：

在N*M的方格上，求(a,b) -> (c,d)的哈密尔顿路径条数

算法：

在上面的基础上加入独立插头

 

用4进制表示，论文中没有说，这里补充一下，关于插头的准确定义

独立插头：往上延伸没有回到轮廓线上

左括号：往上延伸并回到轮廓线上，且对应点在这个位置的右边

右括号：与左括号相反

主要有以下几种转移

只有起点和终点允许的：

1．两边都没有插头 -> 往任一方向开一独立插头

2.只有一边有括号 –> 把这个括号消掉，并把这个括号原来对应的括号改成独立插头

3.只有一个独立插头过来 –>消掉这个独立插头

非起点终点或(n,m)允许的：

一边是括号，一边是独立插头 -> 合并这个括号和独立插头，同时把与那个括号配对的括号改成独立插头

(n,m)非起点终点时只有(n,m) 允许的：

合并两个独立插头

 

另外，所有点都不允许合并左右插头

所以这样就可以了

代码：

 

program syj; const maxs=9000; maxm=9; type state=array[1..maxm]of longint; var i,j,k,l,t,n,m,old,now,sx,sy,tx,ty:longint; a,b:state; s:array[0..maxs]of state; h:array[0..1 shl (maxm*2+1)]of longint; g:array[0..maxs,0..maxm,0..16]of longint; mp:array[0..maxm,0..maxm]of longint; o:array[0..2,0..16]of boolean; f:array[0..1,0..maxs]of longint; function hash:longint; var i:longint; begin hash:=0; for i:=1 to m+1 do hash:=hash*4+a[i]+3; end; procedure dfs(i,k,l:longint); begin if i>m+1 then begin if k=0 then begin inc(t); s[t]:=a; h[hash]:=t; end; exit; end; a[i]:=0; dfs(i+1,k,l); a[i]:=1; dfs(i+1,k+1,l); if k>0 then begin a[i]:=-1; dfs(i+1,k-1,l); end; if l<2 then begin a[i]:=2; dfs(i+1,k,l+1); end; end; procedure find(kk:longint); var k:longint; begin if j=m then begin if a[m+1]<>0 then exit; move(a[1],a[2],m*4); a[1]:=0; end; k:=h[hash]; g[i,j,kk]:=k; end; function pp(i,z:longint):longint; var k:longint; begin k:=z; while k<>0 do begin i:=i+z; if a[i]<2 then k:=k+a[i]; end; pp:=i; end; begin assign(input,'travel.in');reset(input); assign(output,'travel.out');rewrite(output); readln(n,m); dfs(1,0,0); for i:=1 to t do begin b:=s[i]; for j:=1 to m do begin if (b[j]=0)and(b[j+1]=0) then begin a:=b; a[j]:=1;a[j+1]:=-1; find(0); a:=b; a[j]:=0;a[j+1]:=2; find(1); a:=b; a[j]:=2;a[j+1]:=0; find(2); end else if b[j]=0 then begin a:=b; find(3); a:=b; a[j]:=a[j+1];a[j+1]:=0; find(4); if b[j+1]<>2 then begin a:=b; a[j]:=b[j];a[j+1]:=0; a[pp(j+1,b[j+1])]:=2; find(5); end else begin a:=b; a[j]:=0;a[j+1]:=0; find(12); end; end else if b[j+1]=0 then begin a:=b; find(6); a:=b; a[j+1]:=a[j];a[j]:=0; find(7); if b[j]<>2 then begin a:=b; a[j]:=0;a[j+1]:=b[j+1]; a[pp(j,b[j])]:=2; find(8); end else begin a:=b; a[j]:=0;a[j+1]:=0; find(13); end; end else if (b[j]=-1)and(b[j+1]=1) then begin a:=b; a[j]:=0;a[j+1]:=0; find(9); end else if (b[j]=1)and(b[j+1]=1) then begin a:=b; a[j]:=0;a[j+1]:=0; a[pp(j+1,1)]:=1; find(10); end else if (b[j]=-1)and(b[j+1]=-1) then begin a:=b; a[j]:=0;a[j+1]:=0; a[pp(j,-1)]:=-1; find(11); end else if (b[j]=2)and(b[j+1]=2) then begin a:=b; a[j]:=0;a[j+1]:=0; find(14); end else if b[j]=2 then begin a:=b; a[j]:=0;a[j+1]:=0; a[pp(j+1,b[j+1])]:=2; find(15); end else if b[j+1]=2 then begin a:=b; a[j]:=0;a[j+1]:=0; a[pp(j,b[j])]:=2; find(16); end; end; end; readln(sx,sy,tx,ty); mp[n,m]:=2; mp[sx,sy]:=1;mp[tx,ty]:=1; o[0,0]:=true;o[0,3]:=true;o[0,4]:=true;o[0,6]:=true;o[0,7]:=true; o[0,9]:=true;o[0,10]:=true;o[0,11]:=true;o[0,15]:=true;o[0,16]:=true; o[1,1]:=true;o[1,2]:=true;o[1,5]:=true;o[1,8]:=true;o[1,12]:=true; o[1,13]:=true; o[2,14]:=true; f[0,1]:=1; old:=0;now:=1; for i:=1 to n do for j:=1 to m do begin fillchar(f[now],sizeof(f[now]),0); for k:=1 to t do for l:=0 to 16 do if g[k,j,l]>0 then if o[mp[i,j],l] then inc(f[now,g[k,j,l]],f[old,k]); old:=now;now:=1-now; end; writeln(f[old,1]); close(input);close(output); end.  

ps：加入独立插头后状态数激增。。

 

【画圈圈】

今年集训队的题目，大意就略了吧，ms衡阳八中OJ上有

题目主要就是要确定每个点是圈圈内部还是外部，这里其实只要统计包含它的括号对数即可，基数为内部，偶数为外部

代码：

{$inline on} program syj; const mo=123456791; maxs=45000; type state=array[0..13]of longint; var n,m,i,j,k,l,t,old,now,z:longint; cc:char; a,b:state; ok:boolean; s:array[0..maxs]of state; h:array[0..2000000]of longint; ss:array[0..maxs,0..13]of longint; f:array[0..1,0..maxs]of longint; g,gg:array[0..maxs,0..13,0..1]of longint; function hash:longint; inline; var i:longint; begin hash:=0; for i:=1 to m+1 do hash:=hash*3+a[i]+1; end; procedure dfs(i,k:longint); inline; begin if i>m+1 then begin if k=0 then begin inc(t); j:=0; s[t]:=a; for l:=1 to m+1 do begin ss[t,l]:=ss[t,l-1]+a[l]; j:=j*3+a[l]+1; end; h[j]:=t; end; exit; end; a[i]:=0; dfs(i+1,k); a[i]:=1; dfs(i+1,k+1); if k>0 then begin a[i]:=-1; dfs(i+1,k-1); end; end; procedure find(kk:longint); inline; var k,l:longint; begin if j=m then begin if a[m+1]<>0 then exit; move(a[1],a[2],m*4); a[1]:=0; end; g[i,j,kk]:=h[hash]; gg[i,j,kk]:=z; end; procedure add(var x:longint;y:longint); inline; begin inc(x,y); if x>=mo then dec(x,mo); end; begin assign(input,'input.txt');reset(input); assign(output,'output.txt');rewrite(output); readln(n,m); dfs(1,0); for i:=1 to t do begin b:=s[i]; for j:=1 to m do if (b[j]=0)and(b[j+1]=0) then begin a:=b; z:=1; find(0); z:=0; a[j]:=1;a[j+1]:=-1; find(1); end else if b[j]=0 then begin a:=b; find(0); a[j]:=a[j+1];a[j+1]:=0; find(1); end else if b[j+1]=0 then begin a:=b; find(0); a[j+1]:=a[j];a[j]:=0; find(1); end else if b[j]<>b[j+1] then begin a:=b; a[j]:=0;a[j+1]:=0; find(0); end else if (b[j]=1)and(b[j+1]=1) then begin a:=b; a[j]:=0;a[j+1]:=0; k:=1; for l:=j+2 to m+1 do begin inc(k,a[l]); if k=0 then break; end; a[l]:=1; find(0); end else begin a:=b; a[j]:=0;a[j+1]:=0; k:=-1; for l:=j-1 downto 1 do begin inc(k,a[l]); if k=0 then break; end; a[l]:=-1; find(0); end; end; f[0,1]:=1; old:=0;now:=1; for i:=1 to n do begin for j:=1 to m do begin read(cc); fillchar(f[now],sizeof(f[now]),0); for k:=1 to t do for l:=0 to 1 do if (g[k,j,l]>0)and not( (cc='.')and(gg[k,j,l]=1) or (cc='*')and( (gg[k,j,l]=0)or(not odd(ss[k,j])) ) or (cc='#')and( (gg[k,j,l]=0)or(odd(ss[k,j])) ) ) then add(f[now,g[k,j,l]],f[old,k]); old:=now;now:=1-now; end; readln; end; writeln(f[old,1]); close(input);close(output); end.  

 

这种题目还是挺麻烦的！另外，NOI2010那道神题暂时还没弄懂，以后再说吧