强连通分量tarjan缩点——POJ2186 Popular Cows
这里的Tarjan是基于DFS,用于求有向图的强联通分量。
运用了一个点dfn时间戳和low的关系巧妙地判断出一个强联通分量,从而实现一次DFS即可求出所有的强联通分量。
dfn[u]为节点u搜索的次序编号(时间戳),即首次访问u的时间
low[u]为u或u的子树能够追溯到的最早的栈中节点的次序号,即它所在的某个强联通分量(可能是自己一个点,也可能是一棵树)对应的搜索子树的根的访问时间
不难知道,当dfn(u)=low(u)时,以u为根的搜索子树上所有节点是一个强连通分量。
伪代码如下
tarjan(u) { DFN[u]=Low[u]=++time // 为节点u设定次序编号和Low初值 Stack.push(u) // 将节点u压入栈中 for each (u, v) in E // 枚举以u为起点的每一条边 if (v is not visted) // 如果节点v未被访问过 tarjan(v) // 继续向下找 Low[u] = min(Low[u], Low[v]) else if (v in S) // 如果节点v还在栈内,则形成一个强联通分量,low值为根的dfn Low[u] = min(Low[u], DFN[v])//注意到环可能不止一个,有可能是环套环,所以我们的环是要最大的,也就是访问该子树根的时间要尽可能的早 if (DFN[u] == Low[u]) // 如果节点u是强连通分量的根 repeat v = S.pop // 将v退栈,为该强连通分量中一个顶点
print v until (u==v)
}
至于缩点,我们可以强行重新建一个图,但有时我们只需要新图的一些性质,所以有时我们可以不需要重新建图,只统计我们需要的数据就可以了。一个例子可以看下面的例题。
POJ2186 Popular Cows
题目描述
每头奶牛都梦想成为牛棚里的明星。被所有奶牛喜欢的奶牛就是一头明星奶牛。所有奶
牛都是自恋狂,每头奶牛总是喜欢自己的。奶牛之间的“喜欢”是可以传递的——如果A喜
欢B,B喜欢C,那么A也喜欢C。牛栏里共有N 头奶牛,给定一些奶牛之间的爱慕关系,请你
算出有多少头奶牛可以当明星。
输入输出格式
输入格式:
第一行:两个用空格分开的整数:N和M
第二行到第M + 1行:每行两个用空格分开的整数:A和B,表示A喜欢B
输出格式:
第一行:单独一个整数,表示明星奶牛的数量
输入输出样例
输入样例#1:
3 3
1 2
2 1
2 3
输出样例#1:
1
说明
只有 3 号奶牛可以做明星
【数据范围】
10%的数据N<=20, M<=50
30%的数据N<=1000,M<=20000
70%的数据N<=5000,M<=50000
100%的数据N<=10000,M<=50000
对于形成环的点,他们之间能相互到达,所以我们可以把它们缩成一个点。
这样我们通过tarjan缩点之后就形成了一个DAG,分析可以知道一个其他所有点都能到达的点的出度一定为0(因为缩点后不存在环,所以一旦某个点有出度,则出度所指向的点无法再返回到原来的点了),所以如果有一个出度为0的点,这个点所对应的牛的数量就是答案明星奶牛的数量,但是如果有两个以上的出度为0的点,某两个出度为0的点必定无法相互到达,故不存在题目所说的明星奶牛,为0.
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #include<algorithm> 5 const int tvt=10005; 6 const int tot=50005; 7 using namespace std; 8 int dfn[tvt],low[tvt],zhan[tvt],out[tvt],head[tot],next[tot],to[tot],belong[tvt],total[tvt],inde,n,m,len,ans,num; //out负责统计点的出度,total负责统计某个强联通分量的点的个数,inde就是时间(戳),ans统计出度为的个数 9 void add(int u,int v){ 10 num++; 11 next[num]=head[u]; 12 to[num]=v; 13 head[u]=num; 14 } 15 void tarjan(int u){ 16 dfn[u]=low[u]=++inde; 17 zhan[++len]=u; 18 for (int i=head[u];i!=0;i=next[i]){ 19 int v=to[i]; 20 if (!dfn[v]){ 21 tarjan(v); 22 low[u]=min(low[v],low[u]); 23 } 24 else low[u]=min(low[u],dfn[v]); 25 } 26 if (dfn[u]==low[u]) 27 do{int v=zhan[len--]; 28 belong[v]=u; 29 out[v]=0; 30 total[u]++; 31 }while(zhan[len+1]!=u); //即(v==u) 32 } 33 int main(){ 34 scanf("%d%d",&n,&m); 35 num=0; 36 memset(next,0,sizeof(next)); 37 memset(zhan,0,sizeof(zhan)); 38 memset(belong,0,sizeof(belong)); 39 memset(total,0,sizeof(total)); 40 memset(head,0,sizeof(head)); 41 memset(dfn,0,sizeof(dfn)); 42 memset(low,0,sizeof(low)); 43 memset(out,0,sizeof(out)); 44 memset(to,0,sizeof(to)); 45 for (int i=1;i<=m;i++){ 46 int x=0,y=0; 47 scanf("%d%d",&x,&y); 48 add(x,y);//添边 49 } 50 len=0; 51 for (int i=1;i<=n;i++) 52 belong[i]=i; 53 inde=0; 54 for (int i=1;i<=n;i++) 55 if (!dfn[i]) tarjan(i); 56 for (int i=1;i<=n;i++) 57 for (int j=head[i];j!=0;j=next[j]){ 58 int v=to[j]; 59 if (belong[i]!=belong[v]) out[belong[i]]++;//缩点,不在同一个强联通分量的话就视为两个点;把在同一个强联通分量的点的出度都统计到一个点上,达到“缩点效果” 60 } 61 ans=0; 62 for (int i=1;i<=n;i++){ 63 if ((belong[i]==i)&&(!out[i])) ans++; 64 if (ans>=2) break; 65 } 66 if (ans==1) for (int i=1;i<=n;i++) if ((belong[i]==i)&&(!out[i])){ 67 cout<<total[i]<<endl; 68 return 0; 69 } 70 else; 71 else cout<<'0'<<endl; 72 return 0; 73 }
tarjan后,我们只需要点的出度,那么我们进行“缩点”,考察缩点后点的性质,对于在同一个强联通分量的点它们之间连线产生的出度就不计,就只计不在同一个强联通分量的点的出度,并且对同一个强联通分量的所有点的出度都集中到这个分量的某个点上,这样我们就达到了我们需要的效果的“缩点”。belong数组能够告诉我们两个点是不是在同一个强联通分量里
所以,缩点实际上是根据缩点后的我们需要的性质是怎么变化的,我们就针对那个性质去统计。