转自      http://www.cppblog.com/imky/articles/123427.html
  
一个无向连通网络,去掉一个边集可以使其变成两个连通分量则这个边集就是割集;最小割集当然就权和最小的割集。

可以用最小切割最大流定理:

1.min=MAXINT,确定一个源点

2.枚举汇点

3.计算最大流,并确定当前源汇的最小割集,若比min小更新min

4.转到2直到枚举完毕

5.min即为所求输出min

不难看出复杂度很高:枚举汇点要O(n),最短增广路最大流算法求最大流是O((n^2)m)复杂度,在复杂网络中O(m)=O(n^2),算法总复杂度就是O(n^5);哪怕采用最高标号预进流算法求最大流O((n^2)(m^0.5)),算法总复杂度也要O(n^4)

所以用网络流算法求解最小割集复杂度不会低于O(n^4)。

---------

prim算法不仅仅可以求最小生成树,也可以求“最大生成树”。最小割集Stoer-Wagner算法就是典型的应用实例。

求解最小割集普遍采用Stoer-Wagner算法,不提供此算法证明和代码,只提供算法思路:

1.min=MAXINT,固定一个顶点P

2.从点P用“类似”prim的s算法扩展出“最大生成树”,记录最后扩展的顶点和最后扩展的边

3.计算最后扩展到的顶点的切割值(即与此顶点相连的所有边权和),若比min小更新min

4.合并最后扩展的那条边的两个端点为一个顶点(当然他们的边也要合并,这个好理解吧?)

5.转到2,合并N-1次后结束

6.min即为所求,输出min

prim本身复杂度是O(n^2),合并n-1次,算法复杂度即为O(n^3)

如果在prim中加堆优化,复杂度会降为O((n^2)logn)

