hdu 3002 全局最小割模板题
2017-03-10 22:40
274 查看
输入注意处理一下
代码:
#include <iostream>
#include <cstdio>
#include <string.h>
using namespace std;
const int maxn=110;
const int inf=0x3f3f3f3f;
int n,m,res,v[maxn],mat[maxn][maxn],dis[maxn];
bool vis[maxn];
int Stoer_Wagner(int n){
int i, j;
int res = inf;
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(){//374MS 1620K
while (scanf("%d%d",&n,&m)==2){
memset(mat,0,sizeof (mat));
int x,y,z;
while(m--){
scanf("%d%d%d",&x,&y,&z);
mat[x+1][y+1] += z;
mat[y+1][x+1] += z;
}
printf("%d\n",Stoer_Wagner(n));
}
}
代码:
#include <iostream>
#include <cstdio>
#include <string.h>
using namespace std;
const int maxn=110;
const int inf=0x3f3f3f3f;
int n,m,res,v[maxn],mat[maxn][maxn],dis[maxn];
bool vis[maxn];
int Stoer_Wagner(int n){
int i, j;
int res = inf;
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(){//374MS 1620K
while (scanf("%d%d",&n,&m)==2){
memset(mat,0,sizeof (mat));
int x,y,z;
while(m--){
scanf("%d%d%d",&x,&y,&z);
mat[x+1][y+1] += z;
mat[y+1][x+1] += z;
}
printf("%d\n",Stoer_Wagner(n));
}
}
相关文章推荐
- HDU 3002 King of Destruction(全局(无向图)最小割(SW))
- hdu 3002 King of Destruction(全局最小割,StoerWagner算法)
- 最小生成树 HDU 各种畅通工程的题,prim和kru的模板题
- 最小割集 hdu 3002 ( 模板)
- HDU 3691 Nubulsa Expo(全局最小割Stoer_Wagn)
- hdu-5251(矩形最小覆盖模板)
- 最小生成树,POJ和HDU几道题目的解题报告(基于自己的模板)
- HDU 3691 全局最小割
- 模板题 + KMP + 求最小循环节 --- HDU 3746 Cyclic Nacklace
- hdu 1814 2-sat 输出字典最小的和任意序列的 模板题
- hdu 1814 2-sat 输出字典最小的和任意序列的 模板题
- 全局最小割模板
- HDU 3002 King of Destruction(SW求最小割的值)
- 最小生成树,POJ和HDU几道题目的解题报告(基于自己写的模板)
- poj 2914&&hdu 3002 全局最小割Stoer-Wagner算法模板
- HDU 2121 Ice_cream’s world II 最小树形图 模板
- HDU 3691 Nubulsa Expo(全局最小割)(10福州B题)#by zh
- HDU 3691 Nubulsa Expo(全局最小割Stoer-Wagner算法)
- hdu 4009 最小树形图模板题朱刘算法
- hdu 1879 继续畅通工程 最小生成树入门模板题