题目描述
朋友圈里有很多人。每只人的梦想是成为最受欢迎的朋友。有 N 个人,有 M对二元关系 (A,B),
告诉你A 认为 B 是受欢迎的。如果 A 认为 B 是受欢迎的,B 认为 C 是受欢迎的,则 A 也认为
C 是受欢迎的。你的任务是计算被其余人都认为是受欢迎的朋友数量。
输入格式:第一行两个正整数 N 和 M,分别表示一共有 N 个人和 M 对二元关系
(1≤N≤10000,1≤M≤50000)。接下来 M 行,每行两个整数 A 和 B,表示 A 认为 B 是受欢
迎的 (1≤A,B≤N)。
输出格式: 一行一个整数,表示答案。
样例输入
4 6
1 2
1 3
1 4
2 3
2 4
3 4
样例输出
1
预备知识
tarjan算法
1,任选一顶点开始进行深度优先搜索(若深度优先搜索结束后仍有未访问的顶点,则再
从中任选一点再次进行)。搜索过程中已访问的顶点不再访问。搜索树的若干子树构
成了图的强连通分量。 顶点按照被访问的顺序存入栈中。从搜索树的子树返回至一个结点时,检查该结点是
2,否是某一强连通分量的根结点(见下)并将其从栈中删除。如果某结点是强连通分量
的根,则在它之前出栈且还不属于其他强连通分量的结点构成了该结点所在的强连通
分量。
解题思路
这道题如果用Floyd算法来做的话最多只能过一部分数据,当数据大的时候二维数组就会爆内存。而且时间复杂度也为O(V^3)比较大。所以我们用tarjan算法来做。首先我们知道题目问题的意思就是求这个图里有多少个点是其他所有点都能到达改点的。很容易就能想到如果我们用缩点来优化的话(及将每个强连通分量看成“一个点”),就会优化许多点。但是即使是这样之后再怎么做呢。是不是还是用Floyd算法来判断两个缩点之间有没有路呢。这样的话虽然会快一点但是如果碰到比较刁钻的数据还是很可能卡住。我们用到强连通之后的图的性质来做将会很简单。
那么在经过tarjan算法建立缩点图之后有什么性质呢。
1,所有的建立在缩点之间的边都是单向边且图中没有环(因为如果有双向边或者环的话在压缩时就会成为一个缩点,所以不存在)
2,满足问题的点集(缩点)不能有出度,如果这个缩点V有出度,说明能到达另一个缩点U。因为V以外的所有缩点都能到达缩点V。所以U也含有到达缩点V的路径。这样的话就会形成双向边或者环。与性质1冲突。所以满足题目的缩点肯定是没有出度的。
3,满足问题要求的缩点只有1个。如果有2个或者2个以上的没有出度(该缩点里面的所有点的出度和为0),那么很显然这些点连出度都没有,自然不可能存在从一个缩点到另一个缩点的路径。也就没有任何一个满足题目的要求。
分析完这三点性质之后就已经很明了了。我们首先用tarjan算法去记录每个点属于哪个缩点组。然后再遍历每条边,更新每个组(缩点)的出度。最后再遍历缩点,记录出度为0的个数,和标号。最后统计属于该标号的缩点里的点的个数。(如果出度个数为0的值不是1的话就说明没有满足题目要求的点集)
直接上代码
#include <iostream>
#include <string.h>
using namespace std;
const int MAX_N=10001;
const int MAX_M=50001;
int ans=0,css=0,sum_css=0,n,m,answer=0;
//css记录每个点属于哪个缩点
//sum_css记录一共有多少个缩点
//answer记录满足题目要求的点的个数
int belong[MAX_N];//每个点属于哪个缩点
int head[MAX_N];
int dfn[MAX_N];//第一次出现的时间点
int low[MAX_N];//可以追寻该点或该点的子树在栈中最早出现的时间点
bool visit[MAX_N];//记录是否访问过
//手动模拟栈
int top=0;
int stack[MAX_N];
bool in_stack[MAX_N];
int outdegree[MAX_N];//记录每个缩点的出度
int zero_outdegree;//记录出度为0的缩点的编号
int number_zero=0;//记录出度为0的缩点的个数
int idx=0;
struct edge{
int u;
int to;
int next;
}eid[MAX_M];
void insert(int u,int v){
eid[ans].u=u;
eid[ans].to=v;
eid[ans].next=head[u];
head[u]=ans++;
}
int read(){
int w=1;
int s=0;
char ch=getchar();
while(!isdigit(ch)){
if(ch=='-') w*=-1;
ch=getchar();
}
while(isdigit(ch)){
s=(s<<1)+(s<<3)+ch-48;
ch=getchar();
}
return w*s;
}
void tarjan(int u){
low[u]=dfn[u]=++idx;
visit[u]=true;
stack[top++]=u;
in_stack[u]=true;
for(int i=head[u];i!=-1;i=eid[i].next){
int v=eid[i].to;
if(!visit[v]){
tarjan(v);
low[u]=min(low[u],low[v]);
}
else if(in_stack[v]){
low[u]=min(low[u],dfn[v]);
}
}
if(low[u]==dfn[u]){
css++;
sum_css++;
do{
belong[stack[--top]]=css;
in_stack[stack[top]]= false;
}while(u!=stack[top]);
}
}
void search(){
for(int i=0;i<m;++i){
int a=belong[eid[i].u];
int b=belong[eid[i].to];
if(a!=b){
outdegree[a]++;
}
}
for(int i=1;i<=sum_css;++i){
if(outdegree[i]==0){
zero_outdegree=i;
number_zero++;
}
}
if(number_zero==1){
for(int i=1;i<=n;++i){
if(belong[i]==zero_outdegree)
answer++;
}
}
}
int main() {
n=read();
m=read();
memset(head,-1, sizeof(head));
memset(dfn,0, sizeof(dfn));
memset(low,0, sizeof(low));
memset(outdegree,0, sizeof(outdegree));
for(int i=0;i<m;++i){
int u,v;
u=read();v=read();
insert(u,v);
}
for(int i=1;i<=n;++i){
if(!dfn[i])
tarjan(i);
}
search();
cout<<answer;
return 0;
}