1513:受欢迎的牛(DFS&tarjan&强连通分量)

信息学奥赛一本通(C++版)在线评测系统

解析:

tarjan算法,将有向图缩点,转换为拓扑图,并且统计出度为0的强连通分量;

如果这样的强连通分量只有一个,那么答案即为这个强连通分量中点的个数,否则答案为0(因为如果存在多个强连通分量的出度为0,那么其显然不可能被所有牛欢迎)

cpp 复制代码
#include<bits/stdc++.h>
using namespace std;
const int N=1e4+5,M=5e4+5;

int n,m;
int h[N],e[M],ne[M],idx;
int dfn[N],low[N],timestamp;
int stk[N],top;
bool in_stk[N];
int id[N],scc_cnt,sz[N];
int dout[N];	//出度 

void add(int a,int b){
	e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}

void tarjan(int u){
	dfn[u]=low[u]=++timestamp;
	stk[++top]=u;
	in_stk[u]=true;
	for(int i=h[u];~i;i=ne[i]){
		int j=e[i];
		if(!dfn[j]){
			tarjan(j);
			low[u]=min(low[u],low[j]);
		}
		else if(in_stk[j]) low[u]=min(low[u],dfn[j]);
	}
	if(dfn[u]==low[u]){
		int y;
		++scc_cnt;
		do{
			y=stk[top--];
			in_stk[y]=false;
			id[y]=scc_cnt;
			sz[scc_cnt]++;
		}while(y!=u);
	}
}
signed main(){
	scanf("%d%d",&n,&m);
	memset(h,-1,sizeof h);
	for(int i=1;i<=m;i++){
		int a,b;
		scanf("%d%d",&a,&b);
		add(a,b);
	}
	for(int i=1;i<=n;i++) if(!dfn[i]) tarjan(i);
	for(int i=1;i<=n;i++){
		for(int j=h[i];~j;j=ne[j]){
			int k=e[j];
			int a=id[i],b=id[k];
			if(a!=b) dout[a]++;
		}
	}
	int res=0,sum=0;
	for(int i=1;i<=scc_cnt;i++){
		if(!dout[i]){
			res++;
			sum+=sz[i];
			if(res>1){
				sum=0;
				break;
			}
		}
	}
	cout<<sum;
	return 0;
}
相关推荐
小李小李快乐不已1 小时前
二叉树理论基础
数据结构·c++·算法·leetcode
仰泳的熊猫1 小时前
1149 Dangerous Goods Packaging
数据结构·c++·算法·pat考试
_OP_CHEN1 小时前
【算法基础篇】(三十七)图论基础之多源最短路:Floyd 算法吃透所有点对最短路径!
算法·蓝桥杯·图论·算法竞赛·floyd算法·acm/icpc·多源最短路
Web极客码1 小时前
如何选择最适合的内容管理系统(CMS)?
java·数据库·算法
程序员三明治1 小时前
【动态规划】01背包与完全背包问题详解,LeetCode零钱兑换II秒解,轻松解力扣
算法·leetcode·动态规划·java后端·01背包·完全背包·零钱兑换
子木鑫1 小时前
【练习】循环
c语言
自由生长20241 小时前
大数据计算框架-流式计算的Join
算法
IT猿手1 小时前
融合DWA的青蒿素优化算法(Artemisinin Optimization Algorithm, AOA)求解无人机三维动态避障路径规划,MATLAB代码
算法·matlab·无人机
H_z___1 小时前
Codeforces Global Round 31 (Div. 1 + Div. 2) A ~ E
数据结构·算法
Dylan的码园1 小时前
队列与queue
java·数据结构·链表