这个Stoer-Wagner算法可以参见这篇paper(http://docs.google.com/fileview?id=0BwxLvD9mcDNtMjk3MWVkMTAtZjMzNi00ZWE3LTkxYjQtYTQwNzcyZTk3Njk2&hl=en), 其核心思想是迭代缩小规模, 算法基于这样一个事实:

对于图中任意两点s和t, 它们要么属于最小割的两个不同集中, 要么属于同一个集.

如果是后者, 那么合并s和t后并不影响最小割. 基于这么个思想, 如果每次能求出图中某两点之间的最小割, 然后更新答案后合并它们再继续求最小割, 就得到最终答案了. 算法步骤如下:

1. 设最小割cut=INF, 任选一个点s到集合A中, 定义W(A, p)为A中的所有点到A外一点p的权总和.

2. 对刚才选定的s, 更新W(A,p)(该值递增).

3. 选出A外一点p, 且W(A,p)最大的作为新的s, 若A!=G(V), 则继续2.

4. 把最后进入A的两点记为s和t, 用W(A,t)更新cut.

5. 新建顶点u, 边权w(u, v)=w(s, v)+w(t, v), 删除顶点s和t, 以及与它们相连的边.

6. 若|V|!=1则继续1.

看起来很简单, 每次像做最大生成树一样选最大"边"(注意, 这里其实不是边, 而是已经累计的权值之和, 就当是加权的度好了), 然后把最后进入的两个点缩到一块就可以了. 合并点最多有n-1次, 而不加堆优化的prim是O(n^2)的, 所以最终复杂度O(n^3), 要是你有心情敲一大坨代码, 还可以在稀疏图上用Fibonacci Heap优化一下, 不过网上转了一圈, 大多都是说能用Fibonacci Heap优化到怎样怎样的复杂度, 真正能自己写出来的恐怕也没几个, 看看uoregon(俄勒冈大学)的一大坨代码就有点寒. (http://resnet.uoregon.edu/~gurney_j/jmpc/fib.html)

特别注意几个地方, 网上的好几个Stoer-Wagner版本都存在一些小错误:

1. 算法在做"最大生成树"时更新的不是普通意义上的最大边, 而是与之相连的边的权值和, 当所有边都是单位权值时就是累计度.

2. "最后进入A的两点记为s和t", 网上对s有两种解释, 一是在t之前一个加进去的点, 二是t的前趋节点, 也就是最后选择的那条边的另一端. 正解是第一种!

3. 对于稠密图, 比如这题, 我用堆, 映射二分堆, 或者STL的优先队列都会TLE, 还不如老老实实O(n^3).


 1#include <iostream>
 2#include <cmath>
 3using namespace std;
 4int mat[600][600];
 5int res;
 6//Stoer-Wagner算法,加了自己看得懂的备注
 7//无向图全局最小割,用求prim类似方法o(n^3),学习了一个下午……
 8//一开始用枚举源点汇点的最大流求解,复杂度o(n^5) 超时
 9int min(int a,int b)
10{
11    if(a<b) return a;
12    else return b;
13}
14
15void Mincut(int n) 
16{
17    int node[600], dist[600];
18    bool visit[600];
19    int i, prev, j, k;
20    for (i = 0; i < n; i++)
21        node[i] = i;
22    while (n > 1) 
23    {
24        int maxj = 1;
25        for (i = 1; i <n; i++) //初始化到已圈集合的割大小
26    { 
27            dist[node[i]] = mat[node[0]][node[i]];
28            if (dist[node[i]] > dist[node[maxj]])
29                maxj = i;
30        }
31        prev = 0;
32        memset(visit, false, sizeof (visit));
33        visit[node[0]] = true;
34        for (i = 1; i < n; i++)
35    {
36                        if (i == n - 1)//只剩最后一个没加入集合的点,更新最小割
37            { 
38                                res = min(res, dist[node[maxj]]);
39                                for (k = 0; k < n; k++) //合并最后一个点以及推出它的集合中的点
40                {
41                    mat[node[k]][node[prev]]=(mat[node[prev]][node[k]]+=mat[node[k]][node[maxj]]);
42                }
43                                node[maxj] = node[--n]; //缩点后的图
44                        }
45                        visit[node[maxj]] = true;
46                        prev = maxj;
47                        maxj = -1;
48                        for (j = 1; j < n; j++)
49            {
50                                if (!visit[node[j]]) //将上次求的maxj加入集合,合并与它相邻的边到割集
51                { 
52                                      dist[node[j]] += mat[node[prev]][node[j]];
53                                      if (maxj == -1 || dist[node[maxj]] < dist[node[j]])
54                                      maxj = j;
55                                }
56            }
57        }
58
59    }
60    return;
61}
62
63int main()
64 {
65    int n, m, a, b, v;
66    while (scanf("%d%d", &n, &m) != EOF)
67    {
68        res = (1 << 29);
69        memset(mat, 0, sizeof (mat));
70        while (m--)
71    {
72            scanf("%d%d%d", &a, &b, &v);
73            mat[a][b] += v;
74            mat[b][a] += v;
75        }
76        Mincut(n);
77        printf("%d\n", res);
78    }
79    return 0;
80}
81

无向图全局最小边割集相关推荐

  1. Stoer-Wagner无向图全局最小割(hduoj 3691 Nubulsa Expo)

    题意: 一张n个点m条边的无向图,每个点都有一个最大流量,给定起点S,问如何选取终点T,使得最大流最大 翻译一下: 一张n个点m条边的无向联通图,你要割掉一些边使整张图不连通,但是割掉每条边都需要一定 ...

  2. 【HDU - 3002】King of Destruction(无向图全局最小割,SW算法,模板题)

    题干: Zhou xingxing is the successor of one style of kung fu called "Karate Kid".he is falli ...

  3. POJ - 2914 Minimum Cut(全局最小割-Stoer_Wagner)

    题目链接:点击查看 题目大意:给出一张无向图,要求将其分为两个集合,使得最小割最小 题目分析:算法学习自:https://blog.csdn.net/dingdi3021/article/detail ...

  4. POJ 1966 Cable TV Network【无向图点连通度 最小割 E-K算法求最大流】

    题目描述: 给你一个无向图,问你最少删掉几个点,使这个图成不连通. 解题报告: 概念 (1)一个具有 N 个顶点的图,在去掉任意 k-1 个顶点后 (1<=K<=N) 所得的子图仍连通, ...

  5. poj 2914(stoer_wanger算法求全局最小割)

    题目链接:http://poj.org/problem?id=2914 思路:算法基于这样一个定理:对于任意s, t   V ∈ ,全局最小割或者等于原图的s-t 最小割,或者等于将原图进行 Cont ...

  6. 【HDU - 6081】度度熊的王国战略(SW算法,全局最小割)

    题干: Problem Description 度度熊国王率领着喵哈哈族的勇士,准备进攻哗啦啦族. 哗啦啦族是一个强悍的民族,里面有充满智慧的谋士,拥有无穷力量的战士. 所以这一场战争,将会十分艰难. ...

  7. hdu3491 最小点割集(无向图求最小点割集通用方法)

    无向图最小点割集,确定起点S,终点T.每个点都有自己的点权值vi,求最小点权和的割点集,使得S无法到达T. 解法:将每个点拆分为两个点v和v',之间的权值为vi(单向边),将原图中的每条边赋权值为IN ...

  8. 全局最小割(Stoer-Wagner)

    P5632 [模板]Stoer-Wagner算法 板子收集中 #include <bits/stdc++.h> #define inf 0x7fffffff #define ll long ...

  9. 【网络流】学习笔记:一次理解网络流!

    一.从概念入手 网络流用于解决流量问题 网络流:所有弧上流量的集合f={f(u,v)},称为该容量网络的一个网络流. 定义:带权的有向图G=(V,E),满足以下条件,则称为网络流图(flow netw ...

最新文章

  1. 反杀人类、拯救机器狗,被虐士兵机器人化身终结者!这是“波士屯动力”的最新力作...
  2. Python初学问题-if/else的运用
  3. bzoj 1031 [JSOI2007]字符加密Cipher 后缀数组
  4. jsf标签_多个动态包含一个JSF标签
  5. 【渝粤题库】广东开放大学 民事诉讼法 形成性考核
  6. 前端学习(547):node的系统模块fs
  7. C++ mutable 的用法
  8. wps xml转换表格_wps手机版下载-WPS Office 安卓版v12.9.2
  9. JavaCV的摄像头实战之六:保存为mp4文件(有声音)
  10. 微信小程序开发工具基本介绍入门级(备忘)
  11. netware php_服务器_如何在 Netware 服务器中安装多块网卡,如果网络在扩大时服务器只装 - phpStudy...
  12. Deformable DETR进行目标检测,解决size mismatch问题
  13. 如何将多个excel合并成一个
  14. SUPER-VLAN
  15. R可视化会出现的问题1
  16. 计算机电脑照片大小,电脑上怎么调整照片kb
  17. Cocos2d-x 窗口大小调整
  18. 内燃机设计 程序2.0
  19. 【计算机网络】谢希仁笔记 运输层
  20. PDF如何在线转换成PPT呢?

热门文章

  1. 从刘老师的进化的力量到有感,疫情阶段如何弯道超车
  2. Leetcode 1144
  3. MX25上SD卡的插拨检测机制
  4. 自制时代,七种自动化的基础架构即代码工具(IaC)
  5. robotframework调用python类方法_【RF基础】RF调用Python函数基础
  6. 电脑休眠后无法唤醒的解决办法!
  7. 第四篇:ROS常用命令行指令【重点】
  8. HTML+CSS制作人物介绍卡片效果
  9. 量化投资学习——股票价格估值模型
  10. 你知道数据在内存中是如何存储的嘛?