pku3084 最小割

恩，又是网络流，一个sap完事。。

发现网络流这东西确实很厉害，很多东西很多模型都需要用到它。。

题目大意：

  一堆房子，有的房子之间有门，有的房子有侵入者，有的房子没有侵入者，有一个房间需要保护。I代表没有侵入者，NI代表有侵入者，一开始的N,M代表的n个房间和m房间需要保护。NI或者I后边跟一个k，然后后边k个数表示从第i个房间可以到达这些个房间不需要钥匙，而这k个房间到达第i个房间只能到达一次或者从第i个房间过去再过来。。。

简略题解：

  仔细一看就是一个最小割。

  需要保护的房间m就是汇点，新增源点s，如果这个房间没有侵入者，那么从源点向这个点连一条流量为inf的边。对于其他的门，如果是i到j的门，那么从i连一条到j的流量为inf的边，从j到i连一条流量为1的边，这样就确保了从j到i只能进去一次或者进来再出去。

  然后就是个最大流就可以了，如果最后的答案比inf还要大，那么肯定存在一条直接通往被保护房间的路径，那么无解。否则答案就是最小割。

附程序：

program pku_3084; const inf=10000000; var pre,other,c:array[0..10000] of longint; last,d,num:array[0..22] of longint; n,testcase,i,l,st,ed:longint; {---------------------------} procedure grow(x,y,z:longint); begin inc(l); pre[l]:=last[x];last[x]:=l;other[l]:=y;c[l]:=z; inc(l); pre[l]:=last[y];last[y]:=l;other[l]:=x;c[l]:=0; end; {---------------------------} procedure init; var i,j,k,p:longint; flag:boolean; t:char; begin l:=0; fillchar(last,sizeof(last),0); readln(n,ed); st:=n+1;inc(ed); for i:=1 to n do begin read(t); flag:=true; if t<>'I' then begin read(t);flag:=false; end; read(k); for j:=1 to k do begin read(p); inc(p); grow(i,p,inf); grow(p,i,1); end; if flag then grow(st,i,inf); readln; end; end; {---------------------------} function min(x,y:longint):longint; begin if x<y then exit(x) else exit(y); end; {---------------------------} function dfs(u,flow:longint):longint; var i,t:longint; begin if u=ed then exit(flow); dfs:=0; i:=last[u]; while i<>0 do begin if (c[i]>0) and (d[u]=d[other[i]]+1) then begin t:=dfs(other[i],min(flow-dfs,c[i])); dec(c[i],t); if odd(i) then inc(c[i+1]) else inc(c[i-1]); inc(dfs,t); if dfs=flow then exit(dfs); end; i:=pre[i]; end; if d[st]>=st then exit; dec(num[d[u]]); if num[d[u]]=0 then d[st]:=st; inc(d[u]); inc(num[d[u]]); end; {---------------------------} procedure main; var flow,ans:longint; begin ans:=0; fillchar(d,sizeof(d),0); fillchar(num,sizeof(num),0); num[0]:=st; while d[st]<st do begin flow:=dfs(st,inf); inc(ans,flow); end; if ans>=inf then writeln('PANIC ROOM BREACH') else writeln(ans); end; {---------------------------} begin readln(testcase); for i:=1 to testcase do begin init; main; end; end.