费用流 模板 洛谷3381

题目描述

如题,给出一个网络图,以及其源点和汇点,每条边已知其最大流量和单位流量费用,求出其网络最大流和在最大流情况下的最小费用。

输入输出格式

输入格式:
第一行包含四个正整数N、M、S、T,分别表示点的个数、有向边的个数、源点序号、汇点序号。

接下来M行每行包含四个正整数ui、vi、wi、fi,表示第i条有向边从ui出发,到达vi,边权为wi(即该边最大流量为wi),单位流量的费用为fi。

输出格式:
一行,包含两个整数,依次为最大流量和在最大流量情况下的最小费用。

输入输出样例

输入样例#1: 复制
4 5 4 3
4 2 30 2
4 3 20 3
2 3 20 1
2 1 30 9
1 3 40 5
输出样例#1: 复制
50 280
说明

时空限制:1000ms,128M

(BYX:最后两个点改成了1200ms)

数据规模:

对于30%的数据:N<=10,M<=10

对于70%的数据:N<=1000,M<=1000

对于100%的数据:N<=5000,M<=50000

样例说明:

这里写图片描述

如图,最优方案如下:

第一条流为4–>3,流量为20,费用为3*20=60。

第二条流为4–>2–>3,流量为20,费用为(2+1)*20=60。

第三条流为4–>2–>1–>3,流量为10,费用为(2+9+5)*10=160。

故最大流量为50,在此状况下最小费用为60+60+160=280。

故输出50 280。

Edmonds-Karp 再加 spfa跑最短路。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>

using namespace std;
const int MAXN = 5005;
const int MAXM = 50005;
const int inf = 0x3f3f3f3f;

struct Edge{
    int nxt,to,v,cost;
}edge[MAXM*2];

int n,m,cnt=1,ans,S,T,incf[MAXN],fa[MAXN];
int head[MAXN],dis[MAXN],maxflow,pre[MAXN];
bool vis[MAXN];

inline void init(){
    memset(vis,0,sizeof(vis));
    memset(dis,0x3f,sizeof(dis));
    memset(incf,0,sizeof(incf));
    dis[S]=0;incf[S]=inf;vis[S]=1;
}

inline void add(int bg,int ed,int val,int c){
    edge[++cnt].to=ed;
    edge[cnt].v=val;
    edge[cnt].cost=c;
    edge[cnt].nxt=head[bg];
    head[bg]=cnt;
}

inline bool spfa(){
    init();
    queue<int> q;
    q.push(S);
    while(q.size()){
        int x=q.front();q.pop();
        for(int i=head[x];i;i=edge[i].nxt){
            int u=edge[i].to;
                if(edge[i].v && dis[u]>dis[x]+edge[i].cost){
                dis[u]=dis[x]+edge[i].cost;
                incf[u]=min(incf[x],edge[i].v);
                pre[u]=i;
                fa[u]=x;
                if(!vis[u]) vis[u]=1,q.push(u);
            } 
        } 
        vis[x]=0;
    }
    if(dis[T]==inf) return false;
    return true;
}

inline void update(){
    int x=T;
    while(x!=S){
        int i=pre[x];
        edge[i].v-=incf[T];
        edge[i^1].v+=incf[T];
        x=fa[x];
    }
    maxflow+=incf[T];
    ans+=incf[T]*dis[T];
}

int main() {
    scanf("%d%d%d%d",&n,&m,&S,&T);
    for(register int i=1;i<=m;i++){
        int a,b,c,d;
        scanf("%d%d%d%d",&a,&b,&d,&c);
        add(a,b,d,c);
        add(b,a,0,-c);
    }
    while(spfa()) update();
    printf("%d %d",maxflow,ans);
    return 0;
}
posted @ 2018-05-01 21:54  Monster_Qi  阅读(103)  评论(0编辑  收藏  举报