[luoguP2680] 运输计划(lca + 二分 + 差分)
暴力做法 50 ~ 60
枚举删边,求最大路径长度的最小值。
其中最大路径长度运用到了lca
我们发现,求lca的过程已经不能优化了,那么看看枚举删边的过程能不能优化。
先把边按照权值排序,然后。。。
然而并没有什么卵用。
题解给出了一种二分
二分答案。
判断依据:
比当前答案大的路径长度如果有一个公共边,并且最大的路径减去这条公共边小于等于当前答案,那么当前答案就满足。
如果当前答案不满足,那么比他小的答案更不可能满足,因为都不小于等于当前答案,比当前答案小的更不可能。
接下来就是如何判断,num[i] 表示点 i 指向父亲的边被几条路径所共有。
那么只需要 num[x]++, num[y]++, num[lca(x, y)] -= 2
最后dfs一遍求树上前缀和即可。
然后在 (num[i] == 比当前答案大的边的条数) 中找最大的边
然后判断就可以了
但是就是一个点超时,无语,把倍增改成tarjan也是超时,蛋疼。
——代码
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define max(x, y) ((x) > (y) ? (x) : (y))
#define MAXN 300001 int n, m, cnt, qcnt;
int num[MAXN], dis[MAXN], fv[MAXN], fa[MAXN], f[MAXN];
int head[MAXN], to[MAXN << ], val[MAXN << ], next[MAXN << ], qhead[MAXN], qnext[MAXN << ], qto[MAXN << ];
//int f[MAXN][21], deep[MAXN]; struct node
{
int qx, qy, lca, v;
}q[MAXN]; inline int read()
{
int x = , f = ;
char ch = getchar();
for(; !isdigit(ch); ch = getchar()) if(ch == '-') f = -;
for(; isdigit(ch); ch = getchar()) x = (x << ) + (x << ) + ch - '';
return x * f;
} inline void add(int x, int y, int z)
{
to[cnt] = y;
val[cnt] = z;
next[cnt] = head[x];
head[x] = cnt++;
} inline void addq(int x, int y)
{
qto[qcnt] = y;
qnext[qcnt] = qhead[x];
qhead[x] = qcnt++;
} /*inline void dfs(int u)
{
int i, v;
deep[u] = deep[f[u][0]] + 1;
for(i = 0; f[u][i]; i++) f[u][i + 1] = f[f[u][i]][i];
for(i = head[u]; i ^ -1; i = next[i])
{
v = to[i];
if(!deep[v]) fv[v] = val[i], dis[v] = dis[u] + val[i], f[v][0] = u, dfs(v);
}
} inline int Lca(int x, int y)
{
int i;
if(deep[x] < deep[y]) x ^= y ^= x ^= y;
for(i = 18; i >= 0; i--)
if(deep[f[x][i]] >= deep[y])
x = f[x][i];
if(x == y) return x;
for(i = 18; i >= 0; i--)
if(f[x][i] ^ f[y][i])
x = f[x][i], y = f[y][i];
return f[x][0];
}*/ inline int find(int x)
{
return x == fa[x] ? x : fa[x] = find(fa[x]);
} inline void dfs(int u)
{
int i, v, x;
fa[u] = u;
for(i = head[u]; i ^ -; i = next[i])
{
v = to[i];
if(f[u] ^ v) f[v] = u, dis[v] = dis[u] + val[i], fv[v] = val[i], dfs(v);
}
for(i = qhead[u]; i ^ -; i = qnext[i])
{
v = qto[i];
if(f[x = u ^ q[v].qx ^ q[v].qy] || x == ) q[v].lca = find(x);
}
fa[u] = f[u];
} inline void dfs1(int u)
{
int i, v;
for(i = head[u]; i ^ -; i = next[i])
{
v = to[i];
if(v ^ f[u]) dfs1(v), num[u] += num[v];
}
} inline bool pd(int sum)
{
int i, x = , y = m, mid, ans, value = ;
memset(num, , sizeof(num));
while(x <= y)
{
mid = (x + y) >> ;
if(q[mid].v <= sum) x = mid + ;
else ans = mid, y = mid - ;
}
for(i = ans; i <= m; i++) num[q[i].qx]++, num[q[i].qy]++, num[q[i].lca] -= ;
dfs1();
for(i = ; i <= n; i++)
if(num[i] == m - ans + )
value = max(value, fv[i]);
if(!value) return ;
return q[m].v - value <= sum;
} inline bool cmp(node x, node y)
{
return x.v < y.v;
} int main()
{
int i, x, y, z, mid, ans = ;
n = read();
m = read();
memset(head, -, sizeof(head));
memset(qhead, -, sizeof(qhead));
for(i = ; i < n; i++)
{
x = read();
y = read();
z = read();
add(x, y, z);
add(y, x, z);
}
//dfs(1);
for(i = ; i <= m; i++)
{
q[i].qx = read();
q[i].qy = read();
addq(q[i].qx, i);
addq(q[i].qy, i);
//q[i].lca = Lca(q[i].qx, q[i].qy);
//q[i].v = dis[q[i].qx] + dis[q[i].qy] - (dis[q[i].lca] << 1);
}
dfs();
x = y = ;
for(i = ; i <= m; i++) q[i].v = dis[q[i].qx] + dis[q[i].qy] - (dis[q[i].lca] << ), y = max(y, q[i].v);
std::sort(q + , q + m + , cmp);
while(x <= y)
{
mid = (x + y) >> ;
if(pd(mid)) ans = mid, y = mid - ;
else x = mid + ;
}
printf("%d\n", ans);
return ;
}
还留有倍增的痕迹。
[luoguP2680] 运输计划(lca + 二分 + 差分)的更多相关文章
- BZOJ 4326:NOIP2015 运输计划(二分+差分+lca)
NOIP2015 运输计划Description公元 2044 年,人类进入了宇宙纪元.L 国有 n 个星球,还有 n−1 条双向航道,每条航道建立在两个星球之间,这 n−1 条航道连通了 L 国的所 ...
- [luogu2680] 运输计划 (lca+二分+树上差分)
传送门 Description Input Output 一个整数,表示小 P 的物流公司完成阶段性工作所需要的最短时间. Sample Input 6 3 1 2 3 1 6 4 3 1 7 4 3 ...
- luoguP2680 运输计划 题解(二分答案+树上差分)
P2680 运输计划 题目 这道题如果是看的我的树上差分来的,那么肯定一看题目就可以想到树上差分. 至于这是怎么想到的,一步一步来: 1.n有300000,不可能暴力枚举每一条边 2.因为我们要使运 ...
- NOIP2015 运输计划(二分+LCA+差分)
4326: NOIP2015 运输计划 Time Limit: 30 Sec Memory Limit: 128 MBSubmit: 308 Solved: 208[Submit][Status] ...
- BZOJ 4326 NOIP2015 运输计划(树上差分+LCA+二分答案)
4326: NOIP2015 运输计划 Time Limit: 30 Sec Memory Limit: 128 MB Submit: 1388 Solved: 860 [Submit][Stat ...
- [NOIP2015]运输计划(树上差分+LCA+二分)
Description 公元 2044 年,人类进入了宇宙纪元. L 国有 n 个星球,还有 n−1 条双向航道,每条航道建立在两个星球之间,这 n−1 条航道连通了 L 国的所有星球. 小 P 掌管 ...
- bzoj4326: NOIP2015 运输计划(二分+LCA+树上差分)
题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=4326 题目大意:有一颗含有n个顶点的树,每两个点之间有一个边权,现在有m个运输计划,每个 ...
- 洛谷P2680 运输计划 [LCA,树上差分,二分答案]
题目传送门 运输计划 Description 公元 2044 年,人类进入了宇宙纪元.L 国有 n 个星球,还有 n?1 条双向航道,每条航道建立在两个星球之间, 这 n?1 条航道连通了 L 国的所 ...
- LOJ2425 NOIP2015 运输计划 【二分+LCA+树上差分】*
LOJ2425 NOIP2015 运输计划 LINK 题意:给你一颗树,可以将任意一条边的权值变成0,然后求m条路径的长度的最小值 思路: 先二分最后的距离ans,然后我们把路程大于ans的所有路径拿 ...
- BZOJ_4326_[NOIP2015]_运输计划_(二分+LCA_树链剖分/Tarjan+差分)
描述 http://www.lydsy.com/JudgeOnline/problem.php?id=4326 给出一棵带有边权的树,以及一系列任务,任务是从树上的u点走到v点,代价为u到v路径上的权 ...
随机推荐
- bzoj 1611: [Usaco2008 Feb]Meteor Shower流星雨【BFS】
t记录每个格子最早被砸的时间,bfs(x,y,t)表示当前状态为(x,y)格子,时间为t.因为bfs,所以先搜到的t一定小于后搜到的,所以一个格子搜一次就行 #include<iostream& ...
- Manacher HDOJ 3068 最长回文
题目传送门 关于求解最长回文子串,有dp做法,也有同样n^2的但只用O(1)的空间,还有KMP,后缀数组?? int main(void) { ) == ) { ); memset (dp, fals ...
- 思维题 URAL 1409 Two Gangsters
题目传送门 /* 思维题:注意题目一句话:At some moment it happened so that they shot one and the same can. 如果两个人都有射中的话, ...
- Javascript对象基础讲解
1.Object对象详解 javascript 里最基本的数据类型是对象. avaScript里的对象其实是一个无序的属性集合,属性又是一个个的名-值对. 除了字符串,数字,true,false,nu ...
- Android 性能优化(27)*zipalign让apk数据对齐,运行更快。
1.zipalign 简介 zipalign is an archive alignment tool that provides important optimization to Android ...
- 【IIS7.5】Asp文件上传限制,加载页面大小限制
运行环境 window server 2008 R2 X64.IIS7.5.应用程序池.Net4.0 经典模式 分析 IIS7.5默认有两处上传限制: 第一处在,选择左侧的站点,然后找到后侧的管理—— ...
- 专题四:自定义Web浏览器
前言: 前一个专题介绍了自定义的Web服务器,然而向Web服务器发出请求的正是本专题要介绍的Web浏览器,本专题通过简单自定义一个Web浏览器来简单介绍浏览器的工作原理,以及帮助一些初学者揭开浏览器这 ...
- 关于Python多线程condition变量的应用
''' 所谓条件变量,即这种机制是在满足了特定的条件后,线程才可以访问相关的数据. 它使用Condition类来完成,由于它也可以像锁机制那样用,所以它也有acquire方法和release方法,而且 ...
- win7下安装MySQL 5.7.19(解压缩版)
1.官网下载地址:https://downloads.mysql.com/archives/community/ 下载后,得到压缩包: 2.解压,我的解压目录为:E:\mysql-5.7.19\mys ...
- 用纯函数式思维在Java8下写的一段奇葩程序
首先说一下什么是纯函数式.在我的理解,"纯函数式"用一句话就可以描述:Anything is value.--我的理解不一定准确,但我就是这么理解的. 就是所有的东西都是值--没有 ...