[caioj 1115] 网络流入门1 --- dinic最大流

时间:2022-11-05 04:30:16

【题目描述】

一个有向图, M条有向边和N个点,求点1 到点N的最大流。

【输入格式】
第1行:2个整数M (2 <= M <= 200) 和N (0 <= N <= 200)。(注意给出的格式M N)
下来M行: 每行有三个整数:x,y,c。表示一条从点x到点y的有向边,流量为c (0<= c <= 10,000,000)。
【输出格式】
输出一个整数,即最大流量。
SAMPLE INPUT
5 4
1 2 40
1 4 20
2 4 20
2 3 30
3 4 10
SAMPLE OUTPUT
50

提供一个理论上可以卡最短路的数据:

【样例输入】
4 4
1 2 10
2 4 5
2 3 20
3 4 10

【样例输出】

10

分析

dinic(LRJ蓝书模版)

代码

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <vector>
#include <queue>

#define IL inline
#define INF 0x7f7f7f7f

#define open(s) freopen(s".in","r",stdin); freopen(s".out","w",stdout);
#define close fclose(stdin); fclose(stdout);

using namespace std;

inline int read()
{
char c=getchar();
int sum=0,k=1;
for(;'0'>c || c>'9';c=getchar())
if(c=='-') k=-1;
for(;'0'<=c && c<='9';c=getchar()) sum=sum*10+c-'0';
return sum*k;
}

struct Edge
{
int from, to, cap, flow;
};

int n, m;
vector<Edge> edge;
vector<int> G[205];
bool vis[205];
int dep[205];
int cur[205];

IL int min_(int x, int y) { return x < y ? x : y; }

IL void add(int u, int v, int c, int k)
{
edge.push_back((Edge){u, v, c, 0});
edge.push_back((Edge){v, u, 0, 0});
G[u].push_back(k);
G[v].push_back(k|1);
}

IL bool Bfs()
{
memset(vis, 0, sizeof(vis));
queue<int>q;
dep[1] = 0;
vis[1] = 1;
q.push(1);
for(int t = 1, u, v; t;)
{
u= q.front(); q.pop(); --t;
for(int i = 0, s = G[u].size(); i < s; ++i)
{
Edge &e = edge[G[u][i]];
if(!vis[e.to] && e.cap > e.flow)
{
vis[e.to] = 1;
dep[e.to] = dep[u] + 1;
q.push(e.to); ++t;
}
}
}
return vis[n];
}

IL int Dfs(int u, int a)
{
if(u == n || !a) return a;
int flow = 0, f;
for(int &i = cur[u], s = G[u].size(); i < s; ++i)
{
Edge &e = edge[G[u][i]];
if(dep[u] + 1 == dep[e.to] && (f = Dfs(e.to, min_(a, e.cap - e.flow))) > 0)
{
e.flow += f;
edge[G[u][i] ^ 1].flow -= f;
flow += f;
a -= f;
if(!a) break;
}
}
return flow;
}

int main()
{
open("1115");

m = read(); n = read();

for(int i = 0, u, v, c; i < m; ++i)
{
u = read(); v = read(); c = read();
add(u, v, c, i << 1);
}

int flow = 0;
for(; Bfs();)
{
memset(cur, 0, sizeof(cur));
flow += Dfs(1, INF);
}

printf("%d\n",flow);


close;
return 0;
}