P2604 [ZJOI2010]网络扩容

题目描述

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

输入格式

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

输出格式

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

---

题解：
第一问直接最大流即可，第二问设置超级源点连向1，容量为K，费用为0，表示限制最大容量拓展了K，然后利用残余网络，在每条边之间再建一条边，容量为INF，费用为co[i]。

因为要跑最小费用最大流，而原图是无花费的最大流，所以只需要再残余网络上继续增广，就可以保证抛出最小费用了

---

AC代码：

#pragma GCC optimize(2)
#include<bits/stdc++.h>
#include<ext/rope>
using namespace std;
using namespace __gnu_cxx;
#define LL long long
#define pii pair<int,int>
#define mp(a,b) make_pair(a,b)
const int MAXN = 2e5+50;
const int MAXM = 2e6+50;
const int MOD = 1e9+7;
const int INF = 0x3f3f3f3f;
int n,m,k,s,t,tot=1,head[MAXN],to[MAXM],w[MAXM],co[MAXM],nxt[MAXM],h[MAXN];
int a[MAXN],b[MAXN],c[MAXN],d[MAXN],vis[MAXN],dis[MAXN],flow[MAXN],pre[MAXN];
inline void ade(int u,int v,int ww,int cost){
    to[++tot]=v; w[tot]=ww; co[tot]=cost; nxt[tot]=head[u]; head[u]=tot;
}
inline void add(int u,int v,int w,int cost){ ade(u,v,w,cost); ade(v,u,0,-cost); }
inline int bfs(){
    queue<int> que; que.push(s); memset(h,0,sizeof(h)); h[s]=1;
    while(!que.empty()){
        int u=que.front(); que.pop();
        for(int i=head[u];i;i=nxt[i]){
            if(w[i] && !h[to[i]]){
                h[to[i]]=h[u]+1; que.push(to[i]);
            }
        }
    }
    return h[t];
}
inline int dfs(int x,int f){
    if(x==t) return f; int fl=0;
    for(int i=head[x];i;i=nxt[i]){
        if(w[i] && h[to[i]]==h[x]+1){
            int mi=dfs(to[i],min(f,w[i]));
            w[i]-=mi; w[i^1]+=mi; fl+=mi; f-=mi;
        }
    }
    if(!fl) h[x]=-1;
    return fl;
}
inline int dinic(){
    int res=0;
    while(bfs()) res+=dfs(s,INF);
    return res;
}
inline int spfa(){
    for(int i=0;i<=t;i++) dis[i]=INF;
    queue<int> que; que.push(s); dis[s]=0; vis[s]=1; flow[s]=INF;
    while(!que.empty()){
        int u=que.front(); que.pop(); vis[u]=0;
        for(int i=head[u];i;i=nxt[i]){
            if(w[i] && dis[to[i]]>dis[u]+co[i]){
                dis[to[i]]=dis[u]+co[i]; pre[to[i]]=i;
                flow[to[i]]=min(flow[u],w[i]);
                if(!vis[to[i]]) vis[to[i]]=1,que.push(to[i]);
            }
        }
    }
    return dis[t]!=INF;
}
inline int EK(){
    int res=0;
    while(spfa()){
        res+=flow[t]*dis[t]; int x=t;
        while(x!=s){
            int i=pre[x];
            w[i]-=flow[t]; w[i^1]+=flow[t]; x=to[i^1];
        }
    }
    return res;
}
signed main(){
#ifndef ONLINE_JUDGE
    freopen("C:\\Users\\Administrator\\Desktop\\in.txt","r",stdin);
#endif // ONLINE_JUDGE
    scanf("%d%d%d",&n,&m,&k); s=1; t=n;
    for(int i=1;i<=m;i++) scanf("%d%d%d%d",&a[i],&b[i],&c[i],&d[i]);
    for(int i=1;i<=m;i++) add(a[i],b[i],c[i],0);
    int res1=dinic(); s=0; add(s,1,k,0);
    for(int i=1;i<=m;i++) add(a[i],b[i],1e8,d[i]);
    int res2=EK();
    printf("%d %d\n",res1,res2);
    return 0;
}