HDU 3691 Nubulsa Expo(SW算法求最小割)
2013-04-23 11:12
330 查看
题目:给出一个网络流图,给定一个起点,汇点自己定,然后求这个图的最小割。
分析:枚举汇点,时间是三秒,超时!只能用这个办法去求最小割了,我还不是不很懂这个算法,只是理解了一些皮毛,在网上扒个代码,先体验一下,这个算法和prim本质上是两回事,但是实现是很像的,时间复杂度如果不用优先队列优化的话是O(n^3),优化的话是O((n^2)logn)。具体算法的理解我会在晚些时候写总结,下面先把这题的代码贴出来!
代码:
#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
#include <algorithm>
using namespace std;
const int N = 350;
const int INF = 0x7fffffff;
int n, m, S;
int g
, dis
;
bool vis
, d
;
int SW()
{
int p = n, ans = INF, mi, s, t, k, i;
memset( d, 0, sizeof(d) );
while ( --p ) {
memset( vis, 0, sizeof(vis) );
memset( dis, 0, sizeof(dis) );
i = 1;
while ( d[i] ) i++;
vis[i] = 1;
for ( int j = 1; j <= n; ++j )
if ( !d[j] && !vis[j] ) dis[j] = g[i][j];
t = s = i;
for ( ; i <= n; ++i ) {
mi = 0;
for ( int j = 1; j <= n; ++j )
if ( !d[j] && !vis[j] && mi < dis[j] )
mi = dis[k=j];
if ( !mi ) break;
vis[k] = true;
for ( int j = 1; j <= n; ++j )
if ( !d[j] && !vis[j] ) dis[j] += g[k][j];
s = t, t = k;
}
if ( ans > dis[t] ) ans = dis[t];
d[t] = 1;
for ( int j = 1; j <= n; ++j )
if ( !d[j] )
{
g[s][j] += g[t][j];
g[j][s] += g[j][t];
}
}
return ans;
}
int main()
{
while (scanf("%d%d%d", &n, &m, &S) != EOF && ( n || m || S ) ) {
memset ( g, 0, sizeof(g) );
for ( int i = 0, s, e, c; i < m; ++i ) {
scanf("%d%d%d", &s, &e, &c);
g[s][e] += c;
g[e][s] += c;
}
printf("%d\n", SW());
}
}
分析:枚举汇点,时间是三秒,超时!只能用这个办法去求最小割了,我还不是不很懂这个算法,只是理解了一些皮毛,在网上扒个代码,先体验一下,这个算法和prim本质上是两回事,但是实现是很像的,时间复杂度如果不用优先队列优化的话是O(n^3),优化的话是O((n^2)logn)。具体算法的理解我会在晚些时候写总结,下面先把这题的代码贴出来!
代码:
#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
#include <algorithm>
using namespace std;
const int N = 350;
const int INF = 0x7fffffff;
int n, m, S;
int g
, dis
;
bool vis
, d
;
int SW()
{
int p = n, ans = INF, mi, s, t, k, i;
memset( d, 0, sizeof(d) );
while ( --p ) {
memset( vis, 0, sizeof(vis) );
memset( dis, 0, sizeof(dis) );
i = 1;
while ( d[i] ) i++;
vis[i] = 1;
for ( int j = 1; j <= n; ++j )
if ( !d[j] && !vis[j] ) dis[j] = g[i][j];
t = s = i;
for ( ; i <= n; ++i ) {
mi = 0;
for ( int j = 1; j <= n; ++j )
if ( !d[j] && !vis[j] && mi < dis[j] )
mi = dis[k=j];
if ( !mi ) break;
vis[k] = true;
for ( int j = 1; j <= n; ++j )
if ( !d[j] && !vis[j] ) dis[j] += g[k][j];
s = t, t = k;
}
if ( ans > dis[t] ) ans = dis[t];
d[t] = 1;
for ( int j = 1; j <= n; ++j )
if ( !d[j] )
{
g[s][j] += g[t][j];
g[j][s] += g[j][t];
}
}
return ans;
}
int main()
{
while (scanf("%d%d%d", &n, &m, &S) != EOF && ( n || m || S ) ) {
memset ( g, 0, sizeof(g) );
for ( int i = 0, s, e, c; i < m; ++i ) {
scanf("%d%d%d", &s, &e, &c);
g[s][e] += c;
g[e][s] += c;
}
printf("%d\n", SW());
}
}
相关文章推荐
- HDU 3691 Nubulsa Expo(全局最小割Stoer-Wagner算法)
- HDU 3691 Nubulsa Expo (全局最小割)
- Hdu 3691 Nubulsa Expo(无向图最小割)
- HDU 3691 Nubulsa Expo(全局最小割)
- HDU 3691 Nubulsa Expo(全局最小割Stoer_Wagn)
- HDU 3691 Nubulsa Expo(全局最小割)(10福州B题)#by zh
- hdu 1162 Eddy's picture(最小生成树算法)
- hdu1102 Constructing Roads (最小生成树 prim 算法)
- 最小生成树 Prim 算法HDU 2122 Ice_cream’s world III
- HDU - 3691 Nubulsa Expo sw算法求全局最小割
- hdu1054匈牙利算法/最小覆盖点
- HDU 2063 过山车(匈牙利算法 二分图的最小顶点覆盖 二分图最大匹配)
- HDU - 1150 Machine Schedule 最小点覆盖(最大二分图匹配-匈牙利算法)邻接表写法
- HDU 4966 GGS-DDU (最小树形图-朱刘算法)
- HDU 1150 Machine Schedule (最小覆盖,匈牙利算法)
- hdu1150匈牙利算法最小点覆盖
- hdu1102 Constructing Roads(两种基础最小生成树算法)
- HDU 2121 Ice_cream’s world II 无固定点的最小树形图 朱刘算法
- HDU 3002 King of Destruction(全局(无向图)最小割(SW))
- hdu1498匈牙利算法/最小点覆盖