[ZJU 2676]Network Wars(分数规划+最小割)

【题目大意】：

求无向图边权平均值最小的割集~

【题目分析】：

一看到平均值就应该敏锐的觉察到与分数规划有关，这个问题也不例外。

这还是一个0/1分数规划问题，因为每个边只能取或是不取（废话……）

然后东西就和上一篇的东西一样，还是把边权改下，然后求最小割。

这里有几个小技巧需要注意：

1、改完了的负权边直接进割集，因为这样一定更小~~（P.s.分数规划这点就是悲剧，就是二分完了答案之后还不算，还要最优~）

2、网络流不必要拆点，看了cai0715神牛拆了点，但是后来发现没有必要。

3、最好每一次都重建图，避免错误

4、网络流的算法中只有加减，所以坚决不能考虑精度问题！！

5、改用精度的东西也要用，东西要开就开double，real是残废……

实际上写出来并不是很麻烦，就是需要注意的地方很多，又是一个浪费了我一下午的题啊……

【代码】：

program ZJU_2676; const inf=1e8; type ntype=record data,op,next:longint; c:double; end; etype=record x,y,w:longint; end; var dis,vh,g:array[0..101] of longint; edges:array[1..2000] of etype; v,vis:array[0..2000] of boolean; node:array[0..2000] of ntype; ou:array[0..2000] of longint; n,m,i,x,y,w,tot,all:longint; s,t,mid,ans,now:double; function min2(x,y:double):double; begin if x<y then exit(x) else exit(y); end; function max2(x,y:double):double; begin if x>y then exit(x) else exit(y); end; procedure link(x,y:longint;w:double); begin inc(tot); 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].data:=x; node[tot].c:=0; node[tot].op:=tot-1; node[tot].next:=g[y]; g[y]:=tot; end; procedure init(x:double); var i:longint; begin now:=0; fillchar(node,sizeof(node),0); fillchar(g,sizeof(g),0); fillchar(v,sizeof(v),0); tot:=0; for i:=1 to m do if edges[i].w<=x then begin now:=now+edges[i].w-x; v[i]:=true; end else begin link(edges[i].x,edges[i].y,edges[i].w-x); link(edges[i].y,edges[i].x,edges[i].w-x); end; end; function aug(x:longint;flow:double):double; var bak,tmp:double; min,p:longint; begin if x=n then exit(flow); min:=n-1; bak:=flow; 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)); node[p].c:=node[p].c-tmp; node[node[p].op].c:=node[node[p].op].c+tmp; bak:=bak-tmp; if dis[1]>=n 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=0 then begin dec(vh[dis[x]]); if vh[dis[x]]=0 then dis[1]:=n; dis[x]:=min+1; inc(vh[dis[x]]); end; exit(flow-bak); end; procedure dfs(x:longint); var p:longint; begin vis[x]:=true; p:=g[x]; while p<>0 do begin if not vis[node[p].data] and (node[p].c>0) then dfs(node[p].data); p:=node[p].next; end; end; procedure print; begin fillchar(vis,sizeof(vis),0); dfs(1); for i:=1 to m do if vis[edges[i].x]<>vis[edges[i].y] then v[i]:=true; all:=0; for i:=1 to m do if v[i] then begin inc(all); ou[all]:=i; end; writeln(all); for i:=1 to all do write(ou[i],' '); writeln; writeln; end; begin assign(input,'a.in'); assign(output,'a.out'); reset(input); rewrite(output); while not eof do begin readln(n,m); fillchar(g,sizeof(g),0); fillchar(ou,sizeof(ou),0); fillchar(edges,sizeof(edges),0); s:=1e7; t:=1; for i:=1 to m do begin readln(x,y,w); edges[i].x:=x; edges[i].y:=y; edges[i].w:=w; s:=min2(s,w); t:=max2(t,w); end; while s<t do begin mid:=(s+t)/2; init(mid); ans:=0; fillchar(dis,sizeof(dis),0); fillchar(vh,sizeof(vh),0); vh[0]:=n; while dis[1]<n do ans:=ans+aug(1,inf); if abs(ans+now)<1e-7 then begin print; break; end; if ans+now>1e-7 then s:=mid else t:=mid; end; end; close(input); close(output); end.