BZOJ2157 旅游 【树剖 或 LCT】

最新推荐文章于 2024-04-24 17:27:11 发布
最新推荐文章于 2024-04-24 17:27:11 发布
阅读量260 收藏
点赞数
分类专栏： 题解 数据结构 树链剖分 LCT 线段树
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_38678604/article/details/78954581
版权
题解 同时被 3 个专栏收录
169 篇文章 1 订阅
订阅专栏
数据结构
31 篇文章 0 订阅
订阅专栏
线段树
8 篇文章 0 订阅
订阅专栏

题目

Ray 乐忠于旅游，这次他来到了T 城。T 城是一个水上城市，一共有 N 个景点，有些景点之间会用一座桥连接。为了方便游客到达每个景点但又为了节约成本，T 城的任意两个景点之间有且只有一条路径。换句话说， T 城中只有N − 1 座桥。Ray 发现，有些桥上可以看到美丽的景色，让人心情愉悦，但有些桥狭窄泥泞，令人烦躁。于是，他给每座桥定义一个愉悦度w，也就是说，Ray 经过这座桥会增加w 的愉悦度，这或许是正的也可能是负的。有时，Ray 看待同一座桥的心情也会发生改变。现在，Ray 想让你帮他计算从u 景点到v 景点能获得的总愉悦度。有时，他还想知道某段路上最美丽的桥所提供的最大愉悦度，或是某段路上最糟糕的一座桥提供的最低愉悦度。

输入格式

输入的第一行包含一个整数N，表示T 城中的景点个数。景点编号为 0…N − 1。接下来N − 1 行，每行三个整数u、v 和w，表示有一条u 到v，使 Ray 愉悦度增加w 的桥。桥的编号为1…N − 1。|w| <= 1000。输入的第N + 1 行包含一个整数M，表示Ray 的操作数目。接下来有M 行，每行描述了一个操作，操作有如下五种形式： C i w，表示Ray 对于经过第i 座桥的愉悦度变成了w。 N u v，表示Ray 对于经过景点u 到v 的路径上的每一座桥的愉悦度都变成原来的相反数。 SUM u v，表示询问从景点u 到v 所获得的总愉悦度。 MAX u v，表示询问从景点u 到v 的路径上的所有桥中某一座桥所提供的最大愉悦度。 MIN u v，表示询问从景点u 到v 的路径上的所有桥中某一座桥所提供的最小愉悦度。测试数据保证，任意时刻，Ray 对于经过每一座桥的愉悦度的绝对值小于等于1000。

输出格式

对于每一个询问(操作S、MAX 和MIN)，输出答案。

输入样例

3

0 1 1

1 2 2

8

SUM 0 2

MAX 0 2

N 0 1

SUM 0 2

MIN 0 2

C 1 3

SUM 0 2

MAX 0 2

输出样例

3

2

1

-1

5

3

提示

一共有10 个数据，对于第i (1 <= i <= 10) 个数据， N = M = i * 2000。

题解

凭着对树剖的熟练度，20min打完180+行直接1A【只会打板】
学了LCT发现其实树剖的题目都可以用LCT刷刷
LCT明天补上【如果我查得完错】
UPDATE：已补LCT
果然树剖快挺多
这里写图片描述

