ZJOI2010网络扩容

无限orz hzwer神牛……

第一问很简单，按数据建图，然后一遍最大流算法即可。

    第二问则需要用最小费用最大流算法，主要是建图，那么可以从第一问的残留网络上继续建图，对残留网络上的每一条边建一条容量是∞费用是w的边（反向弧容量为0，费用为-w），然后建一个超级源点，从超级源向1建一条容量为k，费用为0的边，对这个图进行最小费用最大流算法。

    最小费用最大流操作：

    1.首先要对于这道题的图来说，有的边需要花费费用，而有的又不用，而不用扩容的边费用为0，需要扩容的边费用为w，容量无限，这就是本题这样建图的原因。

    2.对于残留网络进行费用最短路SPFA算法，不用扩容的边一定会选费用为0的边，然后记录路径，找最小容量对可行路进行增流，更新ans

——hzwer

 uses math;
 const inf=maxlongint;
 type node=record
      from,go,v,c,t,next:longint;
      end;
 var  i,n,m,k,t,u,v,w,c,s,cnt,ans:longint;
      q,first,from,h,cur,d:array[..] of longint;
      e:array[..] of node;
 procedure ins(u,v,w,c:longint);
  begin
    inc(cnt);
    e[cnt].go:=v;e[cnt].from:=u;
    e[cnt].v:=w;e[cnt].t:=c;
    e[cnt].next:=first[u];first[u]:=cnt;
  end;
 procedure insert(u,v,w,c:longint);
  begin
    ins(u,v,w,c);ins(v,u,,-c);
  end;
 procedure ins2(u,v,w,c:longint);
  begin
    inc(cnt);
    e[cnt].go:=v;e[cnt].from:=u;
    e[cnt].v:=w;e[cnt].c:=c;
    e[cnt].next:=first[u];first[u]:=cnt;
  end;
 procedure insert2(u,v,w,c:longint);
  begin
    ins2(u,v,w,c);ins2(v,u,,-c);
  end;
 function bfs:boolean;
  var head,tail,i,x,y:longint;
    begin
      head:=;tail:=;fillchar(h,sizeof(h),);
      q[]:=s;h[s]:=;
      while head<tail do
       begin
         inc(head);
         x:=q[head];
         i:=first[x];
         while i<> do
          begin
            y:=e[i].go;
            if (e[i].v<>) and (h[y]=) then
             begin
               h[y]:=h[x]+;
               inc(tail);
               q[tail]:=y;
             end;
            i:=e[i].next;
          end;
       end;
    exit(h[t]<>);
    end;
 function dfs(x,f:longint):longint;
  var  i,tmp,used,y:longint;
    begin
    if x=t then exit(f);
    tmp:=;used:=;
    i:=cur[x];
    while i<> do
      begin
       y:=e[i].go;
       if (e[i].v<>) and (h[y]=h[x]+) then
        begin
          tmp:=dfs(y,min(f-used,e[i].v));
          dec(e[i].v,tmp);
          inc(e[i xor ].v,tmp);
          inc(used,tmp);
          if e[i].v<> then cur[x]:=i;
          if used=f then exit(f);
        end;
       i:=e[i].next;
      end;
    if used= then h[x]:=-;
    exit(used);
    end;
 procedure dinic;
  begin
    while bfs do
     begin
       for i:= to n do cur[i]:=first[i];
       inc(ans,dfs(s,inf));
     end;
  end;
 procedure build;
  var tmp:longint;
  begin
  tmp:=cnt;
  for i:= to cnt do
    if i and = then insert2(e[i].from,e[i].go,inf,e[i].t);
  end;
 function spfa:boolean;
  var i,x,y,head,tail:longint;
       v:array[..] of boolean;
    begin
     head:=;tail:=;
     for i:= to n do d[i]:=inf;
     fillchar(v,sizeof(v),false);
     q[]:=;d[]:=;v[i]:=true;
     while head<tail do
      begin
       inc(head);
       x:=q[head]; v[x]:=false;
       i:=first[x];
       while i<> do
        begin
         y:=e[i].go;
         if (e[i].v>) and (d[x]+e[i].c<d[y]) then
          begin
            d[y]:=d[x]+e[i].c;
            from[y]:=i;
            if not(v[y]) then
             begin
               inc(tail);
               q[tail]:=y;
               v[y]:=true;
             end;
          end;
         i:=e[i].next;
        end;
      end;
    exit(d[n]<>inf);
    end;
 procedure mcf;
  var i,x:longint;
    begin
    x:=inf;
    i:=from[n];
    while i<> do
      begin
       x:=min(x,e[i].v);
       i:=from[e[i].from];
      end;
    i:=from[n];
    while i<> do
      begin
       dec(e[i].v,x);
       inc(e[i xor ].v,x);
       i:=from[e[i].from];
      end;
    inc(ans,x*d[n]);
    end;
 procedure main;
  begin
    cnt:=;
    readln(n,m,k);
    for i:= to m do
     begin
       readln(u,v,w,c);
       insert(u,v,w,c);
     end;
    ans:=;
    s:=;t:=n;
    dinic;
    write(ans,' ');
    ans:=;
    build;
    ins(,,k,);
    while spfa do mcf;
    writeln(ans);
  end;
 begin
   main;
 end.