当前位置: 首页 > news >正文

网络流 最小割 Dinic算法

image

标准模板

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=210,M=5e3+10;
int n,m,s,t,d[N],cur[N],vis[N];
int h[N],e[M<<1],ne[M<<1],id=1;//从2,3开始配对
LL c[M<<1];
void add(int u,int v,LL w){++id;e[id]=v;ne[id]=h[u];c[id]=w;h[u]=id;
}LL dfs(int u,LL mf){//多路增广if(u==t)return mf;LL sum=0;for(int i=cur[u];i;i=ne[i]){cur[u]=i;//当前弧优化int v=e[i];if(d[v]==d[u]+1&&c[i]){LL f=dfs(v,min(mf,c[i]));sum+=f;//累加u的流出流量c[i]-=f;c[i^1]+=f;//更新残留网mf-=f;//减少u的剩余流量if(mf==0)//余量优化break;}}if(sum==0)//残枝优化d[u]=0;return sum;
}bool bfs(){//对点分层,找增广路memset(d,0,sizeof(d));queue<int> q;q.push(s);d[s]=1;while(q.size()){int u=q.front();q.pop();for(int i=h[u];i;i=ne[i]){int v=e[i];if(d[v]==0&&c[i]){d[v]=d[u]+1;q.push(v);if(v==t)return true;}}}return false;
}LL dinic(){//累加可行流LL flow=0;while(bfs()){memcpy(cur,h,sizeof(h));flow+=dfs(s,1e9);}return flow;
}void mincut(int u){vis[u]=1;for(int i=h[u];i;i=ne[i]){int v=e[i];if(!vis[v]&&c[i])mincut(v);}
}void reset_edges(){//每条边是成对存储,正边id为偶数,反边id为奇数for(int i=2;i<=id;i++){if(i&1){c[i]=0;}else{if(c[i])c[i]=1e9;elsec[i]=1;}}
}
int main(){cin.tie(nullptr)->sync_with_stdio(false);cin>>n>>m>>s>>t;for(int i=1;i<=m;i++){LL u,v,w;cin>>u>>v>>w;add(u,v,w);add(v,u,0);}//最小割cout<<dinic()<<endl;//最小割划分mincut(s);//访问过的点为S集合for(int i=1;i<=n;i++){if(vis[i])cout<<i<<' ';}cout<<endl;//没访问过的点为T集合for(int i=1;i<=n;i++){if(!vis[i])cout<<i<<' ';}cout<<endl;//求最小割的最少边数//重建边时应当把第一遍dinic中剩余容量为0的正向边的边权设为1,//其他正向边设为无穷大,反向边都设为零,//因为只有流满的边才是最小割中的边。reset_edges();cout<<dinic()<<endl;return 0;
}

模板题:洛谷p1344

