转载地址:http://blog.csdn.net/xdu_truth/article/details/8104721

题意:题给出一个无向图和一个源点,让你求从这个点出发到某个点最大流的最小值。由最小割最大流定理,最全局最小割就可以了,但是因为我没有模版所以这题就悲剧了……之前找到了一个模版改半天没改好,后来发现只要改一个地方就行了,AC了发现速度还挺快的,代码也不长,嗯,这个模版要收藏

这个代码跑了600+ms

下面那个模板跑了1400+ms

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int MAXN = 305;
const int MAXV = 0x3F3F3F3F;
int n,m,v[MAXN],mat[MAXN][MAXN],dis[MAXN];
bool vis[MAXN];
int res,s;
int Stoer_Wagner(int n)
{
int i, j;
int res = MAXV;
for (i = 0; i < n; i++)
v[i] = i+1;//初始化第i个结点就是i
while (n > 1)
{
int maxp = 1,prev = 0;
for (i = 1; i < n; i++) //初始化到已圈集合的割大小,并找出最大距离的顶点
{
dis[v[i]] = mat[v[0]][v[i]];
if (dis[v[i]] > dis[v[maxp]])
maxp = i;
}
memset(vis, 0, sizeof(vis));
vis[v[0]] = true;
for (i = 1; i < n; i++)
{
if (i == n - 1) //只剩最后一个没加入集合的点,更新最小割
{
res = min(res,dis[v[maxp]]);
for (j = 0; j < n; j++) //合并最后一个点以及推出它的集合中的点
{
mat[v[prev]][v[j]] += mat[v[j]][v[maxp]];
mat[v[j]][v[prev]] = mat[v[prev]][v[j]];
}
v[maxp] = v[--n];//第maxp个节点去掉,第n个节点变成第maxp个
}
vis[v[maxp]] = true;
prev = maxp;
maxp = -1;
for (j = 1; j < n; j++)
if (!vis[v[j]]) //将上次求的maxp加入集合,合并与它相邻的边到割集
{
dis[v[j]] += mat[v[prev]][v[j]];
if (maxp == -1 || dis[v[maxp]] < dis[v[j]])
maxp = j;
}
}
}
return res;
}
int main()
{
//freopen("input.txt","r",stdin);
while (scanf("%d%d%d", &n, &m, &s),n||m||s)
{
memset(mat,0,sizeof (mat));
int x,y,z;
while (m--)
{
scanf("%d%d%d",&x,&y,&z);
mat[x][y] += z;
mat[y][x] += z;
}
printf("%d\n",Stoer_Wagner(n));
}
}

因为他从1点出发,1这个点必定属于一个集合,最小割相当于割掉一些边,分成两个不同的集合,所以用全局最小割可解

#include<stdio.h>
#include<string.h>
#include<iostream>
#define inf 0x3fffffff
#define N 510
int n,ma[N][N],combine[N];
int seach(int &s,int &t) {
int vis[N],i,j,tm,maxx,w[N];
memset(vis,0,sizeof(vis));
memset(w,0,sizeof(w));
tm=10000;
for(i=0;i<n;i++) {
maxx=-inf;
for(j=0;j<n;j++)
if(!vis[j]&&!combine[j]&&maxx<w[j]) {//找到最大值
maxx=w[j];
tm=j;
}
if(t==tm) {return w[t];}//最后t和tm相等
vis[tm]=1;
s=t;t=tm;
for(j=0;j<n;j++)//
if(!vis[j]&&!combine[j])
w[j]+=ma[t][j];
}
return w[t];
}
int mincut() {
int mi=inf,ans,i,s,t,j;
memset(combine,0,sizeof(combine));
for(i=0;i<n-1;i++) {//只需要找n-1次
s=-1;t=-1;
ans=seach(s,t);//找到t和t的前一个点s
combine[t]=1;//移除
if(ans<mi)mi=ans;
for(j=0;j<n;j++) {//将t合并到s点
ma[s][j]+=ma[t][j];
ma[j][s]+=ma[j][t];
}
}
return mi;
}
int main() {
int m,i,j,k;
while(scanf("%d%d%d",&n,&m,&k),n||m||k) {
memset(ma,0,sizeof(ma));
while(m--) {
scanf("%d%d%d",&i,&j,&k);
i--;j--;
ma[i][j]+=k;ma[j][i]+=k;
}
printf("%d\n",mincut());
}
return 0;
}

