BZOJ3355
3355: [Usaco2004 Jan]有序奶牛
Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 37 Solved: 19
[Submit][Status][Discuss]
Description
Input
Output
Sample Input
3 5
4 2
5 2
2 1
3 1
4 1
Sample Output
2 1
3 5
4 2
5 2
HINT
3在1前,4在1前可推.输出的每一行,不能被其他推出.
分析:对于一条边x->y,若去掉之后x不能到达y,那么它是必需的。我们首先需要求出拓扑序列,然后我们按照拓扑序列反向加边,这样的目的就是为了在已知顺序的情况下维护哪些边可以互相到达,然后用一个bieset去判断一下哪些边是必须的。感谢claris和付大神对我的帮助
#include "iostream"
#include "cstdio"
#include "cstring"
#include "bitset"
#include "algorithm"
using namespace std;
const int maxn=+;
const int maxm=+;
int n,m,ed,ed2;
int ind[maxn]; //入度
int adj[maxm]; //栈顶边对应的顶点
int next1[maxm]; //上一条边的编号
int tail[maxm]; //栈顶边的编号
int vis[maxn]; //存储拓扑排序序列
bitset<maxn>f[maxn];
void add(int x,int y){
ind[y]++;
adj[++ed]=y;
next1[ed]=tail[x];
tail[x]=ed;
}
int adj2[maxm];
int next2[maxm];
int tail2[maxm];
void add2(int x,int y){
adj2[++ed2]=y;
next2[ed2]=tail2[x];
tail2[x]=ed2;
}
typedef pair<int,int> P;
P p[maxm];
int main()
{
//freopen("a.txt", "r", stdin);
//freopen("b.txt", "w", stdout);
scanf("%d%d",&n,&m);
for(int i=;i<=n;i++) f[i][i]=;
ed=;
while(m--){ //数组式邻接表建图
int x,y;
scanf("%d%d",&x,&y);
add(x,y);
} //构造拓扑排序序列
int cnt=;
for(int i=;i<=n;i++)
if(!(ind[i]))
vis[++cnt]=i; int h=;
ed2=;
while(h<=cnt){
int a;
for(int i=tail[a=vis[h++]];i;add2(adj[i],a),i=next1[i])
if(!(--ind[adj[i]]))
vis[++cnt]=adj[i];
} int ans=;
for(int i=;i<=n;i++){
int a;
for(int j=tail2[a=vis[i]];j;f[a]|=f[adj2[j]],j=next2[j])
if(!f[a][adj2[j]])
p[++ans]=P(adj2[j],a);
}
sort(p+,p++ans);
printf("%d\n",ans);
for(int i=;i<=ans;i++)
printf("%d %d\n",p[i].first,p[i].second);
return ;
}
BZOJ3355的更多相关文章
- BZOJ3355 : [Usaco2004 Jan]有序奶牛
对于一条边x->y,若去掉之后x不能到达y,那么它是必需的. 首先拓扑排序求出拓扑序,然后按照终点拓扑序为第一关键字,起点拓扑序为第二关键字从小到大加边. 对于每个点,维护一个bitset,表示 ...
随机推荐
- 深入理解iPhone数据持久化(手把手教你iphone开发 – 基础篇)
在所有的移动开发平台数据持久化都是很重要的部分:在j2me中是rms或保存在应用程序的目录中,在symbian中可以保存在相应的磁盘目录中和数据库中.symbian中因为权限认证的原因,在3rd上大多 ...
- android 按两次物理返回键退出程序
<?xml version="1.0" encoding="utf-8"?> <!-- 定义当前布局的基本LinearLayout --> ...
- jinjia2模板学习
http://docs.jinkan.org/docs/jinja2/templates.html#
- 一起学习CMake – 01
一起学习CMake – 01 本节介绍CMake里最常用的三个命令,分别是cmake_minimum_required; project; add_executable等. CMake是个好东西,在使 ...
- weexapp 开发流程(三)其他页面创建
1.首页 (1)轮播图 步骤一:创建 轮播图 组件(Slider.vue) src / assets / components / Slider.vue <!-- 轮播图 组件 --> & ...
- plsql 无需配置客户端连接.
plsql 无需配置客户端连接. (DESCRIPTION=(ADDRESS_LIST=(ADDRESS=(PROTOCOL=TCP)(HOST=192.168.0.5)(PORT=1521)))(C ...
- Django-权限信息中间件操作
# 在当前app下新建一个middleware的文件夹,然后就可以尽情的写中间件了,只能是这个名字,切记~@!import re from django.shortcuts import redire ...
- 01背包+卡精度 Hdu 2955
<span style="color:#3333ff;">/* ---------------------------------------------------- ...
- const& 的东西
class_name ( class_name const & source ); 是拷贝构造函数的标准声明. 它和如下声明是一个意思 class_name ( const class_nam ...
- 2016/07/05 zend optimizer
Zend Optimizer是由PHP核心引擎“Zend” http://www.zend.com 创建者Zend技术公司所开的免费PHP优化软件.据Zend公司透露使用这个软件某些情况下至少可以提高 ...