poj2186Popular Cows【强连通】【debug两小时】

这个题我调试两个小时，快调哭了都

写一下心得

大意：告诉你一群牛  有n头 （n<10^4)

然后告诉你m对羡慕关系  如  a牛羡慕b牛

并且羡慕可以传递  如a羡慕b  b羡慕c  则a羡慕c

现在问有多少只牛是被所有的牛羡慕的

分析：我的第一个思路是遍历每个牛 然后从每个牛建立反向边  然后dfs看能否扫到所有的牛  但是这样时间复杂度O(n*(n + e) )  承受不住

后来想到可以像flody传递闭包那样但是同样时间复杂度承受

后来想到这个可以强连通一下子

我们首先想这么一个问题
对于一个有向无环图

满足一个牛被所有牛都羡慕的条件是  ：  1  图连通  2  该牛的出度为零（因为他是被所有牛都羡慕的那只  然后如果他再羡慕其他的牛  就会形成环）

对于本题是有环的   用强连通  进行缩点  变成一个有向无环图  再根据那俩条件进行判定即可

但是写的时候  首先是我把图中有孤立点这种情况忽略了

然后就是  脑残的break之后……想起来就心痛啊

代码：

  1 #include <iostream>
  2 #include <cstdio>
  3 #include <cstring>
  4 using namespace std;
  5
  6 const int maxn = 50005;
  7
  8 int n, m;
  9 int tot, cnt;
 10 struct Node {
 11     int to, next;
 12 }e[maxn];
 13 int head[maxn];
 14
 15 void add(int u, int v) {
 16     e[tot].to = v;
 17     e[tot].next = head[u];
 18     head[u] = tot++;
 19 }
 20 int s[maxn];
 21 int ins[maxn];
 22 int sc;
 23 int ans;
 24 int num[maxn];
 25
 26 int low[maxn], dfn[maxn];
 27
 28 void init() {
 29     tot = 1;
 30     memset(head, 0, sizeof(head));
 31     memset(low, 0, sizeof(low));
 32     memset(dfn, 0, sizeof(dfn));
 33     memset(ins, 0, sizeof(ins));
 34     memset(num, 0, sizeof(num));
 35     cnt = 1;
 36     sc = 0;
 37     ans = 0;
 38 }
 39
 40 void dfs(int u) {
 41     s[sc++] = u;
 42     ins[u] = 1;
 43     for(int i = head[u]; i; i = e[i].next) {
 44         int v = e[i].to;
 45         if(!dfn[v]) {
 46             dfn[v] = low[v] = cnt++;
 47             dfs(v);
 48             low[u] = min(low[u], low[v]);
 49         } else if(ins[v] && !num[v]) {
 50             low[u] = min(low[u], dfn[v]);
 51         }
 52     }
 53     if(low[u] == dfn[u]) {
 54         ans++;
 55         int x;
 56         do {
 57             x = s[--sc];
 58             num[x] = ans;
 59             ins[x] = 0;
 60         } while(x != u);
 61     }
 62
 63 }
 64 int outed[maxn];
 65
 66 int main() {
 67     int u, v;
 68     while(EOF != scanf("%d %d",&n, &m) ) {
 69         init();
 70         for(int i = 1; i <= m; i++) {
 71             scanf("%d %d",&u, &v);
 72             add(u, v);
 73         }
 74         for(int i = 1; i <= n; i++) {
 75             if(!dfn[i]) {
 76                 dfn[i] = low[i] = cnt++;
 77                 dfs(i);
 78             }
 79         }
 80     //    for(int i = 1; i <= n; i++) {
 81     //        printf("%d ", num[i]);
 82     //    } puts("");
 83         int ans_num = -1;
 84         memset(outed, 0, sizeof(outed));
 85         for(int i = 1; i <= n; i++) {
 86 //            printf("i = %d\n", i);
 87             if(num[i] == 0) {
 88                 ans_num = 0;
 89                 break;
 90             }
 91             for(int j = head[i]; j ; j = e[j].next) {
 92                 int x = e[j].to;
 93                 if(num[i] != num[x]) {
 94                     outed[num[i]]++;
 95                 }
 96             }
 97         }
 98     //    for(int i = 1; i <= n; i++) {
 99     //        printf("*%d ", outed[i]);
100     //    } puts("");
101     //    printf("ans = %d\n", ans);
102         int flag = 0;
103         if(ans_num != 0) {
104             ans_num = 0;
105             int ax = 0;
106             for(int i = 1; i <= ans; i++) {
107                 if(outed[i] == 0) {
108                     flag++;
109                     if(flag >= 2) {
110                         ans_num = 0;
111                         break;
112                     }
113                     ax = i;
114                 }
115             }
116             if(flag < 2) {
117                 for(int i = 1; i <= n; i++) {
118                     if(num[i] == ax) {
119                         ans_num++;
120                     }
121                 }
122             }
123         }
124         printf("%d\n", ans_num);
125     }
126     return 0;
127 }