hdu 3691最小割将一个图分成两部分的更多相关文章

  1. hdu 3657 最大点权独立集变形(方格取数的变形最小割,对于最小割建图很好的题)

    转载:http://blog.csdn.net/cold__v__moon/article/details/7924269 /* 这道题和方格取数2相似,是在方格取数2的基础上的变形. 方格取数2解法 ...

  2. hdu 5076 最小割灵活运用

    这意味着更复杂的问题,关键的事实被抽象出来:每个点,能够赋予既有的值(挑两个一.需要选择,设定ai,bi). 寻找所有和最大.有条件:如果两个点同时满足: 1,:二进制只是有一个不同之处.  2:中的 ...

  3. hdu 3657 最小割的活用 / 奇偶方格取数类经典题 /最小割

    题意:方格取数,如果取了相邻的数,那么要付出一定代价.(代价为2*(X&Y))(开始用费用流,敲升级版3820,跪...) 建图:  对于相邻问题,经典方法:奇偶建立二分图.对于相邻两点连边2 ...

  4. Being a Hero (hdu 3251 最小割 好题)

    Being a Hero Time Limit: 20000/10000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) T ...

  5. bzoj 3232 圈地游戏 —— 01分数规划+最小割建图(最大权闭合子图)

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3232 心烦意乱的时候调这道题真是...越调越气,就这样过了一晚上... 今天再认真看看,找出 ...

  6. hdu 4289(最小割)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4289 思路:求最小花费,最小割应用,将点权转化为边权,拆点,(i,i+n)之间连边,容量为在城市i的花 ...

  7. hdu 3657 最小割(牛逼!!!!)总算理解了

    <strong></strong> 转载:http://blog.csdn.net/me4546/article/details/6662959 加颜色的太棒了!!! 在网上看 ...

  8. [HDU 3521] [最小割] Being a Hero

    题意: 在一个有向图中,有n个点,m条边$n \le 1000 \And \And  m \le 100000$ 每条边有一个破坏的花费,有些点可以被选择并获得对应的金币. 假设一个可以选的点是$x$ ...

  9. Game HDU - 3657(最小割)

    Game Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submis ...

随机推荐

  1. dalvik.system.VMRuntime 隐藏api的迷惑

    [Android UI界面]关于dalvik.system.VMRuntime 的 使用迷惑 我也遇到了相同问题.不知楼主现在解决了没有? 回答1: [Android UI界面]关于dalvik.sy ...

  2. B. Hongcow Solves A Puzzle

    http://codeforces.com/contest/745/problem/B 题目要求的是,给定一个图形,要求里面判断是否有矩形,且仅有一个 就是 XXX.... XXX...X 是不行的, ...

  3. Python学习 Day 3 字符串 编码 list tuple 循环 dict set

    字符串和编码 字符 ASCII Unicode UTF-8 A 1000001 00000000 01000001 1000001 中 x 01001110 00101101 11100100 101 ...

  4. nodejs的学习

    nodejs 就是使用js来编写服务端的程序.它的特性是(单线程   速度快   耗内存多  异步   事件驱动) ( 一些技术的解决方案:默认情况下是 1.不支持多核,可以使用cluster 进行解 ...

  5. day21-4 菱形继承问题(类的查找顺序)

    目录 菱形继承问题 经典类(了解) 新式类 mro方法 菱形继承问题 在Python中子类可以同时继承多个父类,如A(B,C,D) 如果继承关系为非菱形结构,则会按照先找B这一条分支,然后再找C这条分 ...

  6. SpringCloud 微服务框架

    学习资源:https://ke.qq.com/course/280057 知识体系分为以下几点: 1)使用Eureka搭建注册中心,包括 服务生产者.服务消费者(也称服务注册与发现): Zookeep ...

  7. 时钟周期 VS 机器周期

    时钟周期vs机器周期 Clock cycle The speed of a computer processor, or CPU, is determined by the clock cycle, ...

  8. thinkphp5将一条数据以toArray()放入session

    直接将找出的数据赋予变量,列:$admin: 存入session:Session::set('user_info',$admin->toArray()); 如果数据中有时间字段: 如若有时间字段 ...

  9. Jedis集成到项目中

    Jedis整合到项目中,就可以在项目中使用redis了,作为Java程序狗,这个可以会,贴代码了,不截图了,哈哈 一.maven中的pom.xml中添加依赖 <dependency> &l ...

  10. react-native Socket Event 在控制台的输出

    在XCode中运行react-native 的时候,避免不了的要查看日志信息 ,但是react-native中的Socket的日志简直是太多了,往往是刚看到自己想要看到的信息的时候,瞬间就被最新的日志 ...