poj1679The Unique MST

题目大意:

次小生成树。

解题思路:

假设存在最小生成树T, T*是一颗生成树,但不一定是最小的,则T*必定可以由T更换一条边(x,y)的得来,且这条边必定比原来的要长。

要得到次小生成树,就得更换一条边,且这条边要尽可能的接近但不等于原边的长度,这样的话我们可以枚举每一条边,求一个大于但不等于MST的权值,这样的做法可取,但不是最优的。

上一个图



当我们向一颗生成树中添加一条新边的时候就会形成一个环,也就是多了一条连接如上两个集合的路径,那么我们要做的就是删除这两个集合中的一条边,保证集合不会被拆散,删除这条边之后的权值要是最小的。那么就简单了,直接删除一条权值最大的边即可。

代码:

#include<stdio.h>
#include<algorithm>
#define N 109
#define M 20009
const int inf= 1<<31-1;
int n,m;
//并查集部分
int ut[N];
int find(int x)
{
    int y=x,k;
    while(x!=ut[x])
               x=ut[x];
    while(y!=ut[y]){
        k=y;
        y=ut[y];
        ut[k]=x;
    }
    return x;
}
void merge(int x,int y)
{
    if(x!=y)
    {
        ut[x]=ut[y];
    }
}
//边集
struct node
{
   int a,b;
   int w;
   bool select;
}edg[M];
bool cmp(node x,node y)
{
    if(x.w!=y.w)
        return x.w<y.w;
    if(x.a!=y.a)
        return x.a<y.a;
    return x.b<y.b;
};
//链式前向星保存图的两个部分
struct node1
{
   int to;
   int next;
};
node1 link[N];//边数组
int h;
int head[N],end[N];//
int dis[N][N];
//克鲁斯卡尔修改版
void skruskal()
{
    int k=0,i,x,y,w,v;
    for(h=0;h<n;h++)
    {
        link[h].to=h+1;
        link[h].next=head[h+1];
        end[h+1]=h;
        head[h+1]=h;
    }
    std::sort(edg+1,edg+1+m,cmp);
    for(i=1;i<m;i++)
    {
        if(k==n-1)
              break;
        if(edg[i].w<0)
              continue;
        x=find(edg[i].a);
        y=find(edg[i].b);
        if(x!=y)//两个新的不想交的集合
        {
            for(w=head[x];w!=-1;w=link[w].next)
            {
                for(v=head[y];v!=-1;v=link[v].next)
                {
                    dis[link[w].to][link[v].to]=dis[link[v].to][link[w].to]=edg[i].w;
                }
            }
        link[end[y]].next=head[x];
        end[y]=end[x];
        merge(x,y);
        k++;
        edg[i].select=true;
      }
    }
}
void initi()
{
    int i;
    for(i=0;i<N;i++){
    ut[i]=i;
    head[i]=-1;end[i]=-1;
    }
}
int min(int a,int b)
{
    return a<b?a:b;
}
int main()
{
    int t,i,mst,semst;
    scanf("%d",&t);
    while(t--)
    {
        initi();
        scanf("%d%d",&n,&m);
        for(i=1;i<=m;i++)
            scanf("%d%d%d",&edg[i].a,&edg[i].b,&edg[i].w);
        mst=0;
        skruskal();
        for(i=1;i<=m;i++)
        {
            if(edg[i].select)
                mst+=edg[i].w;
        }
        semst=inf;
        for(i=1;i<=m;i++)
        {
            if(!edg[i].select)
            {
                semst=min(semst,mst+edg[i].w-dis[edg[i].a][edg[i].b]);

            }
        }
        if(mst!=semst)
            printf("%d\n",mst);
        else
            puts("Not Unique!");
    }
}


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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值