题意

\(n(2 \le n \le 500000)\)个点\(m(1 \le m \le 1000000)\)条边的有向无环图,找到一个点,使得删掉这个点后剩余图中的最长路径最短。

分析

神题不会做。

题解

首先我们新建个源\(s\)和汇\(t\),连边\(s->i, i->t\),最远距离分别为\(d[i, 0]\)和\(d[i, 1]\),则一个图中的最长链就是\(max(d[u, 0]+d[v, 1]-1, \exists edge(u, v))\),再由于图中任意一个\(s-t\)割都会有最长链的边,那么问题就转化为去掉\(x\)点及其相关的边后,在途中的一个割中找最大值。由于图是拓扑序的,所以,我们只要把某割点的入边拿掉后、拓扑秩在自己前或相同的点出边都加入到一个堆里,那么堆中就可以形成一个割!这个比较显然。

#include <bits/stdc++.h>
using namespace std;
inline int getint() {
int x=0;
char c=getchar();
for(; c<'0'||c>'9'; c=getchar());
for(; c>='0'&&c<='9'; x=x*10+c-'0', c=getchar());
return x;
}
const int N=500005;
int q[N], in[N], d[N][2], n, cnt1, cnt2, ihead1[N], ihead2[N], sz[N<<2];
struct E {
int next, to;
}e1[N<<2], e2[N<<2];
void add(int x, int y) {
e1[++cnt1]=(E){ihead1[x], y}; ihead1[x]=cnt1;
e2[++cnt2]=(E){ihead2[y], x}; ihead2[y]=cnt2;
++in[y];
}
void ins(int p, int d, int l=0, int r=n, int x=1) {
sz[x]+=d;
if(l==r) {
return;
}
int mid=(l+r)>>1;
p<=mid?ins(p, d, l, mid, x<<1):ins(p, d, mid+1, r, x<<1|1);
}
int ask(int l=0, int r=n, int x=1) {
if(l==r) {
return l;
}
int mid=(l+r)>>1, ls=x<<1, rs=ls|1;
return sz[rs]?ask(mid+1, r, rs):ask(l, mid, ls);
}
int main() {
n=getint();
int m=getint(), fr=0, ta=0;
for(int i=1; i<=m; ++i) {
int x=getint(), y=getint();
add(x, y);
}
for(int i=1; i<=n; ++i) {
if(!in[i]) {
q[ta++]=i;
}
}
while(fr!=ta) {
int x=q[fr++];
for(int i=ihead1[x]; i; i=e1[i].next) {
if(!--in[e1[i].to]) {
q[ta++]=e1[i].to;
}
}
}
for(int j=0; j<n; ++j) {
for(int x=q[j], i=ihead2[x]; i; i=e2[i].next) {
d[x][0]=max(d[x][0], d[e2[i].to][0]+1);
}
for(int x=q[n-j-1], i=ihead1[x]; i; i=e1[i].next) {
d[x][1]=max(d[x][1], d[e1[i].to][1]+1);
}
}
for(int i=1, S=n+1, T=S+1; i<=n; ++i) {
add(S, i);
add(i, T);
ins(d[i][1], 1);
}
d[n+1][0]=-1;
d[n+2][1]=-1;
int ans1=-1, ans2=n+1;
for(int j=0; j<n; ++j) {
int x=q[j], t;
for(int i=ihead2[x]; i; i=e2[i].next) {
ins(d[e2[i].to][0]+d[x][1]+1, -1);
}
if(ans2>(t=ask())) {
ans1=x;
ans2=t;
}
for(int i=ihead1[x]; i; i=e1[i].next) {
ins(d[x][0]+d[e1[i].to][1]+1, 1);
}
}
printf("%d %d\n", ans1, ans2);
return 0;
}

