bzoj 3355: [Usaco2004 Jan]有序奶牛（拓扑排序+bitset）

3355: [Usaco2004 Jan]有序奶牛

Time Limit: 10 Sec   Memory Limit: 128 MB
Submit: 56   Solved: 29
[ Submit][ Status][ Discuss]

Description

    约翰的N(1≤N≤1500)头牛排成一行挤奶时，有确定的顺序．牛被编成连续的号码1．.N,他拥有L条关于奶牛顺序的信息，所有的信息都写成“A在B的前面”这样的形式，而且他知道最后一条是多余的．他觉得，有些冗余信息可以由其他信息推出，可以对这些信息进行精减．请帮助约翰删除尽可能多的冗余信息，但要保证能推出原有的顺序．可以保证的是，答案唯一，且最初的信息没有矛盾，如A在B前面，B在A前面．

Input

    第1行：两个整数N和L．

    第2到L+1行：每行两个整数X和Y(1≤X，y≤N)，表示X在Y前．无重复．

Output

    第1行：整数U.

    第2到U+I行：输出精减后的信息，每行2个数字，按第1列数字排序．

Sample Input

5 6

3 5

4 2

5 2

2 1

3 1

4 1

Sample Output

4

2 1

3 5

4 2

5 2

对于边(u, v)，如果删除这条边后u仍然能到达v，那么这条边就可以删掉

这题没有矛盾，也没有重边，也就是说是个有向无环图

考虑暴力添加每一条边，然后看这条边(u, v)是否已经联通

添加边的顺序是什么？按照边终点拓扑序从小到大排序，如果相同就按边起点拓扑序从大到小排序

如何判定是否已经连通？b[u][v]=1表示从v点能到达u点，每次加边暴力更新就好，复杂度O(nm)

其实这样过不去，不过可以bitset优化成O(nm/32)，就可以AC了

#include<stdio.h>
#include<bitset>
#include<queue>
#include<vector>
#include<algorithm>
using namespace std;
vector<int> G[1505];
int in[1505], rk[1505];
bitset<1505> b[1505];
typedef struct Res
{
	int x, y;
	bool operator < (const Res &b) const
	{
		if(rk[y]<rk[b.y] || rk[y]==rk[b.y] && rk[x]>rk[b.x])
			return 1;
		return 0;
	}
}Road;
Road s[1000005], ans[1000005];
int cnt;
bool comp(Road a, Road b)
{
	if(a.x<b.x || a.x==b.x && a.y<b.y)
		return 1;
	return 0;
}
queue<int> q;
int main(void)
{
	int n, m, x, y, i;
	scanf("%d%d", &n, &m);
	for(i=1;i<=m;i++)
	{
		scanf("%d%d", &s[i].x, &s[i].y);
		G[s[i].x].push_back(s[i].y);
		in[s[i].y]++;
	}
	for(i=1;i<=n;i++)
	{
		if(in[i]==0)
			q.push(i);
	}
	while(q.empty()==0)
	{
		x = q.front();
		q.pop();
		for(i=0;i<G[x].size();i++)
		{
			y = G[x][i];
			in[y]--;
			if(in[y]==0)
			{
				rk[y] = rk[x]+1;
				q.push(y);
			}
		}
	}
	sort(s+1, s+m+1);
	for(i=1;i<=n;i++)
		b[i][i] = 1;
	for(i=1;i<=m;i++)
	{
		if(b[s[i].y][s[i].x]==0)
		{
			b[s[i].y] |= b[s[i].x];
			ans[++cnt] = s[i];
		}
	}
	printf("%d\n", cnt);
	sort(ans+1, ans+cnt+1, comp);
	for(i=1;i<=cnt;i++)
		printf("%d %d\n", ans[i].x, ans[i].y);
	return 0;
}