JZOJ.5328【NOIP2017模拟8.22】世界线

Description

 

Input

Output

 

Sample Input

5 5
1 2
1 3
2 3
3 4
4 5

Sample Output

5

 

Data Constraint

 

Hint

样例解释

先拓扑排序，从出度为0的点开始搜，连这个点的前一个点以及之前的点都要和这个点连边，这样算出每个点能到达的其他点再减去已有的边就可以了。

然而需要用bitset优化状态，然而还有可能会爆空间，我们需要分几次搜索.....

#include <iostream>
#include <cstring>
#include <cstdio>
#include <bitset>
#define N 1300000
using namespace std;
bitset<1001>qwq[N];
int du[N/2],n,m,head[N],next[N],to[N],num,q[N/2],t;
long long ans;
void add(int u,int v){
    num++;
    next[num]=head[u];
    to[num]=v;
    head[u]=num;
}
int main(){
    freopen("worldline.in","r",stdin);
    freopen("worldline.out","w",stdout);
    scanf("%d%d",&n,&m);
    ans=-m;
    num=0;
    for (int i=1,u=0,v=0;i<=m;i++){
        scanf("%d%d",&u,&v);
        add(u,v);
        du[v]++;
    }
    t=0;
    for (int i=1;i<=n;i++)
        if (!du[i]) q[++t]=i;
    for (int v=0,i=1;i<=t;i++)
        for (int j=head[q[i]];j;j=next[j]){
            v=to[j];
            du[v]--;
            if (!du[v]) q[++t]=v;
        }
    for (int l=1,r=1000;l<=n;l=r+1,r+=1000){
        for (int i=1;i<=n;i++) qwq[i].reset();
        for (int i=t;i>=1;i--){
            for (int j=head[q[i]];j;j=next[j]){
                qwq[q[i]]|=qwq[to[j]];
            }
            ans+=qwq[q[i]].count();
            if ((q[i]>=l)&&q[i]<=r) qwq[q[i]][q[i]-l]=1; 
        }
    }
    printf("%lld\n",ans);
}