#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define REP(i,n) for (int i = 1; i <= (n); i++)
#define Redge(u) for (int k = h[u]; k != -1; k = ed[k].nxt)
#define ls (u << 1)
#define rs (u << 1 | 1)
using namespace std;
const int maxn = 20005,maxm = 100005,INF = 1000000000;
inline int RD(){
    int out = 0,flag = 1; char c = getchar();
    while (c < 48 || c > 57) {if (c == '-') flag = -1; c = getchar();}
    while (c >= 48 && c <= 57) {out = (out << 3) + (out << 1) + c - '0'; c = getchar();}
    return out * flag;
}
char opt[10];
int h[maxn],ne = 0,N,M;
struct EDGE{int to,nxt,w,i;}ed[maxm];
inline void build(int u,int v,int w,int i){
    ed[ne] = (EDGE){v,h[u],w,i}; h[u] = ne++;
    ed[ne] = (EDGE){u,h[v],w,i}; h[v] = ne++;
}
int fa[maxn],dep[maxn],son[maxn],V[maxn],id[maxn],P[maxn],H[maxn],top[maxn],siz[maxn];
int cnt = 0;
void dfs1(int u){
    siz[u] = 1; int to;
    Redge(u) if ((to = ed[k].to) != fa[u]){
        dep[to] = dep[u] + 1; fa[to] = u; P[ed[k].i] = to; V[to] = ed[k].w;
        dfs1(to);
        siz[u] += siz[to];
        if (!son[u] || siz[to] > siz[son[u]]) son[u] = to;
    }
}
void dfs2(int u,int flag){
    id[u] = ++cnt; H[cnt] = u; top[u] = flag ? top[fa[u]] : u; int to;
    if (son[u]) dfs2(son[u],true);
    Redge(u) if ((to = ed[k].to) != fa[u] && to != son[u]) dfs2(to,false);
}
int sum[4 * maxn],mx[4 * maxn],mn[4 * maxn],tag[4 * maxn],L,R;
void pup(int u){
    sum[u] = sum[ls] + sum[rs];
    mx[u] = max(mx[ls],mx[rs]);
    mn[u] = min(mn[ls],mn[rs]);
}
void pd(int u){
    if (tag[u]){
        sum[ls] = -sum[ls]; sum[rs] = -sum[rs];
        swap(mx[ls],mn[ls]); mx[ls] = -mx[ls]; mn[ls] = -mn[ls];
        swap(mx[rs],mn[rs]); mx[rs] = -mx[rs]; mn[rs] = -mn[rs];
        tag[ls] ^= 1; tag[rs] ^= 1; tag[u] ^=1;
    }
}
void build(int u,int l,int r){
    if (l == r){
        sum[u] = mx[u] = mn[u] = V[H[l]];
        return;
    }
    int mid = l + r >> 1;
    build(ls,l,mid);
    build(rs,mid + 1,r);
    pup(u);
}
void modify(int u,int l,int r,int v){
    if (l == r) {sum[u] = mx[u] = mn[u] = v; return;}
    pd(u);
    int mid = l + r >> 1;
    if (mid >= L) modify(ls,l,mid,v);
    else modify(rs,mid + 1,r,v);
    pup(u);
}
void rever(int u,int l,int r){
    if (l >= L && r <= R){
        sum[u] = -sum[u];
        swap(mx[u],mn[u]); mx[u] = -mx[u]; mn[u] = -mn[u];
        tag[u] ^= 1;
        return;
    }
    pd(u);
    int mid = l + r >> 1;
    if (mid >= L) rever(ls,l,mid);
    if (mid < R) rever(rs,mid + 1,r);
    pup(u);
}
int Query(int u,int l,int r){
    if (l >= L && r <= R) return sum[u];
    pd(u);
    int mid = l + r >> 1;
    if (mid >= R) return Query(ls,l,mid);
    else if (mid < L) return Query(rs,mid + 1,r);
    else return Query(ls,l,mid) + Query(rs,mid + 1,r);
}
int Qmin(int u,int l,int r){
    if (l >= L && r <= R) return mn[u];
    pd(u);
    int mid = l + r >> 1;
    if (mid >= R) return Qmin(ls,l,mid);
    else if (mid < L) return Qmin(rs,mid + 1,r);
    else return min(Qmin(ls,l,mid),Qmin(rs,mid + 1,r));
}
int Qmax(int u,int l,int r){
    if (l >= L && r <= R) return mx[u];
    pd(u);
    int mid = l + r >> 1;
    if (mid >= R) return Qmax(ls,l,mid);
    else if (mid < L) return Qmax(rs,mid + 1,r);
    else return max(Qmax(ls,l,mid),Qmax(rs,mid + 1,r));
}
void solve1(){
    int u = P[RD()],w = RD();
    L = id[u]; modify(1,1,N,w);
}
void solve2(){
    int u = RD() + 1,v = RD() + 1;
    while (top[u] != top[v]){
        if (dep[top[u]] < dep[top[v]]) swap(u,v);
        L = id[top[u]]; R = id[u];
        rever(1,1,N);
        u = fa[top[u]];
    }
    if (u == v) return;
    if (dep[u] > dep[v]) swap(u,v);
    L = id[u] + 1; R = id[v];
    rever(1,1,N);
}
void solve3(){
    int u = RD() + 1,v = RD() + 1,ans = 0;
    while (top[u] != top[v]){
        if (dep[top[u]] < dep[top[v]]) swap(u,v);
        L = id[top[u]]; R = id[u];
        ans += Query(1,1,N);
        u = fa[top[u]];
    }
    if (u != v){
        if (dep[u] > dep[v]) swap(u,v);
        L = id[u] + 1; R = id[v];
        ans += Query(1,1,N);
    }
    printf("%d\n",ans);
}
void solve4(){
    int u = RD() + 1,v = RD() + 1,ans = -INF;
    while (top[u] != top[v]){
        if (dep[top[u]] < dep[top[v]]) swap(u,v);
        L = id[top[u]]; R = id[u];
        ans = max(ans,Qmax(1,1,N));
        u = fa[top[u]];
    }
    if (u != v){
        if (dep[u] > dep[v]) swap(u,v);
        L = id[u] + 1; R = id[v];
        ans = max(ans,Qmax(1,1,N));
    }
    printf("%d\n",ans);
}
void solve5(){
    int u = RD() + 1,v = RD() + 1,ans = INF;
    while (top[u] != top[v]){
        if (dep[top[u]] < dep[top[v]]) swap(u,v);
        L = id[top[u]]; R = id[u];
        ans = min(ans,Qmin(1,1,N));
        u = fa[top[u]];
    }
    if (u != v){
        if (dep[u] > dep[v]) swap(u,v);
        L = id[u] + 1; R = id[v];
        ans = min(ans,Qmin(1,1,N));
    }
    printf("%d\n",ans);
}
int main(){
    memset(h,-1,sizeof(h));
    N = RD(); int u,v,w;
    REP(i,N - 1) u = RD() + 1,v = RD() + 1,w = RD(),build(u,v,w,i);
    dep[1] = 1; dfs1(1); dfs2(1,0); build(1,1,N);
    M = RD();
    while (M--){
        scanf("%s",opt);
        if (opt[0] == 'C') solve1();
        else if (opt[0] == 'N') solve2();
        else if (opt[0] == 'S') solve3();
        else if (opt[1] == 'A') solve4();
        else solve5();
    }
    return 0;
}

