BZOJ 4484: [Jsoi2015]最小表示拓扑排序 bitset

最新推荐文章于 2020-10-08 16:09:49 发布

原创最新推荐文章于 2020-10-08 16:09:49 发布 · 420 阅读

0 ·

CC 4.0 BY-SA版权

—————————图论同时被 2 个专栏收录

67 篇文章

订阅专栏

—————————STL

15 篇文章

订阅专栏

4484: [Jsoi2015]最小表示

Time Limit: 20 Sec Memory Limit: 512 MB
Submit: 120 Solved: 78
[Submit][Status][Discuss]

Description

【故事背景】
还记得去年JYY所研究的强连通分量的问题吗？去年的题目里，JYY研究了对于有向图的“加边”问题。对于图论有着强烈兴趣的JYY，今年又琢磨起了“删边”的问题。
【问题描述】
对于一个N个点（每个点从1到N编号），M条边的有向图，JYY发现，如果从图中删去一些边，那么原图的连通性会发生改变；而也有一些边，删去之后图的连通性并不会发生改变。
JYY想知道，如果想要使得原图任意两点的连通性保持不变，我们最多能删掉多少条边呢？
为了简化一下大家的工作量，这次JYY保证他给定的有向图一定是一个有向无环图（JYY：大家经过去年的问题，都知道对于给任意有向图的问题，最后都能转化为有向无环图上的问题，所以今年JYY就干脆简化一下大家的工作）。

Input

输入一行包含两个正整数N和M。
接下来M行，每行包含两个1到N之间的正整数x_i和y_i，表示图中存在一条从x_i到y_i的有向边。
输入数据保证，任意两点间只会有至多一条边存在。
N<=30,000,M<=100,000

Output

输出一行包含一个整数，表示JYY最多可以删掉的边数。

Sample Input

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

Sample Output

拓扑排序后，从拓扑序大的开始更新

枚举每一个出点，按拓扑序有小到大排序（拓扑序大的有可能被小的更新掉）

依次更新

#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<complex>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<bitset>
#include<string>
#include<queue>
#include<set>
#include<map>
using namespace std;
inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
	return x*f;
}
const int N=30010;
const int M=100100;
int m,n,ecnt,last[N],fa[N],ans[M],q[N],ind[N],topo[N];
struct EDGE{int fr,to,nt,clr;}e[M];
inline void add(int u,int v)
{e[++ecnt]=(EDGE){u,v,last[u]};last[u]=ecnt;}
struct node{int key,pos;}a[N];
bitset<N>b[N];
inline bool cmp(node x,node y)
{return x.key<y.key;}
int main()
{
	n=read();m=read();int u,v;
	for(int i=1;i<=m;i++)
	{u=read();v=read();add(u,v);ind[v]++;}
	int head=1,tail=1;
	for(int i=1;i<=n;i++)if(!ind[i])q[tail++]=i;
	while(head<tail)
	{
		u=q[head++];
		for(int i=last[u];i;i=e[i].nt)
		{
			ind[e[i].to]--;
			if(!ind[e[i].to])q[tail++]=e[i].to;
		}
	}
	int ans=0;
	for(int j=tail-1;j>0;j--)
	{
		u=q[j];int cnt=0;topo[u]=j;b[u][u]=1;
		for(int i=last[u];i;i=e[i].nt)
		{
			a[++cnt]=(node){topo[e[i].to],e[i].to};
		}
		sort(a+1,a+1+cnt,cmp);
		for(int i=1;i<=cnt;i++)
		{
			if(b[u][a[i].pos])ans++;
			b[u]|=b[a[i].pos];
		}
	}
	printf("%d\n",ans);
	return 0;
}
/*
5 6
1 2
2 3
3 5
4 5
1 5
1 3

2
*/