[7.14NOIP模拟4]通讯题解 (Tarjan缩点+贪心)

weixin_30869099

于 2019-07-16 16:28:00 发布

阅读量124

点赞数

原文链接：http://www.cnblogs.com/Rorschach-XR/p/11195741.html

版权

题目描述

“这一切都是命运石之门的选择。”

试图研制时间机器的机关SERN截获了中二科学家伦太郎发往过去的一条短信，并由此得知了伦太郎制作出了电话微波炉（仮）。

为了掌握时间机器的技术，SERN总部必须尽快将这个消息通过地下秘密通讯网络，传达到所有分部。

SERN共有N个部门（总部编号为0），通讯网络有M条单向通讯线路，每条线路有一个固定的通讯花费Ci。

为了保密，消息的传递只能按照固定的方式进行：从一个已知消息的部门向另一个与它有线路的部门传递（可能存在多条通信线路）。我们定义总费用为所有部门传递消息的费用和。

幸运的是，如果两个部门可以直接或间接地相互传递消息（即能按照上述方法将信息由X传递到Y，同时能由Y传递到X），我们就可以忽略它们之间的花费。

由于资金问题（预算都花在粒子对撞机上了），SERN总部的工程师希望知道，达到目标的最小花费是多少。

输入

多组数据，文件以2个0结尾。

每组数据第一行，一个整数N，表示有N个包括总部的部门（从0开始编号）。然后是一个整数M，表示有M条单向通讯线路。

接下来M行，每行三个整数,Xi,Yi,Ci，表示第i条线路从Xi连向Yi，花费为 Ci。

输出

每组数据一行，一个整数表示达到目标的最小花费。

样例输入

样例输出

150
100
50

提示

【样例解释】

第一组数据：总部把消息传给分部1，分部1再传给分部2.总费用：
100+50=150.

第二组数据：总部把消息传给分部1，由于分部1和分部2可以互相传递消息，
所以分部1可以无费用把消息传给2.总费用：100+0=100.

第三组数据：总部把消息传给分部1，最小费用为50.总费用：50.

【数据范围】

对于10%的数据，保证M=N-1

对于另30%的数据，N ≤ 20 ，M ≤ 20

对于100%的数据，N ≤ 50000 ，M ≤ 10^5 ，Ci ≤ 10^5 ，

数据组数 ≤ 5
数据保证一定可以将信息传递到所有部门。

考场上把已经过掉样例的贪心正解一个字一个字删掉还觉得自己贼牛批 23333

大水题只说一下贪心的正确性：数据保证了缩完点联通，所以对于一个点，它的几个父节点(意会一下，不太严谨)都一定能被到达。

所以只考虑这个点的哪条入边代价最小即可缩完点都不用建图直接跑一遍

mlgb打完T3暴力脑子一抽莫名其妙把这玩意证伪了。

似乎有不少人想在有向图上跑Kruskal？hhhh

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<queue>
#define pa pair<int,int>
using namespace std;
typedef long long ll;
const int N=50005,M=1e5+5;
int n,m;
int to[M],nxt[M],head[N],tot;
int len[M];
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9')
    {if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9')
    {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
    return x*f;
}
void add(int x,int y,int z)
{
    to[++tot]=y;
    nxt[tot]=head[x];
    len[tot]=z;
    head[x]=tot;
}
int dfn[N],low[N],ind,st[N],top,bel[N],vis[N],in_edge[N],deg[N],dp[N],scc,root;
void tarjan(int x)
{
    dfn[x]=low[x]=++ind;
    st[++top]=x;vis[x]=1;
    for(int i=head[x];i;i=nxt[i])
    {
        int y=to[i];
        if(!dfn[y])
        {
            tarjan(y);
            low[x]=min(low[x],low[y]);
        }
        else if(vis[y])low[x]=min(low[x],dfn[y]);
    }
    if(dfn[x]==low[x])
    {
        scc++;
        int now=0;
        do
        {
            now=st[top--];
            bel[now]=scc;
            vis[now]=0;
        }
        while(x!=now);
    }
}
void ini()
{
    for(int i=0;i<=m;i++)
        to[i]=nxt[i]=len[i]=0;
    for(int i=0;i<=n;i++)
        head[i]=dfn[i]=bel[i]=st[i]=vis[i]=deg[i]=0,in_edge[i]=0x3f3f3f3f;
    tot=top=ind=scc=root=0;
}
void work()
{
    ini();
    for(int i=1;i<=m;i++)
    {
        int x=read()+1,y=read()+1,z=read();
        add(x,y,z);
    }
    tarjan(1);
    for(int x=1;x<=n;x++)
        for(int i=head[x];i;i=nxt[i])
        {
            int y=to[i];
            if(bel[y]!=bel[x])
                in_edge[bel[y]]=min(in_edge[bel[y]],len[i]);
        }
    ll ans=0;
    for(int i=1;i<=scc;i++)
        if(in_edge[i]!=0x3f3f3f3f)ans+=1LL*in_edge[i];
    cout<<ans<<endl;
}

int main()
{
    while(1)
    {
        n=read();m=read();
        if(!n&&!m)break;
        work();
    }
    return 0;
}