[JZOJ4686] 通讯 [HDU3072] Intelligence System

BAJim_H

于 2016-08-12 21:31:05 发布

阅读量348

点赞数

分类专栏：题解 ————最小生成树 ————强连通分量

本文链接：https://blog.csdn.net/hzj1054689699/article/details/52195447

版权

题解同时被 3 个专栏收录

353 篇文章 0 订阅

订阅专栏

————最小生成树

6 篇文章 0 订阅

订阅专栏

————强连通分量

3 篇文章 0 订阅

订阅专栏

Summary

给出一个带权有向图。

现在要传递一个消息，已知消息的点可以沿边传向未知的点，费用是这条边的边权。

但是，由于某些神奇的原因，如果两个点能互相到达，那么它们之间传递的费用为0。

一开始只有点0已知消息，问至少需多少费用使所有点已知？

Solution

显然，每个强连通分量中传递费用都是0。

然后我们可以用最小树形图（就是有向图的最小生成树）来做，但是怎么判环呢？

可以缩点！

将每个强连通分量缩成一个点，原图就变成了一个 $DAG$ （有向无环图）

然后再来做最小树形图

显然选出来每个点有且只有一条入边。

贪心的想，每个点都选边权最小的入边不就可以了？

这就是最小树形图的做法，由于是 $DAG$ ，就省去了判环的环节

这里注意，我们并不需要重新构图，仅需要将原来的强连通分量染上色，选边的时候特判一下是否两点颜色不同，并且每种颜色只能选一个入边。

Code

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fod(i,a,b) for(i=a,i<=b,i--)
#define MAXN 50005
#define MAXM 100005
using namespace std;
struct node
{
    int x,y,z;
}a[5*MAXM];
int a1[5*MAXN][2],n,m,dfn[5*MAXN],low[5*MAXN],dep,ans[5*MAXN],st[5*MAXN],top,cl[5*MAXN],clm,b1[5*MAXN];
bool bz[5*MAXN];
bool cmp(node x,node y)
{
    return x.x<y.x;
}
void tarjan(int k)
{
    dfn[k]=low[k]=++dep;
    int i;
    if (a1[k][0]!=0) 
    fo(i,a1[k][0],a1[k][1])
    {
        int p=a[i].y;
        if (dfn[p]==0) 
        {
            st[++top]=p;
            bz[p]=1;
            tarjan(p);
            low[k]=min(low[k],low[p]);
        }
        else if (bz[p]==1) low[k]=min(low[k],dfn[p]); 
    }
    if (low[k]==dfn[k])
    {
        clm++;
        while (st[top]!=k)
        {
            bz[st[top]]=0;
            cl[st[top]]=clm;
            st[top--]=0;
        }
        bz[st[top]]=0;
        cl[st[top]]=clm;
        st[top--]=0;
    }
}
int main()
{
    cin>>n>>m;
    int i,j;
    while (n!=0)
    { 
        int m1=m;
        m=0;
        fo(i,1,m1) 
        {
            int x,y,z;
            scanf("%d%d%d",&x,&y,&z);
            if (x!=y)
            {
                a[++m].x=x;
                a[m].y=y;
                a[m].z=z;
            } 
        }
        sort(a+1,a+m+1,cmp);
        fo(i,2,m)
        {
            if (a[i].x!=a[i-1].x)
            {
                a1[a[i].x][0]=i;
                a1[a[i-1].x][1]=i-1; 
            }
        }
        a1[a[1].x][0]=1;
        a1[a[m].x][1]=m;
        dep=0;
        top=1;
        st[1]=0; 
        memset(dfn,0,sizeof(dfn));
        memset(st,0,sizeof(st));
        memset(bz,0,sizeof(bz));
        bz[0]=1;
        tarjan(0);
        memset(ans,107,sizeof(ans));
        long long s=0;
        memset(b1,107,sizeof(b1));
        fo(i,1,m)
        {
            if (cl[a[i].x]!=cl[a[i].y]) b1[cl[a[i].y]]=min(b1[cl[a[i].y]],a[i].z);
        }
        fo(i,1,clm) if (b1[i]!=1802201963&&cl[0]!=i) s+=b1[i];
        cout<<s<<endl;
        scanf("%d%d",&n,&m);
        memset(a1,0,sizeof(a1)); 
    }
}