【noi2009 植物大战僵尸】最大权闭合图、网络流

这题模型还真是。。裸

连我这种网络流无能力者也在几乎刚读完题就想到了最大权闭合图的模型。。

题目中的植物就是点，权值就是score，依赖关系就是每个点依赖于它右边那个点和能罩着它的点

好了，最大权闭合图的具体搞法参考AMBER大神的论文，这题就搞得差不多了

还有一个要注意的地方，这题中的依赖关系会成环，可以想到如果成环了，那这个环上的所有点都选不了了

一开始我想丑了，其实拓扑排序一遍，把剩下的点删掉就可以了（就是把这个点权值改为负无穷就可以了）

我的程序就裸找环，O（点数*边数），至于时间嘛，丑了点，但还是过了

写的有点丑，海涵海涵。。

 

囧，下面的代码实际上是不对的

判环那个地方我的做法是错误的，正确的做法不仅简洁而且好的多。。。

这里实际上只要拓扑排序一边去掉不能拓扑出来的点就可以了~~

因为拓扑不出来的点一定是打不掉的

 

代码：

program syj; {netflow} const maxn=85;maxm=85;oo=maxlongint shr 1; var n,m,tot,s,t,ans,i,j,k,x,y,e:longint; h,d,fir,pre:array[1..maxn*maxm]of longint; v,vv,o:array[1..maxn*maxm]of boolean; q,next,point,bak,w:array[1..maxn*maxm*maxm+maxn*maxm*2]of longint; a,b:array[1..maxn,1..maxm]of longint; f:array[1..maxn*maxm,1..maxn*maxm]of boolean; procedure find(fa,i:longint); var j,k:longint; begin j:=h[i]; while j<>0 do begin if w[j]=oo then begin k:=point[j]; if not v[k] then begin v[k]:=true;f[fa,k]:=true; if k=fa then exit; find(fa,k); if f[fa,fa] then exit; end; end; j:=next[j]; end; end; procedure imp; var i,j:longint; begin i:=pre[t];j:=oo; while i<>0 do begin if w[fir[i]]<=j then begin j:=w[fir[i]];x:=i; end; i:=pre[i]; end; i:=pre[t]; while i<>0 do begin dec(w[fir[i]],j); inc(w[bak[fir[i]]],j); i:=pre[i]; end; dec(ans,j); end; procedure dfs; var i,j,k:longint; ok:boolean; begin for i:=1 to tot do fir[i]:=h[i]; pre[s]:=0; i:=s; while i<>0 do begin j:=fir[i]; ok:=false; while j<>0 do begin k:=point[j]; if not o[k] then begin if (w[j]>0)and(d[k]=d[i]+1) then begin pre[k]:=i; fir[i]:=j; ok:=true; i:=k; if k=t then begin imp;i:=x; end; break; end; end; j:=next[j]; end; if not ok then begin d[i]:=-1;i:=pre[i]; end; end; end; function bfs:boolean; var i,j,k,st,ed:longint; begin fillchar(d,sizeof(d),$FF); st:=0;ed:=1;q[1]:=s;d[s]:=0; while st<ed do begin inc(st);i:=q[st];j:=h[i]; while j<>0 do begin k:=point[j]; if not o[k] then begin if (w[j]>0)and(d[k]=-1) then begin d[k]:=d[i]+1; inc(ed);q[ed]:=k; if k=t then exit(true); end; end; j:=next[j]; end; end; bfs:=false; end; procedure link(x,y,z:longint); begin inc(e); next[e]:=h[x];point[e]:=y;w[e]:=z;h[x]:=e; inc(e); next[e]:=h[y];point[e]:=x;w[e]:=0;h[y]:=e; bak[e]:=e-1;bak[e-1]:=e; end; begin assign(input,'pvz.in');reset(input); assign(output,'pvz.out');rewrite(output); readln(n,m); tot:=0; for i:=1 to n do for j:=1 to m do begin inc(tot); a[i,j]:=tot; end; s:=tot+1; t:=tot+2; tot:=tot+2; for i:=1 to n do for j:=1 to m do begin read(b[i,j]); x:=b[i,j]; read(x); for k:=1 to x do begin read(x,y); inc(x);inc(y);link(a[x,y],a[i,j],oo); end; if j<m then link(a[i,j],a[i,j+1],oo); end; for i:=1 to n*m do begin fillchar(v,sizeof(v),0); find(i,i); end; for i:=1 to n do for j:=1 to m do if f[a[i,j],a[i,j]] then b[i,j]:=-oo; for i:=1 to n do for j:=1 to m do begin x:=b[i,j]; if x>0 then begin inc(ans,x); link(s,a[i,j],x); end else link(a[i,j],t,-x); end; while bfs do dfs; if ans<0 then ans:=0; writeln(ans); close(input); close(output); end.

最慢个点要1.5秒。。囧