【四川04年省选day2——洪水】状态压缩dp

题目：

乡下有n(n<=8)个小村庄，用m条双向道路连接起来，每两个村庄最多用一条道路直接连接。最近村庄闹涝灾，洪水让这些道路各有一定的概率被淹没，而这些道路被淹没与否是相互独立的。你的任务是计算出所有村庄连在一起(即：从任意一个村庄都可以经过未被淹没的道路到达其他所有村庄)的概率有多大。

分析：n=8，显然不是裸搜就是状态压缩dp，裸搜显然也是很难过的，那么本题的方向就出来了

算法：

   记f[i]表示i这个状态中所有的点连通的概率，如果直接把要选的点分成2个集合什么的做法是会算重的

   而为了避免算重，我们需要规定一些东西，比方说可以这样规定：先孤立出所选集合中的第一个点，然后枚举其他点的连通情况，且其他点形成的连通块之间不能直接相连，只能和这个孤立出的点相连，如图

里面需要套一个搜其他点连通情况的dfs，惊奇发现大过程套小过程，小过程还可以调用大过程，而且大过程的局部变量是小过程的全程变量，只是debug的时候监控不了，使我郁闷好久，考完才想到解决办法。。

大部分概率计算都比较飘，这题也不例外，我差不多花了一个多小时才晓得样例是怎么搞的。。。

这题数据规模真是小，我这个O(2^n*n!*n^2)居然全0ms啦，哈哈。

program syj; {dp} var n,m,all,i,j,k,k1,k2:longint; a:array[1..8,1..8]of real; f:array[0..1<<8]of real; bb:array[0..8]of real; g:array[1..8,1..8]of boolean; procedure dfs(i:longint); var jj:longint; b:array[0..8]of longint; procedure find(j,k:longint); var l,ll,lj,lk:longint; t1,t2:real; begin if j=0 then begin t1:=1; for l:=1 to b[0] do bb[l]:=b[l]; for l:=1 to b[0] do begin dfs(b[l]); j:=b[l];ll:=1;t2:=1; while j>0 do begin if odd(j) then t2:=t2*a[jj,ll]; j:=j>>1;inc(ll); end; t1:=t1*(1-t2)*f[b[l]]; end; for l:=1 to b[0]-1 do for ll:=l+1 to b[0] do begin j:=b[l];lj:=1;t2:=1; while j>0 do begin if odd(j) then begin k:=b[ll];lk:=1; while k>0 do begin if odd(k) then t2:=t2*a[lj,lk]; k:=k>>1;inc(lk); end; end; j:=j>>1;inc(lj); end; t1:=t1*t2; end; f[i]:=f[i]+t1; exit; end; while not odd(j) do begin j:=j>>1;inc(k); end; for l:=1 to b[0] do begin ll:=b[l];b[l]:=b[l] or(1<<(k-1)); find(j>>1,k+1); b[l]:=ll; end; inc(b[0]);b[b[0]]:=1<<(k-1); find(j>>1,k+1); b[b[0]]:=0;dec(b[0]); end; begin if f[i]>-1 then exit; f[i]:=0; for jj:=1 to n do if odd(i>>(jj-1)) then break; b[0]:=0; find(i xor(1<<(jj-1)),1); end; begin assign(input,'flood.in');reset(input); assign(output,'flood.out');rewrite(output); readln(n,m); all:=(1<<n)-1; for i:=1 to n do for j:=1 to n do a[i,j]:=1; for i:=1 to m do begin readln(j,k,a[j,k]);a[k,j]:=a[j,k]; // a[j,k]:=1-a[j,k];a[k,j]:=1-a[k,j]; end; for i:=1 to all do f[i]:=-2; for i:=1 to n do f[1<<(i-1)]:=1; dfs(all); {for i:=1 to n do for j:=i+1 to n do if a[i,j]>1e-8 then for k1:=1 to all-1 do if odd(k1>>(i-1)) then for k2:=1 to all-1 do if odd(k2>>(j-1))and(k1 and k2=0) then f[k1 or k2]:=f[k1 or k2]+f[k1]*f[k2]*a[i,j];} writeln(f[all]:0:3); close(input); close(output); end.