做这题主要是为了学习一下tarjan的强连通分量,因为包括桥,双连通分量,强连通分量很多的求法其实都可以源于tarjan的这种方法,通过一个low,pre数组求出来。

题意:给你许多的A->B ,B->C这样的喜欢的关系,A->B ,B->C也意味着A->C,最后问你被全部别的人喜欢的cow有多少个。如果不告诉你用强连通分量,感觉可能会绕的远一些,但是如果知道了这个思路其实是很显然的。

首先是跑出每个强连通分量,在这种情况下,原来的图就变成了一棵树,一棵有有向边的树,然后不难发现,如果这棵树存在一个出度为0的点,那么它就很有可能是答案,为什么是可能呢?因为我们不知道是不是所有别的点都有路径连到它这里,所以判断的方法有这么几种吧。理论上如果只有1个出度为0的点应该就是答案了。还有一个是反向建图,然后从这个出度为0的点深搜,如果每个点都被搜到的话就说明这个点是可行的。

tarjan的算法大白书上有,具体不做过多的介绍,下面的代码有两种判断的,第一种简单,不用建图只需要统计出度,第二种复杂。在这道题了如果用Kosaraju算法会好过第二种,因为它搜出的强连通分量是按拓扑序的,tarjan则是不按的,所以用Kosaraju算法可以很容易求出那个出度为0的点。

