题意
给出一张DAG,要求删除尽量多的边使得连通性不变.(即:若删边前u到v有路径,则删边后仍有路径).点数30000,边数100000.
分析
如果从u到v有(u,v)这条边,且从u到v只有这一条路径,那么这条边必须保留.否则这条边一定可以删除.因为如果有不止一条路径从u到v,必然存在点x(x!=u,x!=v)使得u可到达x,x可到达v.而删边后必然也满足u可到达x,x可到达v,所以直接删掉(u,v)这条边就可以了.
刚才的分析已经给出了一个判定方法.既然如果有不止一条路径从u到v,必然存在点x(x!=u,x!=v)使得u可到达x,x可到达v,那么我们对每条边(u,v),枚举是否存在这样的x即可.这需要我们求出每个点能到达的点的集合,以及能到达这个点的集合.大力压位一波就好了.因为是DAG所以这个集合可以递推.复杂度O(nm/32).其实这题是看内存猜算法系列,榜上清一色的120多兆,不是压位还能是啥我是200多兆
#include<cstdio>
const int mod=1000000007;
const int maxn=30005,maxm=200005;
struct edge{int to,next;
}lst[maxm],lst2[maxm];int len=1,first[maxn],len2=1,first2[maxn];
void addedge(int a,int b){lst[len].to=b;lst[len].next=first[a];first[a]=len++;
}
void addedge2(int a,int b){lst2[len2].to=b;lst2[len2].next=first2[a];first2[a]=len2++;
}
int sz;
int reach[maxn][maxn/32+2],from[maxn][maxn/32+2];
int getbit(int u,int x){return (reach[u][x/32]>>(x&31))&1;
}
void revbit(int u,int x){reach[u][x/32]^=(1<<(x&31));
}
void revbit2(int u,int x){from[u][x/32]^=(1<<(x&31));
}
bool vis[maxn];
void dfs(int x){if(vis[x])return;vis[x]=true;for(int pt=first[x];pt;pt=lst[pt].next){dfs(lst[pt].to);for(int i=0;i<sz;++i)reach[x][i]|=reach[lst[pt].to][i];}revbit(x,x);
}
void dfs2(int x){if(vis[x])return;vis[x]=true;for(int pt=first2[x];pt;pt=lst2[pt].next){dfs2(lst2[pt].to);for(int i=0;i<sz;++i)from[x][i]|=from[lst2[pt].to][i];}revbit2(x,x);
}
int main(){int n,m;scanf("%d%d",&n,&m);sz=(n+31)/32+1;for(int i=1,a,b;i<=m;++i){scanf("%d%d",&a,&b);addedge(a,b);addedge2(b,a);}for(int i=1;i<=n;++i)if(!vis[i])dfs(i);for(int i=1;i<=n;++i)vis[i]=0;for(int i=1;i<=n;++i)if(!vis[i])dfs2(i);for(int i=1;i<=n;++i)revbit(i,i),revbit2(i,i);int ans=0;for(int i=1;i<=n;++i){for(int pt=first[i];pt;pt=lst[pt].next){int y=lst[pt].to;for(int j=0;j<sz;++j){if(from[y][j]&reach[i][j]){ans++;break;}}}}printf("%d\n",ans);return 0;
}