AcWing P164 可达性统计 题解

Analysis

这道题我一开始想到的是传递闭包，但是时间复杂度是n³，也开不下30000*30000的数组，所以我想到了拓扑+状态压缩（bitset），从后往前找，把能到达的点能到哪里用位运算赋到上一个中，最后调用.count()输出就好了。

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<queue>
 #include<bitset>
 #define maxn 30010
 using namespace std;
 inline int read()
 {
     int x=;
     bool f=;
     char c=getchar();
     for(; !isdigit(c); c=getchar()) if(c=='-') f=;
     for(; isdigit(c); c=getchar()) x=(x<<)+(x<<)+c-'';
     if(f) return x;
     return -x;
 }
 inline void write(int x)
 {
     if(x<){putchar('-');x=-x;}
     if(x>)write(x/);
     putchar(x%+'');
 }
 queue<int> q;
 bitset<maxn> s[maxn];
 int n,m,cnt,num;
 int head[*maxn],rd[maxn],path[maxn];
 struct node
 {
     int v,nex;
 }edge[*maxn];
 inline void add(int u,int v)
 {
     cnt++;
     edge[cnt].v=v;
     edge[cnt].nex=head[u];
     head[u]=cnt;
 }
 inline void topo()
 {
     for(int i=;i<=n;i++)
     {
         if(rd[i]==)q.push(i);
     }
     while(!q.empty())
     {
         int t=q.front();
         q.pop();
         num++;
         path[num]=t;
         for(int i=head[t];i;i=edge[i].nex)
         {
             int to=edge[i].v;
             rd[to]--;
             if(rd[to]==)q.push(to);
         }
     }
 }
 int main()
 {
     n=read();m=read();
     for(int i=;i<=m;i++)
     {
         int x,y;
         x=read();y=read();
         add(x,y);
         rd[y]++;
     }
     topo();
     for(int i=num;i>=;i--)
     {
         int x=path[i];
         s[x][x]=;
         for(int j=head[x];j;j=edge[j].nex)
         {
             int to=edge[j].v;
             s[x]|=s[to];
         }
     }
     for(int i=;i<=n;i++)
     {
         write(s[i].count());
         printf("\n");
     }
     return ;
 }

请各位大佬斧正（反正我不认识斧正是什么意思）