最小费用最大流spfa + ek 算法模板

本文介绍了一个基于最小成本最大流算法的实现,通过SPFA算法寻找增广路径,并更新流量及费用,直至无法找到新的增广路径为止。适用于解决带有成本约束的网络流问题。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define rep(i, a, n) for (int i = a; i <= n; i++)
#define per(i, a, n) for (int i = n; i >= a; i--)
#define lowbit(x) ((x) & -(x))
#define lson l, mid, rt << 1
#define rson mid + 1, r, rt << 1 | 1
#define mem(a, b) memset(a, b, sizeof(a))
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
const double PI = acos(-1.0);
const int MAXN = 10000, MAXM = 100000, inf = 0x3f3f3f3f;
struct edge{
    int v,next,cap,flow,cost;
}e[MAXM];
int head[MAXN],tot,pre[MAXN],dis[MAXN];
bool vis[MAXN];
int N;//节点总个数,编号从0 - n
void init(int n){
    N = n;
    tot = 1;
    mem(head,-1);
}

void add(int u,int v,int cap,int cost){
    e[++tot].v = v, e[tot].cap = cap, e[tot].cost = cost, e[tot].flow = 0, e[tot].next = head[u], head[u] = tot;
    e[++tot].v = u, e[tot].cap = 0, e[tot].cost = -cost, e[tot].flow = 0, e[tot].next = head[v], head[v] = tot;
}

bool spfa(int s,int t){
    queue<int>q;
    rep(i,0,N){//初始化
        dis[i] = inf;
        vis[i] = 0;
        pre[i] = -1;
    }
    dis[s] = 0;
    vis[s] = 1;
    q.push(s);
    while(!q.empty()){
        int u = q.front();
        q.pop();
        vis[u] = 0;
        for(int i = head[u]; ~i ; i = e[i].next){
            int v = e[i].v;
            if(e[i].cap > e[i].flow && dis[v] > dis[u] + e[i].cost){
                dis[v] = dis[u] + e[i].cost;
                pre[v] = i;
                if(!vis[v]){
                    vis[v] = 1;
                    q.push(v);
                } 
            }
        }
    }
    if(pre[t] == -1) return 0;
    else return 1;
}

int minCostMaxflow(int s,int t,int &cost){
    int flow = 0;
    cost = 0;    
    while(spfa(s,t)){
        int Min = inf;
        for(int i = pre[t]; ~i; i = pre[e[i^1].v]){
            if(Min > e[i].cap - e[i].flow)
                Min = e[i].cap - e[i].flow;
        }
        for(int i = pre[t]; ~i; i = pre[e[i^1].v]){
            e[i].flow += Min;
            e[i^1].flow -= Min;
            cost += e[i].cost * Min;
        }
        flow += Min;
    }
    return flow;
}
int main()
{
    IOS
    int n,m,s,t;
    cin >> n >> m >> s >> t;
    init(n);
    rep(i,1,m){
        int u,v,cap,cost;
        cin >> u >> v >> cap >> cost;
        add(u,v,cap,cost);
    }   
    int maxflow , cost;
    maxflow = minCostMaxflow(s,t,cost);
    cout << maxflow << " " << cost << "\n";
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值