【BZOJ 1834】 [ZJOI2010]network 网络扩容

Description

给定一张有向图，每条边都有一个容量C和一个扩容费用W。这里扩容费用是指将容量扩大1所需的费用。求： 1、 在不扩容的情况下，1到N的最大流； 2、 将1到N的最大流增加K所需的最小扩容费用。

Input

输入文件的第一行包含三个整数N,M,K，表示有向图的点数、边数以及所需要增加的流量。 接下来的M行每行包含四个整数u,v,C,W，表示一条从u到v，容量为C，扩容费用为W的边。

Output

输出文件一行包含两个整数，分别表示问题1和问题2的答案。

Sample Input

5 8 2
1 2 5 8
2 5 9 9
5 1 6 2
5 1 1 8
1 2 8 7
2 5 4 9
1 2 1 1
1 4 2 1

Sample Output

13 19
30%的数据中，N<=100
100%的数据中，N<=1000，M<=5000，K<=10

 

第一问最大流

第二问

建立超级源和超级汇，超级源向1连边，容量为K，n向超级汇连边，容量同上。原来的每条边再建立一条，容量为inf，费用为W，在原图上增广，跑最小费用最大流

但是！！！新边一定要在求完第一问之后再建。。。

代码量略大。。。不过我都上套板子。。。

 #include<cstdio>
 #include<cstring>
 #include<cmath>
 #include<iostream>
 using namespace std;
 const int N=,inf=;
 struct ee{int to,next,f,w;}e[N*];
 int S,T,cnt=,n,k,ans1,ans2,timer,m,u[N],v[N],w[N],c[N],mid;
 int head[N],dis[N],pre[N],q[N];
 bool inq[N],flag=;
 void ins(int u,int v,int f,int w){
     e[++cnt].to=v,e[cnt].next=head[u],e[cnt].f=f,e[cnt].w=w,head[u]=cnt;
     e[++cnt].to=u,e[cnt].next=head[v],e[cnt].f=,e[cnt].w=-w,head[v]=cnt;
 }
 bool spfa(){
     for (int i=;i<=T;i++) dis[i]=inf;
     int h=,t=;
     q[t]=S;dis[S]=;inq[S]=;
     while (h!=t){
         int now=q[++h];if(h==) h=;
         for (int i=head[now];i;i=e[i].next){
             int v=e[i].to;
             if (dis[v]>dis[now]+e[i].w&&e[i].f){
                 dis[v]=dis[now]+e[i].w;
                 pre[v]=i;
                 if (!inq[v]){
                     q[++t]=v;if (t==) t=;
                     inq[v]=;
                 }
             }
         }
         inq[now]=;
     }
     if (dis[T]==inf) return ;
     return ;
 }

 void updata(){
     int tmp=T,flow=inf;
     while (tmp!=S){
         int l=pre[tmp],v=e[l].to;
         flow=min(flow,e[l].f);
         tmp=e[l^].to;
     }
     tmp=T;
     while (tmp!=S){
         int l=pre[tmp],v=e[l].to;
         e[l].f-=flow;e[l^].f+=flow;
         tmp=e[l^].to;
     }
     mid+=flow;
     if(mid>=k) flag=;
     ans2+=dis[T]*flow;
 }

 bool bfs(){
     for (int i=;i<=T;i++) dis[i]=inf;
     int h=,t=,now;
     q[]=;dis[]=;
     while(h!=t){
         now=q[++h];
         for (int i=head[now];i;i=e[i].next){
             int v=e[i].to;
             if (e[i].f&&dis[now]+<dis[v]){
                 dis[v]=dis[now]+;
                 if (v==n)return ;
                 q[++t]=v;
             }
         }
     }
     if (dis[n]==inf) return ; return ;
 }

 int dinic(int now,int f){
     if (now==n) return f;
     int rest=f;
     for (int i=head[now];i;i=e[i].next){
         int v=e[i].to;
         if (e[i].f&&dis[v]==dis[now]+&&rest){
             int t=dinic(v,min(rest,e[i].f));
             if (!t) dis[v]=;
             e[i].f-=t;
             e[i^].f+=t;
             rest-=t;
             //if(t) printf("%d %d %d\n",now,v,e[i].f);
         }
     }
     return f-rest;
 }

 int main(){
     scanf("%d%d%d",&n,&m,&k);
     S=,T=n+;
     for (int i=;i<=m;i++){
         scanf("%d%d%d%d",&u[i],&v[i],&c[i],&w[i]);
         ins(u[i],v[i],c[i],);
     }
     ins(S,,k,);ins(n,T,k,);
     while(bfs())
     ans1+=dinic(,inf);
     for (int i=;i<=m;i++)ins(u[i],v[i],inf,w[i]);
     while(flag&&spfa())
     updata();
     printf("%d %d",ans1,ans2);
 }