code

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=210,M=5e3+10;
int n,m,s,t,d[N],cur[N],vis[N];
int h[N],e[M<<1],ne[M<<1],id=1;//从2,3开始配对
LL c[M<<1];
void add(int u,int v,LL w){++id;e[id]=v;ne[id]=h[u];c[id]=w;h[u]=id;
}LL dfs(int u,LL mf){//多路增广if(u==t)return mf;LL sum=0;for(int i=cur[u];i;i=ne[i]){cur[u]=i;//当前弧优化int v=e[i];if(d[v]==d[u]+1&&c[i]){LL f=dfs(v,min(mf,c[i]));sum+=f;//累加u的流出流量c[i]-=f;c[i^1]+=f;//更新残留网mf-=f;//减少u的剩余流量if(mf==0)//余量优化break;}}if(sum==0)//残枝优化d[u]=0;return sum;
}bool bfs(){//对点分层,找增广路memset(d,0,sizeof(d));queue<int> q;q.push(s);d[s]=1;while(q.size()){int u=q.front();q.pop();for(int i=h[u];i;i=ne[i]){int v=e[i];if(d[v]==0&&c[i]){d[v]=d[u]+1;q.push(v);if(v==t)return true;}}}return false;
}LL dinic(){//累加可行流LL flow=0;while(bfs()){memcpy(cur,h,sizeof(h));flow+=dfs(s,1e9);}return flow;
}void mincut(int u){vis[u]=1;for(int i=h[u];i;i=ne[i]){int v=e[i];if(!vis[v]&&c[i])mincut(v);}
}void reset_edges(){//每条边是成对存储,正边id为偶数,反边id为奇数for(int i=2;i<=id;i++){if(i&1){c[i]=0;}else{if(c[i])c[i]=1e9;elsec[i]=1;}}
}
int main(){cin.tie(nullptr)->sync_with_stdio(false);// cin>>n>>m>>s>>t;cin>>n>>m;s=1;t=n;for(int i=1;i<=m;i++){LL u,v,w;cin>>u>>v>>w;add(u,v,w);add(v,u,0);}//最小割// cout<<dinic()<<endl;cout<<dinic()<<' ';//最小割划分// mincut(s);//访问过的点为S集合// for(int i=1;i<=n;i++){//     if(vis[i])//         cout<<i<<' ';// }// cout<<endl;//没访问过的点为T集合// for(int i=1;i<=n;i++){//     if(!vis[i])//         cout<<i<<' ';// }// cout<<endl;//求最小割的最少边数//重建边时应当把第一遍dinic中剩余容量为0的正向边的边权设为1,//其他正向边设为无穷大,反向边都设为零,//因为只有流满的边才是最小割中的边。reset_edges();cout<<dinic()<<endl;return 0;
}
http://www.rkmt.cn/news/15018.html

相关文章:

  • 15.VLANIF(2025年9月30日) - 教程
  • Pdfminer-Vulnerability-Research
  • 10.2笔记
  • Spring Boot 内置日志框架 Logback - 以及 lombok 介绍 - 教程
  • CF VP 记录
  • 原来你是这样的claude code aciton:没想象中好
  • FlareOn1 -- 5get_it
  • python语言手势控制音乐播放器代码QZQ
  • 2025 编码器厂家 TOP 企业品牌推荐排行榜,无磁,光学,脉冲,绝对型,伺服,机械多圈,工业,二进制,拉线编码器公司推荐
  • Spark专题-第三部分:性能监控与实战优化(1)-认识spark ui - 指南
  • 2025 年等离子清洗机厂家 TOP 企业品牌推荐排行榜,大气,真空,宽幅,微波,自动化,常压,低温,大腔体,射频,DBD,介质阻挡放电等离子清洗机公司推荐!
  • 完整教程:如何优雅的布局,height: 100% 的使用和 flex-grow: 1 的 min-height 陷阱
  • 2025担保合同律师事务所推荐,专业团队高效解决法律难题!
  • 2025年筒袋磁力泵实力厂家推荐榜:高效耐用与创新技术深度解
  • Android项目实现自动获取手机号一键登录功能
  • Qt编程: 正则表达式分析 - 实践
  • Manim实现渐变填充特效
  • Spring Boot 集成 Redis 全方位详解 - 指南
  • 十月牛气冲天计数题没做
  • datadome 隐私模式 ck设置
  • CPU温度查看(Core Temp)
  • 深入解析:python学智能算法(三十九)|使用PyTorch模块的normal()函数绘制正态分布函数图
  • 2025污水处理设备厂家 TOP 企业品牌推荐排行榜,一体化,生活,工业,养殖,医疗,农村,学校,餐厨,隧洞,高速污水处理设备公司推荐!
  • 详细介绍:告别“下次注意”,用这套结构化事故复盘方案就对了
  • 关于树状数组的一些东西
  • [问题记录] vmagent 增加 aggregation 表达式后,CPU 上升 2.43 倍, 内存上升 3.82 倍
  • CF1081F Tricky Interactor
  • JAVA SE 基础语法 —— A / 初识 - 指南
  • 2025机械加工供货厂家权威口碑排行:实力与服务深度解析!
  • 2025七水硫酸锌厂家权威推荐榜:优质供应与专业定制首选