bzoj4484 [Jsoi2015]最小表示 拓扑排序+bitset

Description

---

对于一个N个点（每个点从1到N编号），M条边的有向图，JYY发现，如果从图中删去一些边，那么原图的连通性会发生改变；而也有一些边，删去之后图的连通性并不会发生改变。
JYY想知道，如果想要使得原图任意两点的连通性保持不变，我们最多能删掉多少条边呢？
为了简化一下大家的工作量，这次JYY保证他给定的有向图一定是一个有向无环图（JYY：大家经过去年的问题，都知道对于给任意有向图的问题，最后都能转化为有向无环图上的问题，所以今年JYY就干脆简化一下大家的工作）。

输入数据保证，任意两点间只会有至多一条边存在。
N<=30,000,M<=100,000

Solution

---

如果对于一条边(x,y)存在z使得x能到z，z能到y，那么这条边显然可以不要
512M内存提示我们两个bitset[x]求出x能到达的点和能到达x的点，&一下统计就没了

Code

---

#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <vector>
#include <bitset>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define fi first
#define se second

typedef std:: pair <int,int> pair;
const int N=30005;

std:: bitset <N> bit1[N],bit2[N];
std:: vector <int> vec[N],rec[N];
std:: vector <pair> edge;

int queue[N],d1[N],d2[N];

int read() {
	int x=0,v=1; char ch=getchar();
	for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
	for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());
	return x*v;
}

void add_edge(int x,int y) {
	vec[x].push_back(y); d1[y]++;
	rec[y].push_back(x); d2[x]++;
}

void top_sort(int n) {
	int head=1,tail=0;
	rep(i,1,n) {
		if (!d1[i]) queue[++tail]=i;
		bit1[i][i]=1;
	}
	for (;head<=tail;) {
		int now=queue[head++];
		for (int i=0;i<vec[now].size();++i) {
			bit1[vec[now][i]]|=bit1[now];
			if (!(--d1[vec[now][i]])) {
				queue[++tail]=vec[now][i];
			}
		}
	}
	head=1,tail=0;
	rep(i,1,n) {
		if (!d2[i]) queue[++tail]=i;
		bit2[i][i]=1;
	}
	for (;head<=tail;) {
		int now=queue[head++];
		for (int i=0;i<rec[now].size();++i) {
			bit2[rec[now][i]]|=bit2[now];
			if (!(--d2[rec[now][i]])) {
				queue[++tail]=rec[now][i];
			}
		}
	}
}

int main(void) {
	int n=read(),m=read();
	rep(i,1,m) {
		int x=read(),y=read();
		add_edge(x,y);
		edge.push_back(pair(x,y));
	}
	top_sort(n);
	int ans=0;
	rep(i,0,m-1) {
		std:: bitset <N> wjp;
		wjp=bit1[edge[i].se]&bit2[edge[i].fi];
		wjp[edge[i].se]=wjp[edge[i].fi]=0;
		if (wjp.count()) ans++;
	}
	printf("%d\n", ans);
	return 0;
}