#pragma warning(disable:4996)
#include<iostream>
#include<cstring>
#include<string>
#include<cstdio>
#include<vector>
#include<cmath>
#include<algorithm>
#include<stack>
#define maxn 10050
using namespace std; vector<int> G[maxn + 50];
int n, m; int out[maxn + 50];
int sccno[maxn + 50]; // 属于哪个强连通分量
int sccSize[maxn + 50];
int scc_cnt; // 强连通分量的数量 int pre[maxn + 50];
int low[maxn + 50];
int dfs_clock;
stack<int> S; int dfs(int u)
{
int lowu = pre[u] = ++dfs_clock;
S.push(u);
for (int i = 0; i < G[u].size(); i++){
int v = G[u][i];
if (!pre[v]){
int lowv = dfs(v);
lowu = min(lowu, lowv);
}
else if (!sccno[v]){
lowu = min(lowu, pre[v]);
}
}
if (lowu == pre[u]){
scc_cnt++;
while (1){
int x = S.top(); S.pop();
sccno[x] = scc_cnt;
sccSize[scc_cnt]++;
if (x == u) break;
}
}
return low[u] = lowu;
} void find_scc()
{
memset(sccno, 0, sizeof(sccno));
memset(pre, 0, sizeof(pre));
memset(sccSize, 0, sizeof(sccSize));
memset(low, 0, sizeof(low));
while (!S.empty()) S.pop();
dfs_clock = scc_cnt = 0;
for (int i = 1; i <= n; i++){
if (!pre[i]) dfs(i);
}
} vector<int> NG[maxn + 50]; void constructNG()
{
memset(out, 0, sizeof(out));
for (int i = 0; i <= scc_cnt; i++) NG[i].clear();
for (int i = 1; i <= n; i++){
for (int j = 0; j < G[i].size(); j++){
int u = i, v = G[i][j];
if (sccno[v] != sccno[u]){
out[sccno[u]]++;
NG[sccno[v]].push_back(sccno[u]);
}
}
}
} void ndfs(int u)
{
pre[u] = 1;
for (int i = 0; i < NG[u].size(); i++){
ndfs(NG[u][i]);
}
} int main()
{
while (cin >> n >> m)
{
for (int i = 0; i <= n; 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();
constructNG();
int foo = -1; int outnum = 0;
for (int i = 1; i <= scc_cnt; i++){
if (out[i] == 0){
foo = i; outnum++;
}
}
bool flag = true;
/*memset(pre, 0, sizeof(pre));
ndfs(foo);
bool flag = true;
for (int i = 1; i <= scc_cnt; i++){
if (!pre[i]){
flag = false; break;
}
}*/
if (outnum != 1) flag = false;
int ans = 0;
if (!flag) {
puts("0");
continue;
}
printf("%d\n", sccSize[foo]);
}
return 0;
}

POJ2186 Popular Cows 强连通分量tarjan的更多相关文章

  1. POJ2186 Popular Cows [强连通分量|缩点]

    Popular Cows Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 31241   Accepted: 12691 De ...

  2. poj 2186 Popular Cows (强连通分量+缩点)

    http://poj.org/problem?id=2186 Popular Cows Time Limit: 2000MS   Memory Limit: 65536K Total Submissi ...

  3. POJ 2186 Popular Cows --强连通分量

    题意:给定一个有向图,问有多少个点由任意顶点出发都能达到. 分析:首先,在一个有向无环图中,能被所有点达到点,出度一定是0. 先求出所有的强连通分支,然后把每个强连通分支收缩成一个点,重新建图,这样, ...

  4. poj2186 Popular Cows --- 强连通

    给一个有向图,问有多少结点是其它全部结点都能够到达的. 等价于,在一个有向无环图上,找出度为0 的结点.假设出度为0的结点仅仅有一个,那么这个就是答案.假设大于1个.则答案是0. 这题有环.所以先缩点 ...

  5. POJ 2186 Popular Cows(强连通分量缩点)

    题目链接:http://poj.org/problem?id=2186 题目意思大概是:给定N(N<=10000)个点和M(M<=50000)条有向边,求有多少个“受欢迎的点”.所谓的“受 ...

  6. POJ 2186 Popular Cows 强连通分量模板

    题意 强连通分量,找独立的块 强连通分量裸题 #include <cstdio> #include <cstdlib> #include <cstring> #in ...

  7. 强连通分量tarjan缩点——POJ2186 Popular Cows

    这里的Tarjan是基于DFS,用于求有向图的强联通分量. 运用了一个点dfn时间戳和low的关系巧妙地判断出一个强联通分量,从而实现一次DFS即可求出所有的强联通分量. §有向图中, u可达v不一定 ...

  8. 洛谷——P2341 [HAOI2006]受欢迎的牛//POJ2186:Popular Cows

    P2341 [HAOI2006]受欢迎的牛/POJ2186:Popular Cows 题目背景 本题测试数据已修复. 题目描述 每头奶牛都梦想成为牛棚里的明星.被所有奶牛喜欢的奶牛就是一头明星奶牛.所 ...

  9. 强连通分量(tarjan求强连通分量)

    双DFS方法就是正dfs扫一遍,然后将边反向dfs扫一遍.<挑战程序设计>上有说明. 双dfs代码: #include <iostream> #include <cstd ...

随机推荐

  1. android 网络_网络图片查看器

    xml <?xml version="1.0"?> -<LinearLayout tools:context=".MainActivity" ...

  2. 3月7日 Maximum Subarray

    间隔2天,继续开始写LeetCodeOj. 原题: Maximum Subarray 其实这题很早就看了,也知道怎么做,在<编程珠玑>中有提到,求最大连续子序列,其实只需要O(n)的复杂度 ...

  3. LLVM language 参考手册(译)(6)

    模块级内联汇编(Module-Level Inline Assembly) 模块包含“module-level inline assembly”块,这与GCC中的“file scope inline ...

  4. System Generator入门

      System generator 安装之后会在Simulin模块库中添加一些Xilinx FPGA专用的模块库,包括Basic Element,Communication,Control Logi ...

  5. css style与class之间的区别

    问题描述:    网页点击[导出]按钮后,将页面table内容另存成excel文件,却发现无法保存表格样式 分析过程: 1.table表格用class,而不是style.导出时并没有导出class定义 ...

  6. Atan2

    在三角函数中,两个参数的函数atan2是正切函数的 一个变种.对于任意不同时等于0的实参数x和y,atan2(y,x)所表达的意思是坐标原点为起点,指向(x,y)的射线在坐标平面上与x轴正方向之间 的 ...

  7. php取得当前访问url文件名的几种方法

    php下获取当前访问的文件名的几种方法.推荐函数:一是PHP获取当前页面的网址: dedecms用的也是这个哦. <?php //获得当前的脚本网址 function GetCurUrl() { ...

  8. android 在标题栏加上按钮

    public void onCreate(Bundle savedInstanceState) { super.onCreate(savedInstanceState); requestWindowF ...

  9. 《.NET简单企业应用》技术路线

    前言 工作三年了,一直从事基于.NET体系的企业应用开发,心得和经验也攒了点:担心时间长了给忘了,所以得给写下来,以便以后回味回味:更重要的是能让知识系统化和体系化. 本系列以一个简单的企业应用系统为 ...

  10. linux c 打印彩色字符

    #include <stdio.h> #include <string.h> int main(int argc, char **argv) { , j = , str_len ...