364 网络流 费用流 EK 算法

视频链接：https://www.bilibili.com/video/BV1SU4y1D77p/

1. //<a href="https://www.luogu.com.cn/problem/P3381" target="_blank" rel="noopener"> Luogu P3381 【模板】最小费用最大流</a>
2. #include <iostream>
3. #include <cstring>
4. #include <algorithm>
5. #include <queue>
6. using namespace std;
8. const int N=5010,M=100010,INF=1e8;
9. int n,m,S,T;
10. struct edge{int v,c,w,ne;}e[M];
11. int h[N],idx=1;//从2,3开始配对
12. int d[N],mf[N],pre[N],vis[N];
13. int flow,cost;
15. void add(int a,int b,int c,int d){
16.   e[++idx]={b,c,d,h[a]};
17.   h[a]=idx;
18. }
19. bool spfa(){
20.   memset(d,0x3f,sizeof d);
21.   memset(mf,0,sizeof mf);
22.   queue<int> q; q.push(S);
23.   d[S]=0, mf[S]=INF, vis[S]=1;
24.   while(q.size()){
25.     int u=q.front(); q.pop();
26.     vis[u]=0;
27.     for(int i=h[u];i;i=e[i].ne){
28.       int v=e[i].v,c=e[i].c,w=e[i].w;
29.       if(d[v]>d[u]+w && c){
30.         d[v]=d[u]+w; //最短路
31.         pre[v]=i;
32.         mf[v]=min(mf[u],c);
33.         if(!vis[v]){
34.           q.push(v); vis[v]=1;
35.         }
36.       }
37.     }
38.   }
39.   return mf[T]>0;
40. }
41. void EK(){
42.   while(spfa()){
43.     for(int v=T;v!=S;){
44.       int i=pre[v];
45.       e[i].c-=mf[T];
46.       e[i^1].c+=mf[T];
47.       v=e[i^1].v;
48.     }
49.     flow+=mf[T]; //累加可行流
50.     cost+=mf[T]*d[T];//累加费用   
51.   }
52. }
53. int main(){
54.   scanf("%d%d%d%d",&n,&m,&S,&T);
55.   int a,b,c,d;
56.   while(m --){
57.     scanf("%d%d%d%d",&a,&b,&c,&d);
58.     add(a,b,c,d);
59.     add(b,a,0,-d);
60.   }
61.   EK();
62.   printf("%d %d\n",flow,cost);
63.   return 0;
64. }

复制代码

 

免责声明：如果侵犯了您的权益，请联系站长，我们会及时删除侵权内容，谢谢合作！