【bieset优化】可达性统计

lamentropetion

于 2023-03-15 23:49:43 发布

阅读量66

点赞数

分类专栏：图论文章标签：算法 Powered by 金山文档

本文链接：https://blog.csdn.net/weixin_62528401/article/details/129573767

版权

图论专栏收录该内容

87 篇文章 0 订阅

订阅专栏

文章介绍了一种利用拓扑排序和bitset进行可达性统计的方法，通过构建图的边并进行拓扑排序，然后用bitset快速计算1的个数，以此解决AcWing题库中的相关问题。

摘要由CSDN通过智能技术生成

还没怎么用过bitset

感觉就是可以快死位运算和统计1的个数的01数组

164. 可达性统计 - AcWing题库

题意：

思路：

拓扑排序把拓扑序跑出来

跑出来之后逆着拓扑序把可达的点的集合都或上，然后可以快速求出1的个数

Code：

#include <bits/stdc++.h>
using namespace std;
const int mxn=3e4+10,mxe=3e4+10;
struct ty{
    int to,next;
}edge[mxe<<1];
vector<int> v;
queue<int> q;
int n,m,x,y,tot=0;
bitset<mxn> f[mxn];
int head[mxn],in[mxn],dp[mxn];
void init(){
    tot=0;
    for(int i=0;i<=n;i++){
        head[i]=-1;
        in[i]=0;
        dp[i]=1;
    }
    while(!q.empty()) q.pop();
}
void add(int u,int v){
    edge[tot].to=v;
    edge[tot].next=head[u];
    head[u]=tot++;
}
void topsort(){
    for(int i=1;i<=n;i++){
        if(!in[i]) q.push(i);
    }
    while(!q.empty()){
        int u=q.front();
        q.pop();
        v.push_back(u);
        for(int i=head[u];~i;i=edge[i].next){
            in[edge[i].to]--;
            if(!in[edge[i].to]) q.push(edge[i].to);
        }
    }
}
int main(){
    ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
    cin>>n>>m;
    init();
    for(int i=1;i<=m;i++){
        cin>>x>>y;
        add(x,y),in[y]++;
    }
    topsort();
    for(int i=v.size()-1;i>=0;i--){
        int u=v[i];
        f[u][u]=1;
        for(int j=head[u];~j;j=edge[j].next){
            f[u]|=f[edge[j].to];
        }
    }
    for(int i=1;i<=n;i++) cout<<f[i].count()<<'\n';
}