题目地址
题目大意
就是让你求出最小割应该割几条边
解题思路
利用dinic算法来计算。首先解题过程是给每个边的权值value变为value*(e+1)+1 (e为图的边数,乘的数字只要比边数大就行),然后跑一遍最大流,求最小割的边数就flow%(e+1)。
解释一下,先说加一是为啥吧!如果一个边是最小割的边,它一定是满流的,所以你容量加了1,它提供的流量肯定也加一,而其他的边肯定是不变的。所以最后原本的flow其实是变为了flow*(e+1)+p,p就是每个最小割边提供的那个1。那么对这个最大流模(e+1)就剩下了p,也就是最小割的边。
在说一下为什么乘(e+1)呢?首先看如果你对于每个边都直接加1,那最后的flow就变为了flow+p,这个p就没有办法求了,最大流也废了。那可能有人就想到了能不能跑两边,我第一遍先求原本的flow,然后给各个边加1再求flow+p然后相减。这个对错我也不太清楚,想尝试的可以试一试。我们有更简便的方法就是对边值乘上一个大于边数的数。这样的话,求出来的是flow*(e+1)+p,我们乘的数字是大于边数的数字,如果所有边都是最小割的边,那么p就等于m, 式子变为了flow*(e+1)+m如果e+1小于边数的话,在取模(e+1)就会影响到m的值。因此必须要乘以一个大于边数的数字才不会影响最小割的值。
AC代码
#include <iostream>
#include <queue>
#include <vector>
#include <cstring>
using namespace std;
const int INF = 0x3f3f3f3f;
const int MAX_M = 1055;
const int MAX_V = 300;
int level[MAX_V];
struct edge
{
int to, cap, rev;
};
vector<edge> G[MAX_V];
void add_edge(int from, int to, int cost)
{
G[from].push_back((edge){to, cost, G[to].size()});
G[to].push_back((edge){from, 0, G[from].size()-1});
}
void bfs(int s)
{
memset(level, -1, sizeof(level));
queue<int> que;
level[s] = 0;
que.push(s);
while (!que.empty())
{
int v = que.front();
que.pop();
for (int i=0; i<G[v].size(); i++)
{
edge &e = G[v][i];
if (e.cap && level[e.to] < 0)
{
level[e.to] = level[v]+1;
que.push(e.to);
}
}
}
}
int dfs(int v, int t, int f)
{
if (v == t)
{
return f;
}
for (int i=0; i<G[v].size(); i++)
{
edge &e = G[v][i];
if (e.cap && level[v]+1==level[e.to])
{
int d = dfs(e.to, t, min(f, e.cap));
if (d > 0)
{
e.cap-=d;
G[e.to][e.rev].cap+=d;
return d;
}
}
}
return 0;
}
int dinic(int s, int t)
{
int flow = 0;
while (1)
{
bfs(s);
if (level[t]<0)
{
return flow;
}
int f;
while ((f = dfs(s, t, INF)) > 0)
flow+=f;
}
}
int main()
{
int t;
scanf("%d", &t);
while (t--)
{
int n, m;
int s, t;
scanf("%d %d", &n, &m);
scanf("%d %d", &s, &t);
for (int i=0; i<=n; i++)
{
G[i].clear();
}
for (int i=0; i<m; i++)
{
int x, y, cost;
scanf("%d %d %d", &x, &y, &cost);
add_edge(x, y, cost*MAX_M+1);
}
printf("%d\n", dinic(s, t)%MAX_M);
}
return 0;
}