[POI2014]Rally
OJ题号:
BZOJ3832、洛谷3573
思路:
建立超级源汇$S$和$T$,DP求出分别以$S$和$T$为源点的最长路$diss$和$dist$。
对于每条边$i$,设定一个权值$w_i=diss_{i.from}+dist_{i.to}-1$。
表示原图中包含这条边的从$S$到$T$的最长路。
然后按照拓扑序删点$x$,用堆维护不包含$x$的最长路长度。
然而一次性不能把所有边放进去,不然会MLE一个点(因为这个调了一个晚上)。
应该在换$x$的时候,把老$x$的出边重新加入,并将新$x$的入边删去。
注意开的数组不能太多,能合并的信息尽量合并,(比如所有边正反边用一个数组存,取值的时候用异或),不然把堆修改以后还是会MLE。
#include<queue>
#include<cstdio>
#include<cctype>
#include<vector>
#include<ext/pb_ds/priority_queue.hpp>
inline int getint() {
register char ch;
while(!isdigit(ch=getchar()));
register int x=ch^'';
while(isdigit(ch=getchar())) x=(((x<<)+x)<<)+(ch^'');
return x;
}
const int inf=0x7fffffff;
const int V=,E=;
struct Edge {
int from,to;
};
Edge e[E];
int w[E];
int s,t;
int n,m;
std::vector<int> eids[V],eidt[V];
int inds[V]={},indt[V]={};
inline void add_edge(const int u,const int v,int *ind,std::vector<int> *eid,const int i) {
eid[u].push_back(i);
ind[v]++;
}
int diss[V]={},dist[V]={};
std::queue<int> top;
inline void Kahn(const int s,std::vector<int> *eid,int *dis,int *ind,const int op=) {
std::queue<int> q;
q.push(s);
while(!q.empty()) {
int x=q.front();
q.pop();
if(op) top.push(x);
for(register unsigned i=;i<eid[x].size();i++) {
int y=e[eid[x][i]].from^e[eid[x][i]].to^x;
dis[y]=std::max(dis[y],dis[x]+);
if(!--ind[y]) q.push(y);
}
}
}
__gnu_pbds::priority_queue<int> q;
__gnu_pbds::priority_queue<int>::point_iterator p[E];
int v,ans=inf;
int cnt=;
inline void solve() {
while(!top.empty()) {
int x=top.front();
top.pop();
for(register unsigned i=;i<eidt[x].size();i++) {
q.erase(p[eidt[x][i]]);
}
if(!q.empty()) {
if((x!=s)&&(x!=t)&&(q.top()<ans)) {
ans=q.top();
v=x;
}
}
for(register unsigned i=;i<eids[x].size();i++) {
p[eids[x][i]]=q.push(w[eids[x][i]]);
}
}
}
int main() {
n=getint(),m=getint();
s=,t=n+;
for(register int i=;i<=n;i++) {
e[cnt].from=s,e[cnt].to=i;
add_edge(s,i,inds,eids,cnt);
add_edge(i,s,indt,eidt,cnt);
cnt++;
}
for(register int i=;i<m;i++) {
int &u=e[cnt].from=getint(),&v=e[cnt].to=getint();
add_edge(u,v,inds,eids,cnt);
add_edge(v,u,indt,eidt,cnt);
cnt++;
}
for(register int i=;i<=n;i++) {
e[cnt].from=i,e[cnt].to=t;
add_edge(i,t,inds,eids,cnt);
add_edge(t,i,indt,eidt,cnt);
cnt++;
}
Kahn(s,eids,diss,inds,);
Kahn(t,eidt,dist,indt);
for(register int i=;i<cnt;i++) {
w[i]=diss[e[i].from]+dist[e[i].to]-;
}
solve();
printf("%d %d\n",v,ans);
return ;
}
[POI2014]Rally的更多相关文章
- 【BZOJ3832】[POI2014]Rally(拓扑排序,动态规划)
[BZOJ3832][POI2014]Rally(拓扑排序,动态规划) 题面 BZOJ,权限题 洛谷 题解 这题好强啊,感觉学了好多东西似的. 首先发现了一个图画的很好的博客,戳这里 然后我来补充一下 ...
- 3832: [Poi2014]Rally
3832: [Poi2014]Rally 链接 分析: 首先可以考虑删除掉一个点后,计算最长路. 设$f[i]$表示从起点到i的最长路,$g[i]$表示从i出发到终点的最长路.那么经过一条边的最长路就 ...
- BZOJ3832[Poi2014]Rally——权值线段树+拓扑排序
题目描述 An annual bicycle rally will soon begin in Byteburg. The bikers of Byteburg are natural long di ...
- 【BZOJ】3832: [Poi2014]Rally
题意 \(n(2 \le n \le 500000)\)个点\(m(1 \le m \le 1000000)\)条边的有向无环图,找到一个点,使得删掉这个点后剩余图中的最长路径最短. 分析 神题不会做 ...
- BZOJ 3832: [Poi2014]Rally
Sol 线段树+拓扑序. 先把图的拓扑序搞出来,然后统计从起点到该点最长链,从该点到终点的最长链,然后建个起点终点,这里跟网络流很像,把它统一到一个有起点的图中,这里也要注意下细节处理.S,T的一个边 ...
- BZOJ3832 : [Poi2014]Rally
f[0][i]为i出发的最长路,f[1][i]为到i的最长路 新建源汇S,T,S向每个点连边,每个点向T连边 将所有点划分为两个集合S与T,一开始S中只有S,其它点都在T中 用一棵线段树维护所有连接属 ...
- BZOJ3832: [Poi2014]Rally(拓扑排序 堆)
题意 题目链接 Sol 最直观的思路是求出删除每个点后的最长路,我们考虑这玩意儿怎么求 设\(f[i]\)表示以\(i\)结尾的最长路长度,\(g[i]\)表示以\(i\)开始的最长路长度 根据DAG ...
- BZOJ:3832: [Poi2014]Rally
题意: 给出$DAG$,询问删掉哪个点之后最长路径最短 思路: 我们令$f[x]$表示从最远的点到达它的距离,$g[x]$表示它能够到达最远的点的距离 那么对于$(x -> y)$一条边来说,它 ...
- BZOJ3832 [Poi2014]Rally 【拓扑序 + 堆】
题目链接 BZOJ3832 题解 神思路orz,根本不会做 设\(f[i]\)为到\(i\)的最长路,\(g[i]\)为\(i\)出发的最长路,二者可以拓扑序后\(dp\)求得 那么一条边\((u,v ...
随机推荐
- jdk1.8.0_45源码解读——HashMap的实现
jdk1.8.0_45源码解读——HashMap的实现 一.HashMap概述 HashMap是基于哈希表的Map接口实现的,此实现提供所有可选的映射操作.存储的是<key,value>对 ...
- html5 canvas 弧形描边渐变
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/ ...
- unp的hello word
最近无聊,再学学unp吧.学学网络编程.下面是一个demo: intro/daytimetcpcli.c #include "unp.h" int main(int argc, c ...
- Asp.Net中索引器的用法
索引器定义类似于属性,但其功能与属性并不相同.索引器提供一种特殊的方法编写get和set访问器.属性可以像访问字段一样访问对象的数据,索引器可以使用户像访问数组一样访问类成员. 一.索引器特性 1.g ...
- Grunt 插件发布过程;
首先: Npm adduser User: xxxx Password: yyyy Email: glootz@gmail.com 进入到插件根目录下: Npm publish即可 Package.j ...
- springMvc + Maven 项目提示 hessian 依赖包 无法下载;
首先 从 https://github.com/alibaba/dubbo/archive/master.zip 下载最新的 dubbo 源码包到本地某个目录, 解压出来: cmd 进入该目录: 执行 ...
- HTML5 移动开发(移动设备检测及对HTML5的支持)
1.如何选择要使用的特性以及所面向的浏览器 2.哪些浏览器支持HTML5 3.如何检测是否支持HTML5 4.如何开发贷容错性的Web应用程序 5.CSS3媒体查询如何增强检测脚本 使用HTML5 ...
- CodeForces 1059B
Description Student Andrey has been skipping physical education lessons for the whole term, and now ...
- 【SVN】命令行忽略不必要的文件和文件夹
SVN命令参考: https://www.cnblogs.com/wlsxmhz/p/5775393.html 我们需要明白命令行设置忽略文件和文件夹是通过设置svn:ignore属性设置的,pr ...
- python3操作数据库 借助pycharm快速连接并操作mysql数据库
1.https://blog.csdn.net/qiushi_1990/article/details/78041299