ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

CH2101 可达性统计

2019-02-22 18:43:17  阅读:244  来源: 互联网

标签:ch int tot rg CH2101 data 可达性 统计 getchar


题意

给定一张N个点M条边的有向无环图,分别统计从每个点出发能够到达的点的数量。N,M≤30000。

分析

有向无环图,可以按拓扑序逆序统计答案。可以用bitset维护可达性。

时间复杂度\(O(N (N+M)/32 )\),空间大小\(N^2/8\)字节。

代码

#include<bits/stdc++.h>
#define rg register
#define il inline
#define co const
template<class T>il T read(){
    rg T data=0,w=1;
    rg char ch=getchar();
    while(!isdigit(ch)){
        if(ch=='-') w=-1;
        ch=getchar();
    }
    while(isdigit(ch))
        data=data*10+ch-'0',ch=getchar();
    return data*w;
}
template<class T>il T read(rg T&x){
    return x=read<T>();
}
typedef long long LL;
using namespace std;

co int N=3e4+1;
int n,m,deg[N],a[N];
int ver[N],Next[N],head[N],tot,cnt;
bitset<N> f[N];
void add(int x,int y){
    ver[++tot]=y,Next[tot]=head[x],head[x]=tot;
    ++deg[y];
}
void topsort(){
    queue<int> q;
    for(int i=1;i<=n;++i)
        if(deg[i]==0) q.push(i);
    while(q.size()){
        int x=q.front();q.pop();
        a[++cnt]=x;
        for(int i=head[x];i;i=Next[i]){
            int y=ver[i];
            if(--deg[y]==0) q.push(y);
        }
    }
}
void calc(){
    for(int i=cnt;i;--i){
        int x=a[i];
        f[x][x]=1;
        for(int i=head[x];i;i=Next[i]){
            int y=ver[i];
            f[x]|=f[y];
        }
    }
}
int main(){
//  freopen(".in","r",stdin);
//  freopen(".out","w",stdout);
    read(n),read(m);
    for(int i=1,x,y;i<=m;++i){
        read(x),read(y);
        add(x,y);
    }
    topsort();
    calc();
    for(int i=1;i<=n;++i) printf("%llu\n",f[i].count());
    return 0;
}

标签:ch,int,tot,rg,CH2101,data,可达性,统计,getchar
来源: https://www.cnblogs.com/autoint/p/10419979.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有