POJ 1679 The Unique MST (两种思路的次小生成树)

https://vjudge.net/problem/POJ-1679

次小生成树模板题

思路1:由于数据比较小,我们把最小生成树每条边记录下来,之后进行k(k = 最小生成树边数)次最小生成树算法,每次都去掉一条边,判断得到的最小是是否等于最小生成树,这样必然会很慢

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<set>
#include<map>
#include<queue>
#include<cmath>
#define ll long long
#define mod 1000000007
#define inf 0x3f3f3f3f
using namespace std;
struct node
{
    int u, v, w;
}a[10005];
int pre[105];
int miny[10005];
int n, m;
void init()
{
    for(int i = 1; i <= n; i ++)
        pre[i] = i;
}
int find(int x)
{
    if(x == pre[x]) return x;
    else return find(pre[x]);
}
int join(int a, int b)
{
    int i = find(a);
    int j = find(b);
    if(i == j) return 0;
    else
    {
        pre[i] = j;
        return 1;
    }
}
bool cmp(node a, node b)
{
    return a.w < b.w;
}
int main()
{
    int t;
    scanf("%d", &t);
    while(t --)
    {
        scanf("%d%d", &n, &m);
        init();
        for(int i = 1; i <= m; i ++)
            scanf("%d%d%d", &a[i].u, &a[i].v, &a[i].w);
        sort(a + 1, a + 1 + m, cmp);
        int ans = 0;
        int tot = 0;
        int cnt = 0;
        for(int i = 1; i <= m; i ++)
        {
            if(join(a[i].u, a[i].v))
            {
                ans += a[i].w;
                miny[++ cnt] = i;
                tot ++;
            }
            if(tot == n - 1) break;
        }
        bool uni = 0;
        int ans1 = ans;
        for(int k = 1; k <= cnt; k ++)
        {
            init();
            tot = 0, ans = 0;
            bool flag = 0;
            for(int i = 1; i <= m; i ++)
            {
                if(miny[k] == i) continue;
                if(join(a[i].u, a[i].v))
                {
                    ans += a[i].w;;
                    tot ++;
                }
                if(tot == n - 1)
                {
                    flag = 1;
                    break;
                }
            }
            if(flag && ans == ans1)
            {
                uni = 1;
                break;
            }
        }
        if(! uni)
        {
            printf("%d\n", ans1);
        }
        else printf("Not Unique!\n");
    }
    return 0;
}

思路二:Kruskal次小生成树算法,和最小生成树同理,只不过我们每次都记录下来每两个点之间的最大距离,这样删去这条边再加上一条边,必然构成生成树(证明略)

对于一棵已经求出的最小生成树
枚举每一条不在最小生成树上的边
并把这条边加入最小生成树
设这条边连接的两个节点为u和v
这时树上u到v的路径会出现回路
所以我们需要删掉u到v路径上的一条边使其重新变成一颗树
这时生成树权值的增量就是新加入的边权-删掉的边权
为了使增量最小,我们要删去的自然是u到v路径中权值最大的

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<set>
#include<map>
#include<queue>
#include<cmath>
#define ll long long
#define mod 1000000007
#define inf 0x3f3f3f3f
using namespace std;
struct node
{
    int u, v, w;
}a[10005];
int pre[105];
int miny[10005];
int maxy[105][105];
vector<int> G[105];
bool vis[10005];
int n, m;
int ans;
void init()
{
    for(int i = 1; i <= n; i ++)
    {
        G[i].clear();
        G[i].push_back(i);
        pre[i] = i;
    }
    memset(vis, 0, sizeof(vis));
    memset(maxy, -1, sizeof(maxy));
    ans = 0;
}
int find(int x)
{
    if(x == pre[x]) return x;
    else return find(pre[x]);
}
bool cmp(node a, node b)
{
    return a.w < b.w;
}
bool Kruskal()
{
    init();
    int tot = 0;
    sort(a + 1, a + 1 + m, cmp);
    for(int i = 1; i <= m; i ++)
    {
        int u = find(a[i].u);
        int v = find(a[i].v);
        if(u != v)
        {
            tot ++;
            ans += a[i].w;
            vis[i] = 1;
            pre[u] = v;
            int sz1 = G[u].size(), sz2 = G[v].size();
            for(int j = 0; j < sz1; j ++)
                for(int k = 0; k < sz2; k ++)
                    maxy[G[u][j]][G[v][k]] = maxy[G[v][k]][G[u][j]] = a[i].w;
            for(int j = 0; j < sz1; j ++)
                G[v].push_back(G[u][j]);
        }
        if(tot == n - 1) break;
    }
    int sec = inf;
    for(int i = 1; i <= m; i ++)
    {
        if(! vis[i])
        {
            sec = min(sec, ans - maxy[a[i].u][a[i].v] + a[i].w);
        }
    }
    if(ans == sec) return 0;
    return 1;
}
int main()
{
    int t;
    scanf("%d", &t);
    while(t --)
    {
        scanf("%d%d", &n, &m);
        init();
        for(int i = 1; i <= m; i ++)
            scanf("%d%d%d", &a[i].u, &a[i].v, &a[i].w);
        if(Kruskal()) printf("%d\n", ans);
        else printf("Not Unique!\n");
    }
    return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值