POJ 2186-Popular Cows (图论-强联通分量Korasaju算法)

时间:2023-03-09 14:31:37
POJ 2186-Popular Cows (图论-强联通分量Korasaju算法)

题目链接:http://poj.org/problem?id=2186

题目大意:有n头牛和m对关系, 每一对关系有两个数(a, b)代表a牛认为b牛是“受欢迎”的,且这种关系具有传递性, 如果a牛认为b牛“受欢迎”, b牛认为c牛“受欢迎”, 那么a牛也认为c牛“受欢迎”。 现在想知道有多少头牛受除他本身外其他所有牛的欢迎?

解题思路:如果有两头或者多头牛受除他本身外其他所有牛的欢迎, 那么在这两头或者多头牛之中, 任意一头牛也受两头或者多头牛中别的牛的欢迎, 即这两头或者多头牛同属于一个强联通分量, 且其他强联通分量都可以到达该强联通分量。那么可以用Korasaju算法进行强联通分量的分解, 然后还能够得到各个强联通分量拓扑排序后的顺序, 那么唯一可以成为解的只有拓扑序最后的那个联通分量, 并且需要检查其他强联通分量是否能全部到达这个强联通分量, 如果能够全部可达, 那么该强联通分量有多少元素即为问题的解, 否则为0;

强联通分量的分解的步骤:

1:对于图G, 深度优先遍历G, 算出每个节点u结束的时间s[u], 起点如何选择无所谓。

2:对于图G的转置图rG, 选择遍历的起点时, 按照节点的结束时间从大到小进行, 遍历过程中, 一遍遍历, 一遍给节点做分类标记, 每找到一个新的起点, 分类标记加一

3:对于第二步骤中,每一个相同标记的点即为一个强联通分量

代码如下:

#include<stdio.h>
#include<string.h>
#include<vector>
#include<algorithm>
using namespace std;
const int N = ; vector<int>G[N], rG[N], vs;
bool used[N];
int v;
int cmp[N]; void add_edge(int a, int b)
{
G[a].push_back(b);
rG[b].push_back(a);
} void dfs(int n)
{
used[n] = true;
for(int i=; i<G[n].size(); ++ i)
{
int v = G[n][i];
if(used[v] == false)
dfs(v);
}
vs.push_back(n);
} void rdfs(int n, int k)
{
used[n] = true, cmp[n] = k;
for(int i=; i<rG[n].size(); ++ i)
{
int v=rG[n][i];
if(used[v] == false)
rdfs(v, k);
}
} int scc()
{
memset(used, false, sizeof(used));
vs.clear();
for(int i=; i<v; ++ i)
if(used[i] == false)
dfs(i); memset(used, false, sizeof(used));
int k=;
for(int i=vs.size()-; i>=; -- i)
{
if(used[vs[i]] == false)
rdfs(vs[i], k++);
}
return k;
} int main()
{
int m;
scanf("%d%d", &v, &m);
for(int i=; i<=m; ++ i)
{
int a, b;
scanf("%d%d", &a, &b);
add_edge(a-, b-);
} int n = scc(); int u = , num = ;
for(int i=; i<v; ++ i)
if(cmp[i] == n-)
{
u = i, num ++;
}
memset(used, false, sizeof(used));
rdfs(u, );
for(int i=; i<v; ++ i)
if(used[i] == false)
{
num = ;
break;
}
printf("%d\n", num);
}

注:新学的强联通分量,那里不对请指出, 万分感谢!