2101 可达性统计

 1 #include <cstdio>
 2 #include <iostream>
 3 #include <algorithm>
 4 #include <cstring>
 5 #include <queue>
 6 #include <bitset>
 7 using namespace std;
 8 
 9 #define res register int
10 const int N=30000+10;
11 int head[N],ver[N<<1],nxt[N<<1],deg[N],n,m,tot(0);
12 int a[N],cnt(0);
13 bitset<N> f[N];
14 
15 inline void add(int x,int y) {
16     ver[++tot]=y; nxt[tot]=head[x]; head[x]=tot;
17 }
18 
19 queue<int> q;
20 inline void toposort()
21 {
22     for(res i=1 ; i<=n ; i++) if(!deg[i]) q.push(i);
23     while(q.size())
24     {
25         int x=q.front(); q.pop(); a[++cnt]=x;
26         for(res i=head[x] ; i ; i=nxt[i])
27         {
28             int y=ver[i];
29             if((--deg[y])==0) q.push(y);
30         }
31         
32     }
33 }
34 
35 int main()
36 {
37     scanf("%d %d",&n,&m);
38     for(res i=1 ; i<=m ; i++)
39     {
40         int x,y; scanf("%d %d",&x,&y);
41         add(x,y); deg[y]++;
42     }
43     toposort();
44     for(res i=cnt ; i>=1 ; i--)
45     {
46         int x=a[i];
47         f[x][x]=1;
48         for(res j=head[x] ; j ; j=nxt[j])
49         {
50             int y=ver[j];
51             f[x]|=f[y];
52         }
53     }
54     for(res i=1 ; i<=n ; i++) printf("%d\n",f[i].count());
55     return 0;
56 }