【JZOJ5328】【NOIP2017提高A组模拟8.22】世界线

Description

Data Constraint

Solution

这题用bitset就能过，直接O($N^2\over32$)草过。但由于空间会爆炸，我们可以每10000位做一次。

Code

#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<set>
#include<bitset>
using namespace std;
const int maxn=1e5+5;
int first[maxn],last[maxn],next[maxn],d[maxn],v[maxn],p[maxn],a[maxn];
int n,m,i,t,j,k,l,x,y,z,num,ans;
bitset<3001>f[maxn];
void lian(int x,int y){
    last[++num]=y;next[num]=first[x];first[x]=num;
}
void bfs(int l,int r){
    int i,j=0;
    for (i=1;i<=n;i++){
        d[i]=p[i];
        if (!d[i]) v[++j]=i;
    }
    i=0;
    while (i<j){
        x=v[++i];
        for (t=first[x];t;t=next[t]){
            d[last[t]]--;
            if (!d[last[t]]) v[++j]=last[t];
        }
    }
    for (i=n;i>=1;i--){
        x=v[i];f[x].reset();
        for (t=first[x];t;t=next[t]){
            y=last[t];
            if (y>l && y<=r){
                f[x].set(y-l);
            }
            f[x]=f[x]|f[y];
        }
    }
    for (i=1;i<=n;i++)
        ans+=f[i].count();
}
void dg(int x){
    int t;d[x]=0;
    for (t=first[x];t;t=next[t]){
        dg(last[t]);
        d[x]+=d[last[t]]+1;
    }
}
int main(){
    freopen("worldline.in","r",stdin);freopen("worldline.out","w",stdout);
    scanf("%d%d",&n,&m);
    for (i=1;i<=m;i++)
        scanf("%d%d",&x,&y),lian(x,y),d[y]++,p[y]++;
    for (i=1;i<=n/3000+1;i++)
        bfs((i-1)*3000,i*3000);
    ans-=m;
    printf("%d\n",ans);
}