[PKU 2125]Destroying the Graph(二分图最小权覆盖集)

【题目大意】：

给你一个有向图，现在要你删除所有的边。可以删除每个点的所有出边或是入边，都有他们相应的费用，求最小费用和方案。

【题目分析】：

首先呢，这个东西的想法就是要参照最小路径覆盖，将点拆成出点和入点，中间连边。这样我们发现实际上我们再加个源和汇，连接出入点的容量为成本的话，这个东西的简单割就是一个方案。而我们要最小的，那就是最小割~~真是个悲剧的构图法……

这个模型的名字叫做二分图的最小权覆盖集。对于二分图最大权独立集，这个就是一个对偶问题。用总的减下就好~~

P.s. 万恶的Special Judge，真的很悲剧……

【代码】：

program PKU_2125; const inf=10000000; type ntype=record data,pre,c,op,next:longint; end; var g,dis,vh:array[0..300] of longint; v:array[0..300] of boolean; node:array[0..50000] of ntype; ou:array[0..300,1..2] of longint; ans,s,t,x,y,nn,i,n,m,tot,all:longint; function min2(x,y:longint):longint; begin if x<y then exit(x) else exit(y); end; procedure link(x,y,w:longint); begin inc(tot); node[tot].pre:=x; node[tot].data:=y; node[tot].c:=w; node[tot].op:=tot+1; node[tot].next:=g[x]; g[x]:=tot; inc(tot); node[tot].pre:=y; node[tot].data:=x; node[tot].c:=0; node[tot].op:=tot-1; node[tot].next:=g[y]; g[y]:=tot; end; function aug(x,flow:longint):longint; var p,bak,min,tmp:longint; begin if x=t then exit(flow); bak:=flow; min:=nn-1; p:=g[x]; while p<>0 do begin if node[p].c>0 then begin if dis[node[p].data]+1=dis[x] then begin tmp:=aug(node[p].data,min2(bak,node[p].c)); dec(bak,tmp); dec(node[p].c,tmp); inc(node[node[p].op].c,tmp); if dis[s]>=nn then exit(flow-bak); if bak=0 then break; end; if dis[node[p].data]<min then min:=dis[node[p].data]; end; p:=node[p].next; end; if bak=flow then begin dec(vh[dis[x]]); if vh[dis[x]]=0 then dis[s]:=nn; dis[x]:=min+1; inc(vh[dis[x]]); end; exit(flow-bak); end; procedure dfs(x:longint); var p:longint; begin v[x]:=true; p:=g[x]; while p<>0 do begin if (node[p].c>0) and not v[node[p].data] then dfs(node[p].data); p:=node[p].next; end; end; begin assign(input,'a.in'); assign(output,'a.out'); reset(input); rewrite(output); readln(n,m); s:=0; t:=n*2+1; nn:=t+1; for i:=1 to n do begin read(x); link(i+n,t,x); end; for i:=1 to n do begin read(x); link(s,i,x); end; for i:=1 to m do begin readln(x,y); link(x,y+n,inf); end; vh[0]:=nn; ans:=0; while dis[s]<nn do inc(ans,aug(s,inf)); dfs(s); all:=0; for i:=1 to tot do if (v[node[i].pre]<>v[node[i].data]) and odd(i) then if node[i].pre=s then begin inc(all); ou[all,1]:=node[i].data; ou[all,2]:=-1; end else if node[i].data=t then begin inc(all); ou[all,1]:=node[i].pre-n; ou[all,2]:=1; end; writeln(ans); writeln(all); for i:=1 to all do begin write(ou[i,1],' '); if ou[i,2]=1 then writeln('+') else writeln('-'); end; close(input); close(output); end.