强连通分量缩点Tarjan算法

题目:
听说lcy帮大家预定了新马泰7日游,Wiskey真是高兴的夜不能寐啊,他想着得快点把这消息告诉大家,虽然他手上有所有人的联系方式,但是一个一个联系过去实在太耗时间和电话费了。他知道其他人也有一些别人的联系方式,这样他可以通知其他人,再让其他人帮忙通知一下别人。你能帮Wiskey计算出至少要通知多少人,至少得花多少电话费就能让所有人都被通知到吗?
Input
多组测试数组,以EOF结束。
第一行两个整数N和M(1<=N<=1000, 1<=M<=2000),表示人数和联系对数。
接下一行有N个整数,表示Wiskey联系第i个人的电话费用。
接着有M行,每行有两个整数X,Y,表示X能联系到Y,但是不表示Y也能联系X。
Output
输出最小联系人数和最小花费。
每个CASE输出答案一行。
题目原地址:HDU1827
题意很好理解,解法常用解法强连通分量缩点Tarjan算法。
Tarjan算法是基于对图深度优先搜索的算法,每个强连通分量为搜索树中的一棵子树。搜索时,把当前搜索树中未处理的节点加入一个堆栈,回溯时可以判断栈顶到栈中的节点是否为一个强连通分量。
理解这个算法,但是不会讲出来,推荐一个B站老师讲解的视频 Tarjanj讲解
谈谈我对代码的理解,对于本题,都在注释中

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=2000+10;
const int maxm=100000+10;
const int INF = 0x3f3f3f3f;
vector<int>g[1010];
 int n,m;
int scc_cnt,in[maxn],scc_be[maxn],dfn[maxn],low[maxn],cnt,S[maxn],top,isin[maxn];
void tarjan(int u){
    dfn[u]=low[u]=++cnt;
    S[++top]=u;	//入栈
    isin[u]=1;						//isin数组是用来判断u是否在栈中  
    for(int i=0;i<g[u].size();i++){
        int v=g[u][i];
        if(!dfn[v]){	//如果没有访问过v,就继续递归
            tarjan(v);
            low[u]=min(low[u],low[v]);	//这个是用来回溯时,判断u的指向v的low值是否比自己小
        }
        else if(isin[v]) low[u]=min(low[u],dfn[v]);	//如果v被访问过且这个点在isin栈中存在,说明找到强连通分量了,也就是环了。
    }
    if(low[u]==dfn[u]){	//这个就是用于所有点访问完,进行回溯时用的
        scc_cnt++;	//用来标记每一强连通分量的标记值
        while(1){
            int x=S[top--];	//退栈
            scc_be[x]=scc_cnt;	//标记
            isin[x]=0;
            if(x==u) break;	//整个强连通分量全部退完就结束了
        }
    }
}
void find_scc(){
    memset(dfn,0,sizeof(dfn));
    memset(in,0,sizeof(in));
    memset(low,0,sizeof(low));
    top=-1;
    scc_cnt=0;
    for(int i=1;i<=n;i++)
        if(!dfn[i]) tarjan(i);
}
int main(){
    while(~scanf("%d%d",&n,&m)){
        int w[1010];
        cnt = 0;
        for(int i=1;i<=n;i++){
            scanf("%d",&w[i]);
            g[i].clear();
        }
        int u,v;
        for(int i=0;i<m;i++){
            scanf("%d%d",&u,&v);
            g[u].push_back(v);
        }
        find_scc();
        //这两层for循环的作用就是对所有的强连通分量进行一个处理
        for(int i=1;i<=n;i++){
            for(int j=0;j<g[i].size();j++){
                int v=g[i][j];
                if(scc_be[v]!=scc_be[i])	//如果i和v不在一个强连通分量里,那v就不必在接下来计算最终中出现了,因为在v所属强连通分量里找会更节省
                	in[scc_be[v]]=1;   
            }
        }
        int res=0,num=0;
        for(int i=1;i<=scc_cnt;i++){ 	//这里就容易理解了
            if(!in[i]){
                num++;
                int mi=INF;
                for(int j=1;j<=n;j++){
                    if(scc_be[j]==i)
                        mi=min(mi,w[j]);
                }
                res+=mi;
            }
        }
        cout << endl;
        printf("%d %d\n",num,res);
    }
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值