poj Minimum( CutStoer Wagner算法)
Minimum Cut
题目:
给出一张图。要求你删除最小割权和图。
算法分析:
//////////////////// 转载 ———
ylfdrib //////////////////////////////////////////////
一个无向连通网络,去掉一个边集能够使其变成两个连通分量则这个边集就是割集;
最小割集当然就权和最小的割集。
能够用最小分割最大流定理:
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.
看起来非常easy, 每次像做最大生成树一样选最大"边"(注意, 这里事实上不是边, 而是已经累计的权值之和, 就当是加权的度好了), 然后把最后进入的两个点缩到一块就能够了. 合并点最多有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).
///////////////////////////////////////////////////////////////////////////////////////////////////
#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring> const int INF = 1 << 30;
const int MAXN = 500 + 10;
int vst[MAXN];
int wet[MAXN];
int combine[MAXN];
int map[MAXN][MAXN]; int S,T,minCut,N; //最大生成树
void Search(){
int i,j,Max,x,y;
memset(vst,0,sizeof(vst));
memset(wet,0,sizeof(wet));
S = T = -1;
for(i = 0;i < N;++i){
Max = -INF;
for(j = 0;j < N;++j){
if(!combine[j] && !vst[j] && wet[j] > Max)
Max = wet[x = j];
}
if(T == x) return; //有环
S = T; T = x;
minCut = Max;
vst[x] = 1;
for(j = 0;j < N;++j){
if(!combine[j] && !vst[j]){
wet[j] += map[x][j];
}
}
}
} int Stoer_Wagner(){
int i,j;
memset(combine,0,sizeof(combine));
int ans = INF;
for(i = 0;i < N - 1;++i){
Search();
if(minCut < ans) ans = minCut; //最小割
if(0 == ans) return 0;
combine[T] = 1;
for(j = 0;j < N;++j){ //合并
if(!combine[j]){
map[S][j] += map[T][j];
map[j][S] += map[j][T];
}
}
}
return ans;
} int main()
{
int a,b,c,M;
while(~scanf("%d%d",&N,&M)){
memset(map,0,sizeof(map));
while(M--){
scanf("%d%d%d",&a,&b,&c);
map[a][b] += c;
map[b][a] += c;
} printf("%d\n",Stoer_Wagner());
}
return 0;
}
poj Minimum( CutStoer Wagner算法)的更多相关文章
- POJ 2914 Minimum Cut Stoer Wagner 算法 无向图最小割
POJ 2914 题意:给定一个无向图 小于500节点,和边的权值,求最小的代价将图拆为两个联通分量. Stoer Wagner算法: (1)用类似prim算法的方法求"最大生成树" ...
- POJ 2914 - Minimum Cut - [stoer-wagner算法讲解/模板]
首先是当年stoer和wagner两位大佬发表的关于这个算法的论文:A Simple Min-Cut Algorithm 直接上算法部分: 分割线 begin 在这整篇论文中,我们假设一个普通无向图G ...
- POJ 1258 Agri-Net(Prim算法求解MST)
题目链接: http://poj.org/problem?id=1258 Description Farmer John has been elected mayor of his town! One ...
- POJ Minimum Cut
Minimum Cut Time Limit: 10000MS Memory Limit: 65536K Total Submissions: 9302 Accepted: 3902 Case ...
- POJ 矩阵相乘 (随机化算法-舍伍德(Sherwood))
周三的算法课,主要讲了随机化算法,介绍了拉斯维加斯算法,简单的理解了为什么要用随机化算法,随机化算法有什么好处. 在处理8皇后问题的时候,穷举法是最费时的,回朔比穷举好点,而当数据量比较大的时候,如1 ...
- poj 2253 Frogger dijkstra算法实现
点击打开链接 Frogger Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 21653 Accepted: 7042 D ...
- poj - 3259 Wormholes (bellman-ford算法求最短路)
http://poj.org/problem?id=3259 农夫john发现了一些虫洞,虫洞是一种在你到达虫洞之前把你送回目的地的一种方式,FJ的每个农场,由n块土地(编号为1-n),M 条路,和W ...
- POJ 3041.Asteroids-Hungary(匈牙利算法)
Asteroids Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 23963 Accepted: 12989 Descr ...
- POJ 1861 Network (Kruskal算法+输出的最小生成树里最长的边==最后加入生成树的边权 *【模板】)
Network Time Limit: 1000MS Memory Limit: 30000K Total Submissions: 14021 Accepted: 5484 Specia ...
随机推荐
- Even Tree 小议
原题链接:https://www.hackerrank.com/challenges/even-tree/problem 思路:把树还原出来,对每个结点,计算以该结点为根的子树结点数.子树结点数为偶数 ...
- EOJ 3242 重复数
重复数 Time limit per test: 1.0 seconds Time limit all tests: 1.0 seconds Memory limit: 256 megabytes 有 ...
- Windows系统安装MySQL
在Windows中安装mysql不够幸运的话,会遇到相当多的坑,当然这也算是一种财富吧,让自己碰到问题去查找解决方案.有时候不是一时半会就可以解决的.有同学说过安装mysql安装两天还没有装上.不用担 ...
- 项目实战2—LVS负载均衡
负载均衡集群企业级应用实战-LVS 实现基于LVS负载均衡集群的电商网站架构 随着业务的发展,网站的访问量越来越大,网站访问量已经从原来的1000QPS,变为3000QPS,网站已经不堪重负,响应缓慢 ...
- oracle数据库管理系统常见的错误(一)
oracle数据库管理系统常见的错误之一如下: Listener refused the connection with the following error:ORA-12519, TNS:no a ...
- vue新手入门——vue-cli搭建
首先说明,以下内容vue官网都有文档,如果觉得麻烦啰嗦,请移步至 安装-vue.js . 准备工作: 1.下载并安装node环境,一般情况下安装好node之后,npm也会安装好.具体安装的话,百度大概 ...
- Getting Started With setuptools and setup.py
https://pythonhosted.org/an_example_pypi_project/setuptools.html http://www.ianbicking.org/docs/setu ...
- MySQL 2006 超时
http://ronaldbradford.com/blog/sqlstatehy000-general-error-2006-mysql-server-has-gone-away-2013-01-0 ...
- TensorFlow简易学习[2]:实现线性回归
上篇介绍了TensorFlow基本概念和基本操作,本文将利用TensorFlow举例实现线性回归模型过程. 线性回归算法 线性回归算法是机器学习中典型监督学习算法,不同于分类算法,线性回归的输出是整个 ...
- [最短路]P1828 香甜的黄油 Sweet Butter
题目描述 农夫John发现做出全威斯康辛州最甜的黄油的方法:糖.把糖放在一片牧场上,他知道N(1<=N<=500)只奶牛会过来舔它,这样就能做出能卖好价钱的超甜黄油.当然,他将付出额外的费 ...