洛谷P4069 [SDOI2016]游戏(李超线段树)
题意
Sol
这题细节好多啊qwq。。稍不留神写出一个小bug就要调1h+。。
思路就不多说了,把询问区间拆成两段就是李超线段树板子题了。
关于dis的问题可以直接维护。
// luogu-judger-enable-o2
/*
李超线段树板子题
*/
#include<bits/stdc++.h>
#define Pair pair<int, int>
#define MP(x, y) make_pair(x, y)
#define fi first
#define se second
#define int long long
#define LL long long
#define Fin(x) {freopen(#x".in","r",stdin);}
#define Fout(x) {freopen(#x".out","w",stdout);}
using namespace std;
const int MAXN = 1e6 + 10;
template <typename A, typename B> inline bool chmin(A &a, B b){if(a > b) {a = b; return 1;} return 0;}
template <typename A, typename B> inline bool chmax(A &a, B b){if(a < b) {a = b; return 1;} return 0;}
inline int read() {
char c = getchar(); int x = 0, f = 1;
while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
return x * f;
}
int N, M, dis[MAXN], son[MAXN], fa[MAXN], top[MAXN], dep[MAXN], siz[MAXN], ID[MAXN], rev[MAXN], root, tot;
LL val[MAXN];
vector<Pair> v[MAXN];
struct Line {
int k, b;
}s[MAXN];
int ls[MAXN], rs[MAXN], ll[MAXN], rr[MAXN], mn[MAXN];
void dfs1(int x, int _fa) {
fa[x] = _fa; dep[x] = dep[_fa] + 1; siz[x] = 1;
for(auto &to : v[x]) {
if(to.fi == _fa) continue;
dis[to.fi] = dis[x] + to.se;
dfs1(to.fi, x);
siz[x] += siz[to.fi];
if(siz[to.fi] > siz[son[x]]) son[x] = to.fi;
}
}
void dfs2(int x, int topf) {
top[x] = topf; ID[x] = ++ID[0]; rev[ID[0]] = x;
if(!son[x]) return ;
dfs2(son[x], topf);
for(auto &to : v[x]) {
if(top[to.fi]) continue;
dfs2(to.fi, to.fi);
}
}
int LCA(int x, int y) {
while(top[x] ^ top[y]) {
if(dep[top[x]] < dep[top[y]]) swap(x, y);
x = fa[top[x]];
}
return dep[x] > dep[y] ? y : x;
}
int calc(Line s, int x) {
return s.k * dis[rev[x]] + s.b;
}
double GetPoint(Line x, Line y) {
return (double) (y.b - x.b) / (x.k - y.k);
}
void update(int k) {
if(s[k].k || s[k].b) chmin(mn[k], min(calc(s[k], ll[k]), calc(s[k], rr[k])));
chmin(mn[k], min(mn[ls[k]], mn[rs[k]]));
}
void IntAdd(int &k, int l, int r, int ql, int qr, Line seg) {
if(!k) k = ++tot, mn[k] = val[0], ll[k] = l, rr[k] = r;
int mid = l + r >> 1;
if(ql <= l && r <= qr) {
int pl = calc(s[k], l), pr = calc(s[k], r), nl = calc(seg, l), nr = calc(seg, r);
if(!s[k].k && !s[k].b) {s[k] = seg; update(k); return ;}
if(pl <= nl && pr <= nr) return ;//这里必须要写等号
if(pl > nl && pr > nr) {s[k] = seg; update(k); return ;}
double xx = GetPoint(s[k], seg);
int m = dis[rev[mid]];
if(pl > nl) {
if(xx > m) IntAdd(rs[k], mid + 1, r, ql, qr, s[k]), s[k] = seg;
else IntAdd(ls[k], l, mid, ql, qr, seg);
} else {
if(xx > m) IntAdd(rs[k], mid + 1, r, ql, qr, seg);
else IntAdd(ls[k], l, mid, ql, qr, s[k]), s[k] = seg;
}
update(k);
return ;
}
if(ql <= mid) IntAdd(ls[k], l, mid, ql, qr, seg);
if(qr > mid) IntAdd(rs[k], mid + 1, r, ql, qr, seg);
update(k);
}
int IntMin(int k, int l, int r, int ql, int qr) {
int ret = val[0];
if(ql <= l && r <= qr)
return mn[k];
int mid = l + r >> 1;
if(s[k].k || s[k].b)
chmin(ret, min(calc(s[k], max(l, ql)), calc(s[k], min(r, qr))));
if(ql <= mid) chmin(ret, IntMin(ls[k], l, mid, ql, qr));
if(qr > mid) chmin(ret, IntMin(rs[k], mid + 1, r, ql, qr));
return ret;
}
void Modify(int x, int y, int k, int b) {
while(top[x] != top[y]) {
if(dep[top[x]] < dep[top[y]]) swap(x, y);
IntAdd(root, 1, N, ID[top[x]], ID[x], {k, b});
x = fa[top[x]];
}
if(dep[x] < dep[y]) swap(x, y);
IntAdd(root, 1, N, ID[y], ID[x], {k, b});
}
void Change(int s, int t, int a, int b) {
int lca = LCA(s, t);
Modify(s, lca, -a, b + a * dis[s]);
Modify(t, lca, a, -a * dis[lca] + b + a * (dis[s] - dis[lca]));
}
LL Query(int x, int y) {
int ans = val[0];
while(top[x] != top[y]) {
if(dep[top[x]] < dep[top[y]]) swap(x, y);
chmin(ans, IntMin(root, 1, N, ID[top[x]], ID[x]));
x = fa[top[x]];
}
if(dep[x] < dep[y]) swap(x, y);
chmin(ans, IntMin(root, 1, N, ID[y], ID[x]));
return ans;
}
signed main() {
N = read(); M = read();
for(int i = 0; i <= N; i++) mn[i] = val[i] = 123456789123456789ll;
for(int i = 1; i <= N - 1; i++) {
int x = read(), y = read(), z = read();
v[x].push_back(MP(y, z));
v[y].push_back(MP(x, z));
}
dfs1(1, 0);
dfs2(1, 1);
int GG =0;
for(int i = 1; i <= M; i++) {
int opt = read(), s = read(), t = read();
if(opt == 1) {
int a = read(), b = read();
Change(s, t, a, b);
} else {
cout << Query(s, t) << '\n';
}
}
return 0;
}
洛谷P4069 [SDOI2016]游戏(李超线段树)的更多相关文章
- 洛谷P4069 [SDOI2016]游戏(李超线段树)
题面 传送门 题解 如果我们把路径拆成两段,那么这个路径加可以看成是一个一次函数 具体来说,设\(dis_u\)表示节点\(u\)到根节点的距离,那么\((x,lca)\)这条路径上每个节点的权值就会 ...
- 【洛谷P4097】Segment 李超线段树
题目大意:维护一个二维平面,给定若干条线段,支持询问任意整数横坐标处对应的纵坐标最靠上的线段的 id,相同高度取 id 值较小的,强制在线. 题解:初步学习了李超线段树.李超线段树的核心思想在于通过标 ...
- 洛谷 P3373 【模板】线段树 2
洛谷 P3373 [模板]线段树 2 洛谷传送门 题目描述 如题,已知一个数列,你需要进行下面三种操作: 将某区间每一个数乘上 xx 将某区间每一个数加上 xx 求出某区间每一个数的和 输入格式 第一 ...
- 【BZOJ-4515】游戏 李超线段树 + 树链剖分 + 半平面交
4515: [Sdoi2016]游戏 Time Limit: 40 Sec Memory Limit: 256 MBSubmit: 304 Solved: 129[Submit][Status][ ...
- 洛谷 P2574 XOR的艺术(线段树 区间异或 区间求和)
To 洛谷.2574 XOR的艺术 题目描述 AKN觉得第一题太水了,不屑于写第一题,所以他又玩起了新的游戏.在游戏中,他发现,这个游戏的伤害计算有一个规律,规律如下 1. 拥有一个伤害串为长度为n的 ...
- 洛谷P3372 【模板】线段树 1
P3372 [模板]线段树 1 153通过 525提交 题目提供者HansBug 标签 难度普及+/提高 提交 讨论 题解 最新讨论 [模板]线段树1(AAAAAAAAA- [模板]线段树1 洛谷 ...
- 洛谷P4891 序列(势能线段树)
洛谷题目传送门 闲话 考场上一眼看出这是个毒瘤线段树准备杠题,发现实在太难调了,被各路神犇虐哭qwq 考后看到各种优雅的暴力AC......宝宝心里苦qwq 思路分析 题面里面是一堆乱七八糟的限制和性 ...
- 洛谷P4344 脑洞治疗仪 [SHOI2015] 线段树+二分答案/分块
!!!一道巨恶心的数据结构题,做完当场爆炸:) 首先,如果你用位运算的时候不小心<<打成>>了,你就可以像我一样陷入疯狂的死循环改半个小时 然后,如果你改出来之后忘记把陷入死循 ...
- Bzoj5294/洛谷P4428 [Bjoi2018]二进制(线段树)
题面 Bzoj 洛谷 题解 考虑一个什么样的区间满足重组之后可以变成\(3\)的倍数.不妨设\(tot\)为一个区间内\(1\)的个数.如果\(tot\)是个偶数,则这个区间一定是\(3\)的倍数,接 ...
随机推荐
- linux防火墙(三)—— iptables语法之匹配条件
一.iptables规则的匹配条件类型有三类 1.通用匹配:可直接使用,不依赖于其他条件或扩展,包括网络协议.IP地址.网络接口等条件 2.隐含匹配:要求以特定的协议匹配作为前提,包括端口.TCP标记 ...
- java项目配置域名(tomcat直接配置 or 使用nginx反向代理)
一: tomcat直接配置域名:https://blog.csdn.net/qq_36330228/article/details/78516160 二: 使用nginx进行反向代理 tomcat服 ...
- centos7 初始化脚本
#!/bin/bash # 时间: 2018-11-21 # 作者: HuYuan # 描述: CentOS 7 初始化脚本 # 加载配置文件 if [ -n "${1}" ];t ...
- POJ 2501
#include<iostream> #include<iomanip> #include<stdio.h> #include<string> #inc ...
- odoo开发笔记--模型字段compute用法
compute属性,实现的主要功能是,前端界面选择某个字段的时候,指定与该字段关联的其他字段可以关联,并联动的显示. 可以和inverse属性同时使用,不加inverse属性的话,前端界面的显示效果只 ...
- 鼠标右键添加"在此处打开命令窗口"
从windows7开始,提供了一个便于从当前文件夹打开cmd命令行窗口的快捷方式: 直接使用 Shift+鼠标右键==>“在此处打开命令行窗口” 有可能以后会遇到的问题,可参考如下.. htt ...
- 整理几个 RPC 框架
gRPChttp://www.grpc.io/https://github.com/grpcMotan支撑微博千亿调用的轻量级RPC框架:Motanhttp://h2ex.com/820WangleW ...
- 删除qq互联
1.删除source/plugin 下的文件qqconnect文件夹 2.删除数据: DELETE FROM `dn_common_plugin` WHERE `identifier` = 'qqco ...
- java SynchronousQueue
SynchronousQueue 内部没有容量,但是由于一个插入操作总是对应一个移除操作,反过来同样需要满足.那么一个元素就不会再SynchronousQueue 里面长时间停留,一旦有了插入线程和移 ...
- springboot-32-使用mvc测试
Spring Boot可以和大部分流行的测试框架协同工作:通过Spring JUnit创建单元测试:生成测试数据初始化数据库用于测试:Spring Boot可以跟BDD(Behavier Driven ...