【BZOJ】3832: [Poi2014]Rally的更多相关文章

  1. 【BZOJ】3524: [Poi2014]Couriers

    [算法]主席树 [题解]例题,记录和,数字出现超过一半就递归查找. 主席树见[算法]数据结构 #include<cstdio> #include<algorithm> #inc ...

  2. 【BZOJ】3835: [Poi2014]Supercomputer

    题意 \(n(1 \le 1000000)\)个点的有根树,\(1\)号点为根,\(q(1 \le 1000000)\)次询问,每次给一个\(k\),每一次可以选择\(k\)个未访问的点,且父亲是访问 ...

  3. 【BZOJ】3526: [Poi2014]Card

    题意 \(n(n \le 200000)\)张卡片,正反有两个数\(a[i], b[i]\).\(m(m \le 1000000)\)次操作,每次交换\(c[i].d[i]\)位置上的卡片.每一次操作 ...

  4. 【BZOJ】3523: [Poi2014]Bricks

    题意 \(n(n \le 1000000)\)个物品,颜色分别为\(a[i]\),现在要求排在一排使得相邻两个砖块的颜色不同,且限定第一个砖块和最后一个砖块的颜色,输出一个合法解否则输出-1. 分析 ...

  5. 【BZOJ】3521: [Poi2014]Salad Bar

    题意 长度为\(n(1 \le n \le 1000000)\)的\(01\)字符串.找一个最长的连续子串\(S\),使得不管是从左往右还是从右往左取,都保证每时每刻已取出的\(1\)的个数不小于\( ...

  6. 【BZOJ】3834: [Poi2014]Solar Panels

    http://www.lydsy.com/JudgeOnline/problem.php?id=3834 题意:求$max\{(i,j)\}, smin<=i<=smax, wmin< ...

  7. BZOJ:3832: [Poi2014]Rally

    题意: 给出$DAG$,询问删掉哪个点之后最长路径最短 思路: 我们令$f[x]$表示从最远的点到达它的距离,$g[x]$表示它能够到达最远的点的距离 那么对于$(x -> y)$一条边来说,它 ...

  8. 【BZOJ】3524 [POI2014] Couriers(主席树)

    题目 传送门:QWQ 传送到洛谷QWQ 分析 把求区间第k大的改一改就ok了. 代码 #include <bits/stdc++.h> using namespace std; ; ], ...

  9. 【BZOJ】3052: [wc2013]糖果公园

    http://www.lydsy.com/JudgeOnline/problem.php?id=3052 题意:n个带颜色的点(m种),q次询问,每次询问x到y的路径上sum{w[次数]*v[颜色]} ...

随机推荐

  1. [Unity] 导出Android APK包出错

    确认Android环境是OK的. 检查 StreamingAssets 目录下是否有中文的文件名 检查其它目录的中文文件名. 移除一些插件再试.

  2. Redis学习笔记7--Redis管道(pipeline)

    redis是一个cs模式的tcp server,使用和http类似的请求响应协议.一个client可以通过一个socket连接发起多个请求命令.每个请求命令发出后client通常会阻塞并等待redis ...

  3. AP(affinity propagation)研究

    待补充…… AP算法,即Affinity propagation,是Brendan J. Frey* 和Delbert Dueck于2007年在science上提出的一种算法(文章链接,维基百科) 现 ...

  4. windows环境下面安装Apache2.4+MySql5.7+PHP5.6

    之前学习PHP一致是只用phpStudy集成开发环境,这对于新手而言无疑是帮助极大的,因为傻瓜式安装与使用方法减少了我们很多不必要的麻烦.但是作为一名合格的PHP开发人员,掌握PHP+MySQL+Ap ...

  5. 第3月第9天 循环引用 block

    一.一个对象没有被引用,那么在函数块完成时就会被dealloc,这种情况因为对象销毁了,block块也永远不会执行. MyNetworkOperation *op = [[MyNetworkOpera ...

  6. Linux下c开发 之 线程通信(转)

    Linux下c开发 之 线程通信(转) 1.Linux“线程” 进程与线程之间是有区别的,不过Linux内核只提供了轻量进程的支持,未实现线程模型.Linux是一种“多进程单线程”的操作系统.Linu ...

  7. Java Native Interface 二 JNI中对Java基本类型和引用类型的处理

    本文是<The Java Native Interface Programmer's Guide and Specification>读书笔记 Java编程里会使用到两种类型:基本类型(如 ...

  8. Qt界面中嵌入其他exe程序的界面,使用Qt5

    下面用一个小例子来演示如何在Qt的界面中嵌入其他exe程序的界面,最终效果如下图所示.本文参考了 http://blog.csdn.net/jiaoyaziyang/article/details/4 ...

  9. IOS单元测试,OCMock

    1.环境搭建 使用官方推荐,手动搭建.. 1.下载所需文件 http://ocmock.org/download/ 根据自己的需要下载不同的版本. 打开iOS library文件夹,里面就有需要用到的 ...

  10. TypeError: datetime.datetime(2016, 9, 25, 21, 12, 19, 135649) is not JSON serializable解决办法

    1.一个简单的方法来修补json模块,这样序列将支持日期时间. import json import datetime json.JSONEncoder.default = lambda self, ...