Template -「网络流 & 二分图」
EK.
很少用到，知道思想即可。
懒得写封装的屑。
queue<int> q;
int Cap[MAXN][MAXN], Flow[MAXN][MAXN], Aug[MAXN], fa[MAXN], n;
void Add_Cap(int u, int v, int C) { Cap[u][v] += C; }
int bfs(int s, int t) {
for(int i = 1; i <= n; i++) {
Aug[i] = 0;
fa[i] = 0;
}
while(!q.empty())
q.pop();
Aug[s] = INF;
q.push(s);
while(!q.empty()) {
int u = q.front();
q.pop();
for(int v = 1; v <= n; v++)
if(Cap[u][v] - Flow[u][v] > 0 && !Aug[v]) {
Aug[v] = Min(Aug[u], Cap[u][v] - Flow[u][v]);
fa[v] = u;
if(v == t)
return Aug[v];
q.push(v);
}
}
return 0;
}
int EK(int s, int t) {
int Delta = bfs(s, t), res = 0, u;
while(Delta) {
res += Delta;
for(int v = t; v != s; v = u) {
u = fa[v];
Flow[u][v] += Delta;
Flow[v][u] -= Delta;
}
Delta = bfs(s, t);
}
return res;
}
Dinic.
时间复杂度玄学的常用最大流板子。
struct Maximum_Flow {
#define Type int
struct edge {
int v, nxt;
edge() {}
edge(int V, int Nxt) {
v = V, nxt = Nxt;
}
} e[MAXM << 1];
int n, cnt, s, t;
Type Cap[MAXM << 1], Flow[MAXM << 1];
int Lab[MAXL], Cur[MAXL], head[MAXL], Col[MAXL];
queue<int> q;
void init(int N, int S, int T) {
for(int i = 0; i <= cnt; i++)
Flow[i] = 0, Cap[i] = 0;
n = N, cnt = 0, s = S, t = T;
for(int i = 1; i <= n; i++)
head[i] = -1, Col[i] = 0;
}
void Add_Edge(int u, int v, Type w) {
Cap[cnt] += w;
e[cnt] = edge(v, head[u]);
head[u] = cnt++;
e[cnt] = edge(u, head[v]);
head[v] = cnt++;
}
bool Lab_Vertex() {
for(int i = 1; i <= n; i++)
Lab[i] = 0;
Lab[t] = 1;
while(!q.empty())
q.pop();
q.push(t);
while(!q.empty()) {
int v = q.front();
q.pop();
for(int i = head[v], u; ~i; i = e[i].nxt) {
u = e[i].v;
if(!Lab[u] && Cap[i ^ 1] - Flow[i ^ 1]) {
Lab[u] = Lab[v] + 1;
q.push(u);
if(u == s)
return Lab[s];
}
}
}
return Lab[s];
}
Type Widen(int u, Type Limit) {
if(u == t)
return Limit;
Type Used = 0, Delta;
for(int i = Cur[u], v; ~i; i = e[i].nxt) {
v = e[i].v;
Cur[u] = i;
if(Lab[v] + 1 != Lab[u] || Cap[i] - Flow[i] <= 0)
continue;
Delta = Widen(v, Min(Limit - Used, Cap[i] - Flow[i]));
Used += Delta, Flow[i] += Delta, Flow[i ^ 1] -= Delta;
if(Used == Limit)
return Used;
}
return Used;
}
Type Dinic() {
Type res = 0;
while(Lab_Vertex()) {
for(int i = 1; i <= n; i++)
Cur[i] = head[i];
res += Widen(s, INF);
if(res < 0)
return INF;
}
return res;
}
void Color(int u) {
Col[u] = 1;
for(int i = head[u]; ~i; i = e[i].nxt)
if(Cap[i] - Flow[i] > 0 && !Col[e[i].v])
Color(e[i].v);
}
#undef Type
} Flow_Graph;
spfa + Dinic.
普通最小费用最大流。（小心被卡？
struct Minimum_Cost_Maximum_Flow {
#define Type int
struct edge {
int v, nxt;
Type Wei, Cap, Flow;
edge() {}
edge(int V, int Nxt, Type C, Type W, Type F) {
v = V, nxt = Nxt, Cap = C, Wei = W, Flow = F;
}
} e[MAXM << 1];
int head[MAXM << 1], cnt = 0;
void Add_Edge(int u, int v, Type c, Type w) {
e[cnt] = edge(v, head[u], c, w, 0);
head[u] = cnt++;
e[cnt] = edge(u, head[v], 0, -w, 0);
head[v] = cnt++;
}
int Cur[MAXN], n, s, t;
bool vis[MAXN], Flag[MAXN];
Type Dist[MAXN], Aug[MAXN], Flow, Cost;
void init(int N, int S, int T) {
n = N, s = S, t = T, cnt = 0;
for(int i = 1; i <= n; i++)
head[i] = -1;
}
bool spfa() {
for(int i = 1; i <= n; i++)
Dist[i] = INF, vis[i] = false, Aug[i] = INF;
Dist[s] = 0, vis[s] = true, Aug[s] = INF;
queue<int> q;
q.push(s);
while(!q.empty()) {
int u = q.front(); q.pop();
vis[u] = false;
for(int i = head[u], v; ~i; i = e[i].nxt) {
v = e[i].v;
if(e[i].Cap - e[i].Flow > 0 && Dist[v] > Dist[u] + e[i].Wei) {
Dist[v] = Dist[u] + e[i].Wei;
Aug[v] = Min(Aug[u], e[i].Cap - e[i].Flow);
if(!vis[v])
vis[v] = true, q.push(v);
}
}
}
return Dist[t] != INF;
}
Type Widen(int u, Type Limit) {
if(u == t)
return Limit;
Flag[u] = true;
Type Used = 0, Delta;
for(int i = Cur[u], v; ~i; i = e[i].nxt) {
v = e[i].v;
if(Flag[v] || Dist[u] + e[i].Wei != Dist[v] || e[i].Cap - e[i].Flow <= 0)
continue;
Cur[u] = i;
Delta = Widen(v, Min(Limit - Used, e[i].Cap - e[i].Flow));
Used += Delta, e[i].Flow += Delta, e[i ^ 1].Flow -= Delta;
if(Used == Limit)
break;
}
Flag[u] = false;
return Used;
}
void Dinic() {
Flow = 0, Cost = 0;
Type Delta;
while(spfa()) {
for(int i = 1; i <= n; i++)
Cur[i] = head[i];
Delta = Widen(s, INF);
Flow += Delta, Cost += Delta * Dist[t];
}
}
#undef Type
} Flow_Graph;
spfa + Dijkstra + EK.
时间复杂度稳定，不会被卡。呃不会被卡吧。
struct Minimum_Cost_Maximum_Flow {
#define Type int
struct edge {
int v, nxt;
Type Wei, Cap, Flow;
edge() {}
edge(int V, int Nxt, Type C, Type W, Type F) {
v = V, nxt = Nxt, Cap = C, Wei = W, Flow = F;
}
} e[MAXM << 1];
int head[MAXM << 1], cnt = 0;
void Add_Edge(int u, int v, Type c, Type w) {
e[cnt] = edge(v, head[u], c, w, 0);
head[u] = cnt++;
e[cnt] = edge(u, head[v], 0, -w, 0);
head[v] = cnt++;
}
struct node {
int x;
Type dis;
node() {}
node(int X, Type Dis) {
x = X, dis = Dis;
}
friend bool operator < (node One, node TheOther) {
return One.dis > TheOther.dis;
}
};
struct Back {
int Pre, id;
Back() {}
Back(int P, int Id) {
Pre = P, id = Id;
}
} Last[MAXN];
int Cur[MAXN], n, s, t;
bool vis[MAXN], Flag[MAXN];
Type Dist[MAXN], Aug[MAXN], h[MAXN], Flow, Cost;
void init(int N, int S, int T) {
n = N, s = S, t = T, cnt = 0;
for(int i = 1; i <= n; i++)
head[i] = -1;
}
void spfa(int s, int t) {
for(int i = 1; i <= n; i++)
h[i] = INF, vis[i] = false;
h[s] = 0, vis[s] = true;
queue<int> q;
q.push(s);
while(!q.empty()) {
int u = q.front(); q.pop();
vis[u] = false;
for(int i = head[u], v; ~i; i = e[i].nxt) {
v = e[i].v;
if(e[i].Cap - e[i].Flow > 0 && h[v] > h[u] + e[i].Wei) {
h[v] = h[u] + e[i].Wei;
if(!vis[v])
vis[v] = true, q.push(v);
}
}
}
}
bool Dijkstra(int s, int t) {
for(int i = 1; i <= n; i++)
Dist[i] = INF, Last[i] = Back(-1, -1), vis[i] = false, Aug[i] = INF;
priority_queue<node> q;
Dist[s] = 0;
q.push(node(s, Dist[s]));
while(!q.empty()) {
int u = q.top().x; q.pop();
if(vis[u])
continue;
vis[u] = true;
for(int i = head[u], v; ~i; i = e[i].nxt) {
v = e[i].v;
if(e[i].Cap - e[i].Flow > 0 && Dist[v] > Dist[u] + e[i].Wei + h[u] - h[v]) {
Last[v] = Back(u, i);
Dist[v] = Dist[u] + e[i].Wei + h[u] - h[v];
Aug[v] = Min(Aug[u], e[i].Cap - e[i].Flow);
q.push(node(v, Dist[v]));
}
}
}
return Dist[t] != INF;
}
void EK() {
Flow = 0, Cost = 0;
spfa(s, t);
while(Dijkstra(s, t)) {
Flow += Aug[t];
Cost += Aug[t] * (Dist[t] + h[t]);
int pos = t;
while(pos != s) {
e[Last[pos].id].Flow += Aug[t];
e[Last[pos].id ^ 1].Flow -= Aug[t];
pos = Last[pos].Pre;
}
for(int i = 1; i <= n; i++)
h[i] += Dist[i];
}
}
#undef Type
} Flow_Graph;
Kuhn-Munkres.
最大匹配板子。当然你也可以直接 Dinic。
struct Bipartite_Graph {
struct edge {
int v, nxt;
edge() {}
edge(int V, int Nxt) {
v = V, nxt = Nxt;
}
} e[MAXM << 1];
int head[MAXN], n, m, cnt;
void Add_Edge(int u, int v) {
e[cnt] = edge(v, head[u]);
head[u] = cnt++;
}
int Mat[MAXN][2], Tim[MAXN], tot;
void init(int N, int M) {
n = N, m = M;
for(int i = 1; i <= n; i++)
head[i] = -1, Tim[i] = 0, Mat[i][0] = 0, Mat[i][1] = 0;
cnt = 0, tot = 0;
}
bool dfs(int u) {
if (Tim[u] == tot)
return false;
Tim[u] = tot;
for (int i = head[u], v; ~i; i = e[i].nxt) {
v = e[i].v;
if (!Mat[v][1] || dfs(Mat[v][1])) {
Mat[v][1] = u, Mat[u][0] = v;
return true;
}
}
return false;
}
int calc() {
int ans = 0;
for (int i = 1; i <= m; i++)
Mat[i][1] = 0;
for (int i = 1; i <= n; i++)
Tim[i] = 0, Mat[i][0] = 0;
for (int i = n; i >= 1; i--) {
tot++;
ans += dfs(i);
}
return ans;
}
} Graph;
KM.
最小权最大匹配。
int n, m;
bool Vis[MAXN];
int Mat[MAXN], fa[MAXN];
LL w[MAXN][MAXN], Slack[MAXN], Val[MAXN][2];
void bfs(int S) {
for(int i = 0; i <= n; i++)
Slack[i] = INF, Vis[i] = false;
LL d;
int pos = 0, p, u;
for(Mat[pos] = S; Mat[pos]; pos = p) {
Vis[pos] = true, u = Mat[pos], d = INF;
for(int v = 1; v <= n; v++) {
if(Vis[v])
continue;
if(Val[u][0] + Val[v][1] - w[u][v] < Slack[v]) {
Slack[v] = Val[u][0] + Val[v][1] - w[u][v];
fa[v] = pos;
}
if(Slack[v] < d) {
d = Slack[v];
p = v;
}
}
for(int v = 0; v <= n; v++)
if(Vis[v])
Val[Mat[v]][0] -= d, Val[v][1] += d;
else
Slack[v] -= d;
}
for(; pos; pos = fa[pos])
Mat[pos] = Mat[fa[pos]];
}
LL KM() {
LL res = 0;
for(int i = 1; i <= n; i++)
bfs(i);
for(int i = 1; i <= n; i++)
res += w[Mat[i]][i];
return res;
}
Template -「网络流 & 二分图」的更多相关文章
- Template -「整体二分」
写的简单.主要是留给自己做复习资料. 「BZOJ1901」Dynamic Rankings. 给定一个含有 \(n\) 个数的序列 \(a_1,a_2 \dots a_n\),需要支持两种操作: Q ...
- Template -「矩阵 - 行列式」
#include <cstdio> int Abs(int x) { return x < 0 ? -x : x; } int Max(int x, int y) { return ...
- [loj #6003]「网络流 24 题」魔术球 二分图最小路径覆盖，网络流
#6003. 「网络流 24 题」魔术球 内存限制:256 MiB时间限制:1000 ms标准输入输出 题目类型:传统评测方式:Special Judge 上传者: 匿名 提交提交记录统计讨论测试数据 ...
- LOJ6000 - 「网络流 24 题」搭配飞行员
原题链接 题意简述 求二分图的最大匹配. 题解 这里写的是匈牙利算法. 表示节点的当前匹配. 为真表示在这一轮匹配中,无法给节点一个新的匹配.所以如果为真就不用再dfs它了,直接continue就好. ...
- LibreOJ #6013. 「网络流 24 题」负载平衡 最小费用最大流 供应平衡问题
#6013. 「网络流 24 题」负载平衡 内存限制:256 MiB时间限制:1000 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: 匿名 提交提交记录统计讨论测试数据 题目描述 ...
- LibreOJ #6008. 「网络流 24 题」餐巾计划 最小费用最大流 建图
#6008. 「网络流 24 题」餐巾计划 内存限制:256 MiB时间限制:1000 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: 匿名 提交提交记录统计讨论测试数据 题目描述 ...
- LibreOJ #6007. 「网络流 24 题」方格取数 最小割 最大点权独立集 最大流
#6007. 「网络流 24 题」方格取数 内存限制:256 MiB时间限制:1000 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: 匿名 提交提交记录统计讨论测试数据 题目描述 ...
- 【刷题】LOJ 6227 「网络流 24 题」最长k可重线段集问题
题目描述 给定平面 \(\text{xoy}\) 上 \(n\) 个开线段组成的集合 \(\text{I}\) ,和一个正整数 \(k\) ,试设计一个算法. 从开线段集合 \(\text{I}\) ...
- 「THUWC 2017」随机二分图
「THUWC 2017」随机二分图 解题思路 : 首先有一个 \(40pts\) 的做法: 前 \(20pts\) 暴力枚举最终的匹配是怎样的,check一下计算方案数,后 \(20pts\) 令 \ ...
随机推荐
- Java 18 新增@snipppet标签，注释中写样例代码更舒适了！
在这次的Java 18中,新增了一个@snipppet标签,主要用于JavaDoc中需要放示例代码的场景.其实在Java 18之前,已经有一个@code标签,可以用于在JavaDoc中编写小段的代码内 ...
- ubuntu 16.04，ros kinetic 使用husy_gazebo
我当前使用的是ubuntu 16.04,ros kinetic ,Gazebo版本为7.0.protoc需要确保版本为2.6.1,而我当前的为3.4.0,因此需要将系统中的protoc替换为2.6.1 ...
- 使用Spring MVC开发RESTful API(续)
使用多线程提高REST服务性能 异步处理REST服务,提高服务器吞吐量 使用Runnable异步处理Rest服务 AsyncController.java @RestController @GetMa ...
- 535. Encode and Decode TinyURL - LeetCode
Question 535. Encode and Decode TinyURL Solution 题目大意:实现长链接加密成短链接,短链接解密成长链接 思路:加密成短链接+key,将长链接按key保存 ...
- Spring Boot整合模板引擎jsp
jsp也算是一种模板引擎吧.整合jsp前,先说一下运行SpringBoot项目的几种方式 1. 运行SpringBoot项目的几种方式 1.1 使用内嵌Tomcat运行项目 在IDE中右键运行启动类, ...
- Spring Boot中的微信支付（小程序）
前言 微信支付是企业级项目中经常使用到的功能,作为后端开发人员,完整地掌握该技术是十分有必要的. logo 一.申请流程和步骤 图1-1 注册微信支付账号 获取微信小程序APPID 获取微信商家的商户 ...
- nazo.io 通关记录
游戏网址 说在前面 答案错误页面 nazo.io/wrong 攻略 第0关 谜.io 纯粹是欢迎你来游戏. 所以他给你的start就是答案. 第1关 欢迎 它用灰体字写了key: welcome 直接 ...
- MTK 平台sensor arch 介绍-scp
架构介绍 路径:vendor/mediatek/proprietary/tinysys/scp 1.[build]编译相关 2.[driver]scp 的driver,I2C,power,eint 3 ...
- unittest框架里的常用断言方法：用于检查数据
1.unittest框架里的常用断言方法:用于检查数据. (1)assertEqual(x,y) 检查两个参数类型相同并且值相等.(2)assertTrue(x) 检查唯一的参数值等于True(3)a ...
- Tensor的组合与分块
>>> a = torch.Tensor([[1,2],[3,4]])>>> atensor([[1., 2.], [3., 4.]]) >>> ...