这是一个全局最小割的题:
#include <iostream>
#include <cstring>
#include <cstdio>
#include <queue>
using namespace std ;
const int INF = 999999999 ;
const int MAXN = 505 ;
int map[MAXN][MAXN] , low[MAXN] , vis[MAXN] , last , pre , N , flag ;
bool bfs() ;
void prim( int ) ;
int stoerwagner( int ) ;
int main()
{
int m , n , start , end , w , ans ;
while( scanf( "%d%d" , & n , & m ) != EOF )
{
N = n ;
memset( map , 0 , sizeof( map ) ) ;
while( m -- )
{
scanf( "%d%d%d" , & start , & end , & w ) ;
map[start][end] += w ;
map[end][start] += w ;
}
flag = bfs() ;
if( flag ){
ans = stoerwagner( n ) ;
cout << ans << endl ;
}
else cout << 0 << endl ;
}
return 0 ;
}
bool bfs()
{
flag = 0 ;
memset( vis , 0 , sizeof( vis ) ) ;
queue < int > q ;
vis[0] = 1 ;
q.push( 0 ) ;
while( !q.empty() )
{
int temp = q.front() ;
q.pop() ;
for( int i = 0 ; i < N ; i ++ )
if( map[temp][i] > 0 && !vis[i] )
{
q.push( i ) ;
vis[i] = 1 ;
}
}
for( int i = 0 ; i < N ; i ++ )
if( !vis[i] ) return false ;
return true ;
}
int stoerwagner( int n )
{
int Min = INF , temp , tmp ;
while( n > 1 )
{
prim( n ) ;
temp = 0 ;
for( int j = 0 ; j < N ; j ++ )
if( map[j][last] ) temp += map[j][last] ;
if( temp < Min && temp > 0 ) Min = temp ;
n -- ;
map[pre][last] = map[last][pre] = 0 ;
temp = min( pre , last ) ;
tmp = max( pre , last ) ;
for( int j = 0 ; j < N ; j ++ )
{
if( map[tmp][j] ) {
map[temp][j] += map[tmp][j] ;
map[tmp][j] = 0 ;
}
if( map[j][tmp] ) {
map[j][temp] += map[j][tmp] ;
map[j][tmp] = 0 ;
}
}
}
return Min ;
}
void prim( int n )
{
memset( vis , 0 , sizeof( vis ) ) ;
memset( low , 0 , sizeof( low ) ) ;
low[0] = INF ;
for( int i = 1 ; i <= n ; i ++ )
{
int Max = 0 , mark ;
for( int j = 0 ; j < N ; j ++ )
if( !vis[j] && low[j] > Max ) {
Max = low[j] ;
mark = j ;
}
vis[mark] = 1 ;
if( i == n - 1 ) pre = mark ;
if( i == n ) last = mark ;
for( int j = 0 ; j < N ; j ++ )
if( !vis[j] ) low[j] += map[mark][j] ;
}
}
分享到:
相关推荐
poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题...
poj 3414解题报告poj 3414解题报告poj 3414解题报告poj 3414解题报告
poj 1012解题报告poj 1012解题报告poj 1012解题报告poj 1012解题报告
poj 2329解题报告poj 2329解题报告poj 2329解题报告poj 2329解题报告
poj 1440解题报告 poj 1440解题报告 poj 1440解题报告 poj 1440解题报告
poj 3083解题报告poj 3083解题报告poj 3083解题报告poj 3083解题报告
poj 1659解题报告poj 1659解题报告poj 1659解题报告poj 1659解题报告
poj 3720解题报告poj 3720解题报告poj 3720解题报告poj 3720解题报告
poj1691解题报告 题目来源:http://acm.pku.edu.cn/JudgeOnline/showproblem?problem_id=1691(POJ No.1691) 解法: 搜索
poj2828解题报告,希望能帮到志同道合的算法爱好者
北大poj解题报告,希望能帮到软件工程的同学,每天一道,持之以恒,熟能生巧,与您共勉!
北大ACM1316解题报告
Problem 1061 青蛙的约会 poj解题报告。有源码。可直接提交C++实现
这是我发的第二篇解题报告,写的不怎么的。呵呵!!!!!
ACM Poj Pku 解题报告答案 打包 下载 600多题 史上最全 不是网上乱传的200多题,更不是100多题就挂着10分才能下的题 下了这个 大家也不要浪费分数去下载其它版本的了,基本上都有 共享 一起进步 中国加油 ACMer...
ACM POJ 解题报告北大POJ 大量解题代码
本人整理的POJ解题报告,一共有250道题
北大ACM在线评测系统POJ的题目解题报告。涵盖各种类型的acm题目。值得参考借鉴。打包下载。
2遍dp poj_3613解题报告 poj_3613解题报告
这是我的第三篇结题报告,欢迎大家指正!!!!!!!!!!!!!