LCT：【代码短但会慢一点】

#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define isrt(u) (!f[u] || (ch[f[u]][0] != u && ch[f[u]][1] != u))
#define isr(u) (ch[f[u]][1] == u)
#define REP(i,n) for (int i = 1; i <= (n); i++)
#define Redge(u) for (int k = h[u]; k != -1; k = ed[k].nxt)
#define ls (ch[u][0])
#define rs (ch[u][1])
using namespace std;
const int maxn = 40005,maxm = 100005,INF = 1000000000;
inline int RD(){
    int out = 0,flag = 1; char c = getchar();
    while (c < 48 || c > 57) {if (c == '-') flag = -1; c = getchar();}
    while (c >= 48 && c <= 57) {out = (out << 3) + (out << 1) + c - '0'; c = getchar();}
    return out * flag;
}
int rev[maxn],Min[maxn],Max[maxn],f[maxn],ch[maxn][2],sum[maxn];
int tag[maxn],V[maxn],N,M,ed[maxn];
char opt[10];
void Rever(int u){
    swap(Min[u],Max[u]); Max[u] = -Max[u]; Min[u] = -Min[u];
    sum[u] = -sum[u]; V[u] = -V[u];
    tag[u] ^= 1;
}
void pup(int u){
    Max[u] = max(Max[ls],Max[rs]);
    Min[u] = min(Min[ls],Min[rs]);
    if (u > N) Max[u] = max(Max[u],V[u]);
    if (u > N) Min[u] = min(Min[u],V[u]);
    sum[u] = sum[ls] + sum[rs] + V[u];
}
void pd(int u){
    if (rev[u]){swap(ls,rs); rev[ls] ^= 1; rev[rs] ^= 1; rev[u] ^= 1;}
    if (tag[u]){
        if (ls) Rever(ls);
        if (rs) Rever(rs);
        tag[u] = 0;
    }
}
void push_down(int u){
    if (!isrt(u)) push_down(f[u]);
    pd(u);
}
void spin(int u){
    int s = isr(u),fa = f[u];
    f[u] = f[fa]; if (!isrt(fa)) ch[f[fa]][isr(fa)] = u;
    ch[fa][s] = ch[u][s ^ 1]; if (ch[u][s ^ 1]) f[ch[u][s ^ 1]] = fa;
    f[fa] = u; ch[u][s ^ 1] = fa;
    pup(fa); pup(u);
}
void splay(int u){
    for (push_down(u); !isrt(u); spin(u))
        if (!isrt(f[u])) spin((isr(u) ^ isr(f[u])) ? u : f[u]);
}
void Access(int u){
    for (int v = 0; u; u = f[v = u])
        splay(u),ch[u][1] = v,pup(u);
}
void Make_rt(int u){
    Access(u); splay(u); rev[u] ^= 1;
}
void Split(int u,int v){
    Make_rt(u); Access(v); splay(v);
}
void Link(int u,int v){
    Make_rt(u); f[u] = v;
}
void solve1(){
    int u = ed[RD()],v = RD();
    splay(u); V[u] = v; pup(u);
}
void solve2(){
    int u = RD() + 1,v = RD() + 1;
    Split(u,v); Rever(v);
}
void solve3(){
    int u = RD() + 1,v = RD() + 1;
    Split(u,v);
    printf("%d\n",sum[v]);
}
void solve4(){
    int u = RD() + 1,v = RD() + 1;
    Split(u,v);
    printf("%d\n",Max[v]);
}
void solve5(){
    int u = RD() + 1,v = RD() + 1;
    Split(u,v);
    printf("%d\n",Min[v]);
}
int main(){
    N = RD();
    int u,v,w,id = N;
    for (int i = 0; i <= N; i++) Max[i] = -INF,Min[i] = INF;
    REP(i,N - 1){
        u = RD() + 1,v = RD() + 1,w = RD();
        ed[i] = ++id; V[id] = Max[id] = Min[id] = sum[id] = w;
        Link(u,id); Link(v,id);
    }
    M = RD();
    while (M--){
        scanf("%s",opt);
        if (opt[0] == 'C') solve1();
        else if (opt[0] == 'N') solve2();
        else if (opt[0] == 'S') solve3();
        else if (opt[1] == 'A') solve4();
        else solve5();
    }
    return 0;
}
Mychael
关注
专栏目录
BZOJ 2157 「国家集训队」旅游（树链剖分，线段树，边权转点权）【BZOJ计划】
繁凡さん的博客
09-22 1388
【BZOJ修复计划 #20】BZOJ 2157 旅游 【国家集训队2011】
洛谷P1505 [国家集训队]旅游（BZOJ2157）
forezxl的博客
06-06 350
树链剖分 洛谷题目传送门 BZOJ题目传送门 操作有点多的裸树剖。套路下放边权后就是打板子。注意操作的时候不要包括LCA。 代码： #include&lt;cctype&gt; #include&lt;cstdio&gt; #include&lt;cstring&gt; #include&lt;algorithm&gt; #define N 20005 #define F inline ...
BZOJ 2157: 旅游( 树链剖分 )
weixin_30549175的博客
06-26 287
树链剖分..样例太大了根本没法调...顺便把数据生成器放上来--------------------------------------------------------------------------------#include<cstdio>#include<cstring>#include<algorithm>#include<iostrea...
【BZOJ】2157: 旅游
02-07 94
http://www.lydsy.com/JudgeOnline/problem.php?id=2157 题解：裸lct不解释.. #include <bits/stdc++.h> using namespace std; struct node *null; struct node { node *c[2], *f; bool flag, rev, t...
bzoj2157
weixin_34194702的博客
12-17 113
LCT板子，打个lazy即可 # include <stdio.h> # include <stdlib.h> # include <iostream> # include <algorithm> # include <string.h> # include <map> # define IL inline ...
BZOJ2157【Link Cut Tree】
Lethelody的专栏
03-21 430
/* I will wait for you */ #include #include #include #include #include #include #include #include #include #include #include #include #include #include typedef long long LL; typedef unsigned long lo
蓝桥杯-3516-景区导游
最新发布
qq_40345909的博客
04-24 699
某景区一共有 N 个景点，编号 1到 N。景点之间共有 N -1条双向的摆渡车线路相连，形成一棵树状结构。在景点之间往返只能通过这些摆渡车进行，需要花费一定的时间。小明是这个景区的资深导游，他每天都要按固定顺序带客人游览其中 K 个景点:41,42,...，AK。今天由于时间原因，小明决定跳过其中一个景点，只带游客按顺序游览其中 区-1个景点。具体来说，如果小明选择跳过 4i，那么他会按顺序带游客游览请你对任意一个 Ai，计算如果跳过这个景点，小明需要花费多少时。
[BZOJ5212][Zjoi2018]历史（树形dp+LCT）
女装的你，如此好看！
09-15 272
Address https://www.luogu.org/problemnew/show/P4338 https://www.lydsy.com/JudgeOnline/problem.php?id=5212 http://uoj.ac/problem/374 https://loj.ac/problem/2434 Solution 简版题意：给定树上每个点 Access 的次数...
[BZOJ3651]网络通信（LCT）
zyf2000
03-28 957
人之所以痛苦，在于追求错误的东西。如果你不给自己烦恼，别人也永远不可能给你烦恼。因为你自己的内心，你放不下。 好好的管教你自己，不要管别人。——送给我自己
bzoj4712: 洪水 动态Dp 树链剖分+线段树 或 LCT维护矩阵乘法
maxtir的博客
07-22 568
bzoj4712: 洪水 Description 小A走到一个山脚下，准备给自己造一个小屋。这时候，小A的朋友（op，又叫管理员）打开了创造模式，然后飞到 山顶放了格水。于是小A面前出现了一个瀑布。作为平民的小A只好老实巴交地爬山堵水。那么问题来了：我们把这 个瀑布看成是一个n个节点的树，每个节点有权值（爬上去的代价）。小A要选择一些节点，以其权值和作为代价将 这些点删除（堵上），使得根...
JAVA程序设计：追逐游戏（LCP 21）
信仰.的博客
09-14 1042
秋游中的小力和小扣设计了一个追逐游戏。他们选了秋日市集景区中的 N 个景点，景点编号为 1~N。此外，他们还选择了 N 条小路，满足任意两个景点之间都可以通过小路互相到达，且不存在两条连接景点相同的小路。整个游戏场景可视作一个无向连通图，记作二维数组 edges，数组中以 [a,b] 形式表示景点 a 与景点 b 之间有一条小路连通。 小力和小扣只能沿景点间的小路移动。小力的目标是在最快时间内追到小扣，小扣的目标是尽可能延后被小力追到的时间。游戏开始前，两人分别站在两个不同的景点 startA 和 sta
题目-负责任的导游
poison
12-07 839
负责任的导游 小红在一家旅游公司做导游。她当前的任务是带一些游客去著名的A景区游览。由于A景区的景点较多，且比较分散，一些景点之间修有双向道路，景区内开通了环保巴士在这些景点间运行。出于安全考虑，环保巴士行驶在不同的道路上都规定了自己的最大乘客数目。小红有一份包含景区景点间道路状况和环保巴士最大载客容量的地图。由于这次小红所带的团体较大，她无法用一辆巴士将全体团员带往目的地景点。虽然景点循环运行
66. Plus One
猫敷雪
01-30 473
题目描述(简单难度) 数组代表一个数字，[ 1, 2, 3 ] 就代表 123，然后给它加上 1，输出新的数组。数组每个位置只保存 1 位，也就是 0 到 9。 解法一 迭代 public class Plus_One { public static int[] plusOne(int[] digits) { //从最低位遍历 for (int i = digits....
ICPC NEAU Programming Contest 2020~~ D. 旅游
Long_hen的博客
06-08 278
皮皮准备去旅游，共有n个景点可以选择，景点编号为1~n，每个景点都有一个“美观值”ai​。共有m次查询，对于第i次查询，皮皮将从xi号景点开始游览，之后他会选择沿着编号递增的顺序选择游览其他景点，但是如果这个景点的美观值不大于他刚刚游览过的景点，他就会跳过这个景点。也就是说，皮皮在游览一个美观值为u的景点v后，他将游览的下一个景点是编号大于v、美观值大于u的，编号最小的景点。皮皮将一共访问yi​个景点，请你输出他最后一个访问的景点编号，如果他不能访问yi​个景点，输出-1。 输入描述 输入第一行一个整.
无向图最小环问题
Mintind的博客
02-08 1232
CODEVS 2611 观光旅游 题目描述 Description 某旅游区里面有Ｎ个景点。两个景点之间可能直接有道路相连，用a[i][j]表示它的长度，否则它们之间没有直接的道路相连。这里所说的道路是没有规定方向的，也就是说，如果从i到j有直接的道路，那么从j到i也有，并且长度与之相等。 旅游区规定：每个游客的旅游线路只能是一个回路（好霸道的规定）。也就是说，游客可以任取一个景点出发，依次经过若...
LeetCode(401-600)
HENRYyangtuo的博客
03-16 860
目录LeetCode 301. 删除无效的括号LeetCode 303. 区域和检索 - 数组不可变LeetCode 304. 二维区域和检索 - 矩阵不可变LeetCode 309. 最佳买卖股票时机含冷冻期LeetCode 322. 零钱兑换LeetCode 328. 奇偶链表LeetCode 336. 回文对LeetCode 344. 反转字符串LeetCode 367. 有效的完全平方数LeetCode 372. 超级次方LeetCode 377. 组合总和 ⅣLeetCode 383. 赎金信Le
【每日一题】8月3日题目精讲—小A的最短路
qq_35975367的博客
08-22 429
来源：牛客网： 文章目录题目描述题解：代码： 时间限制：C/C++ 3秒，其他语言6秒 空间限制：C/C++ 262144K，其他语言524288K 64bit IO Format: %lld 题目描述 小A这次来到一个景区去旅游，景区里面有N个景点，景点之间有N-1条路径。小A从当前的一个景点移动到下一个景点需要消耗一点的体力值。但是景区里面有两个景点比较特殊，它们之间是可以直接坐观光缆车通过，不需要消耗体力值。而小A不想走太多的路，所以他希望你能够告诉它，从当前的位置出发到他想要去的那个地方，他最少
BZOJ3999 [TJOI2015]旅游
neither_nor
12-29 1028
水题……链剖然后维护区间最大值，最小值，和区间从左往右走的答案，还有从右往左走的答案即可 #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #define
【蓝桥杯算法】旅游景点排序问题，配合案例讲解冒泡排序，选择排序，插入排序，快速排序
Java—wang的博客
03-31 993
问题描述: 小明来到一个景区游玩，但是他的时间有限，没办法逛遍所有的景点，所以他从网上下载了每个景点的评分，他希望能够使游览的景点的评分总和最高，于是他希望你帮他对于N个景点排下序。 输入格式: 输入的第一行包含一个正整数N，表示N个景点。 第二行有N个正整数，表示每个景点的评分值。 输出格式:　 　输出一行，包含N个正整数，表示N个景点的评分从大到小的排列 样例输入: 4 3 2 4 1 样例输出: 4 3 2 1 问题分析： 这个问题，其实就是一个简单的排序问题，之所以要写这个简单的题，是因为我
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值