bzoj1834: [ZJOI2010]network 网络扩容

1834: [ZJOI2010]network 网络扩容

Time Limit: 3 Sec  Memory Limit: 64 MB
Submit: 2677  Solved: 1359
[Submit][Status][Discuss]

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

#include<cstdio>
#include<iostream>
#include<cstring>
#define N 102333
#include<queue>
#define inf 2147483647
using namespace std;
inline int read(){
     int x=0; char ch=getchar();
     while(ch<'0'||ch>'9') ch=getchar();
     while(ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar();
     return x;
}
struct node{
     int w,to,next,c,from;
}e[N<<1],E[N<<1];
int head[N],dis[N],from[N],i,j,zs1,zs2,zs3,zs4,tot=1,t,anss,m,n,sum,flow;
bool vis[N];
int tot1,tot2,num[N][2];
inline void ins(int u,int v,int w,int cost){
     e[++tot].to=v; e[tot].next=head[u]; head[u]=tot; e[tot].w=w; e[tot].c=cost; e[tot].from=u;
}
inline bool spfa(){
     for(int i=0;i<=t;i++) dis[i]=inf; vis[0]=1; dis[0]=0; queue<int>q; q.push(0);
     while(!q.empty()) {
          int x=q.front(); q.pop();
          for(int k=head[x];k;k=e[k].next)
            if (dis[x]+e[k].c<dis[e[k].to] && e[k].w>0) {
                    from[e[k].to]=k;
                dis[e[k].to]=dis[x]+e[k].c;
                if(!vis[e[k].to]) {
                     vis[e[k].to]=1; q.push(e[k].to);
                }
            }
          vis[x]=0;
     }
     if(dis[t]==inf) return 0;else return 1;
}
inline bool bfs(){
     for(i=0;i<=t;i++) dis[i]=-1; queue<int>q; q.push(1); dis[1]=0;
     while(!q.empty()) {
          int x=q.front(); q.pop();
          for(int k=head[x];k;k=e[k].next) 
             if(dis[e[k].to]<0 && e[k].w>0) {
                   dis[e[k].to]=dis[x]+1; q.push(e[k].to);
             }
     }
     if(dis[t]>0) return 1;else return 0;
}
int find(int x,int low){
     if(x==t) return low;
     int delta=low,now;
     for(int k=head[x];k;k=e[k].next) 
       if(e[k].w>0 && dis[e[k].to]==dis[x]+1){ 
           now=find(e[k].to,min(e[k].w,delta));
           e[k].w-=now; e[k^1].w-=now;   delta-=now;
           if(!delta) return low;
        } 
     dis[x]=-1;
     return low-delta;
}
inline void run() {
     int x=inf;
     for(int i=from[t];i;i=from[e[i].from]) x=min(x,e[i].w);
     for(int i=from[t];i;i=from[e[i].from]) {
         anss+=e[i].c*x;e[i].w-=x; e[i^1].w+=x;
     }
}
int main () {
     n=read(); int road=read();m=read(); t=n;
     for(i=1;i<=road;i++) {
          zs1=read(),zs2=read(),zs3=read(),zs4=read();
          ins(zs1,zs2,zs3,0); ins(zs2,zs1,0,0); num[++tot1][0]=tot-1,num[tot1][1]=zs3;
          E[++tot2].from=zs1; E[tot2].to=zs2; E[tot2].c=zs4;
     }
     while(bfs()){
          while(sum=find(1,inf)) flow+=sum;
     }
     printf("%d ",flow);
     for(i=1;i<=tot1;i++) {
          e[num[i][0]].w=num[i][1]; e[num[i][0]^1].w=0;
     }
     int u,v,cost;
     ins(0,1,flow+m,0),ins(1,0,0,0);
     for(i=1;i<=tot2;i++) {
          u=E[i].from,v=E[i].to,cost=E[i].c;
          ins(u,v,inf,cost); ins(v,u,0,-cost);
     }
     while(spfa()) run(); 
     printf("%d",anss); 
}