Nubulsa Expo. 全局最小割

最新推荐文章于 2019-11-30 17:34:03 发布

Samaritan_infi

最新推荐文章于 2019-11-30 17:34:03 发布

阅读量274

点赞数

CC 4.0 BY-SA版权

分类专栏：网络流

本文链接：https://blog.youkuaiyun.com/Samaritan_infi/article/details/78245121

网络流专栏收录该内容

9 篇文章

订阅专栏

本文介绍了一个基于C++实现的最小割算法，通过不断合并顶点来寻找网络中最小的割集，适用于解决最大流问题中的最小割求解。文章提供了完整的代码示例，并解释了如何初始化图以及添加边。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include <bits/stdc++.h>
//点标从0-n-1, 开始时先init 复杂度n^3
//对于边(u,v,flow):
//g[u][v]+=flow;
//g[v][u]+=flow;
typedef long long ll;
const int N = 305;
const ll inf = 1e18;
ll g[N][N], w[N];
int a[N], v[N], na[N];
void add_edge(int u , int v , int flow){
    g[u][v] += flow ;
    g[v][u] += flow ;
}
ll mincut(int n) {
	int i, j, pv, zj;
	ll best = inf;
	for(i = 0; i < n; i ++) v[i] = i;

	while(n > 1) {
		for(a[v[0]] = 1, i = 1; i < n; i ++) {
			a[v[i]] = 0;
			na[i-1] = i;
			w[i] = g[v[0]][v[i]];
		}
		for(pv = v[0], i = 1; i < n; i ++) {
			for(zj = -1, j = 1; j < n; j ++)
				if(!a[v[j]] && (zj < 0 || w[j] > w[zj])) zj = j;

			a[v[zj]] = 1;
			if(i == n-1) {
				if(best > w[zj]) best = w[zj];
				for(i = 0; i < n; i ++) {
					g[v[i]][pv] = g[pv][v[i]] += g[v[zj]][v[i]];
				}
				v[zj] = v[--n];
				break;
			}
			pv = v[zj];
			for(j = 1; j < n; j ++) if(!a[v[j]])
				w[j] += g[v[zj]][v[j]];
		}
	}
	return best;
}
void init(int n){
    for(int i = 0; i < n; i ++)
        for(int j = 0; j < n; j ++)
            g[i][j] = g[j][i] = 0;
}

int main()
{
    int n , m , st ;
    while( ~ scanf("%d %d %d" , & n , & m , & st ))
    {
        if(!n && !n && !st) break ;
        init(n) ;
        int a , b , c;
        for(int i = 0 ; i < m ; i ++ ){
            scanf("%d %d %d" , &a , &b , &c) ;
            add_edge(a - 1 , b - 1 , c) ;
        }
        printf("%lld\n" , mincut(n)) ;
    }
    return 0 ;
}