LG1505 [国家集训队]旅游 树链剖分
问题描述
题解
边化点权。
超级多操作的树剖板子...
以后就拿这个当树剖板子复习吧...
\(\mathrm{Code}\)
#include<bits/stdc++.h>
using namespace std;
template <typename Tp>
void read(Tp &x){
x=0;char ch=1;int fh;
while(ch!='-'&&(ch>'9'||ch<'0')) ch=getchar();
if(ch=='-') ch=getchar(),fh=-1;
else fh=1;
while(ch>='0'&&ch<='9') x=(x<<1)+(x<<3)+ch-'0',ch=getchar();
x*=fh;
}
template <typename Tp>
void fr(Tp &x){
char ch=1;
while(ch!='C'&&ch!='N'&&ch!='S'&&ch!='M') ch=getchar();
if(ch=='M'){
ch=1;while(ch!='A'&&ch!='I') ch=getchar();
}
if(ch=='C') x=1;
else if(ch=='N') x=2;
else if(ch=='S') x=3;
else if(ch=='A') x=4;
else x=5;
}
const int maxn=40007;
const int INF=0x3f3f3f3f;
int n;
int Head[maxn],to[maxn],Next[maxn],tot=1,w[maxn];
int size[maxn],fa[maxn],dep[maxn];
int New[maxn],pre[maxn],ind,top[maxn];
int son[maxn];
//edge id i -> 2*i-1,2*i
void add(int x,int y,int z){
to[++tot]=y,Next[tot]=Head[x],Head[x]=tot,w[tot]=z;
}
void dfs1(int x,int f,int dp){
size[x]=1,fa[x]=f,dep[x]=dp;
int mxx=-1;
for(int i=Head[x];i;i=Next[i]){
int y=to[i];
if(y==f) continue;
dfs1(y,x,dp+1);
size[x]+=size[y];
if(size[y]>mxx) mxx=size[y],son[x]=y;
}
}
void dfs2(int x,int tp){
New[x]=++ind,pre[ind]=x,top[x]=tp;
if(!son[x]) return;
dfs2(son[x],tp);
for(int i=Head[x];i;i=Next[i]){
int y=to[i];
if(y==fa[x]||y==son[x]) continue;
dfs2(y,y);
}
}
int val[maxn],f[maxn];
void preprocess(){
for(int i=2;i<=tot;i+=2){//错误笔记:tot一开始等于1,则第一条边从2开始
int x=to[i xor 1],y=to[i];
if(dep[x]>dep[y]) val[x]=w[i],f[i]=f[i xor 1]=x;//错误笔记:边化点权,要记录每个边化到哪个点上了。
else val[y]=w[i],f[i]=f[i xor 1]=y;
}
}
#define lfc (x<<1)
#define rgc ((x<<1)|1)
#define mid ((l+r)>>1)
int mx[maxn<<1],mi[maxn<<1],sum[maxn<<1];
int lazy[maxn<<1];
void pushup(int x){
mx[x]=max(mx[lfc],mx[rgc]);
mi[x]=min(mi[lfc],mi[rgc]);
sum[x]=sum[lfc]+sum[rgc];
}
void pushdown(int x){
if(!lazy[x]) return;
lazy[lfc]^=1,lazy[rgc]^=1;
sum[lfc]=-sum[lfc],sum[rgc]=-sum[rgc];
int lf=mx[lfc],rg=mx[rgc];
mx[lfc]=-mi[lfc],mx[rgc]=-mi[rgc];
mi[lfc]=-lf,mi[rgc]=-rg;
lazy[x]=0;//错误笔记:lazy标记有没有清空
}
void build(int x,int l,int r){
if(l==r){
mx[x]=mi[x]=sum[x]=val[pre[l]];
return;
}
build(lfc,l,mid);build(rgc,mid+1,r);
pushup(x);
}
int L,R,need;
void dd_change(int x,int l,int r){
if(l==r&&L==l){
sum[x]=mx[x]=mi[x]=need;return;
}
pushdown(x);
if(L<=mid) dd_change(lfc,l,mid);
else dd_change(rgc,mid+1,r);
pushup(x);
}
void qj_change(int x,int l,int r){
if(l>R||r<L) return;
if(L<=l&&r<=R){
lazy[x]^=1;int temp=mx[x];
sum[x]=-sum[x],mx[x]=-mi[x],mi[x]=-temp;
return;
}
pushdown(x);
qj_change(lfc,l,mid);qj_change(rgc,mid+1,r);
pushup(x);
}
int sum_query(int x,int l,int r){
if(l>R||r<L) return 0;
if(L<=l&&r<=R) return sum[x];
pushdown(x);
return sum_query(lfc,l,mid)+sum_query(rgc,mid+1,r);
}
int max_query(int x,int l,int r){
if(l>R||r<L) return -INF;
if(L<=l&&r<=R) return mx[x];
pushdown(x);
return max(max_query(lfc,l,mid),max_query(rgc,mid+1,r));
}
int min_query(int x,int l,int r){
if(l>R||r<L) return INF;
if(L<=l&&r<=R) return mi[x];
pushdown(x);
return min(min_query(lfc,l,mid),min_query(rgc,mid+1,r));
}
int T,op;
void cz1(){
int x;
read(x);read(need);
x=New[f[x<<1]],L=x;//错误笔记:双向存边,在数组中编号要*2
dd_change(1,1,n);
}
void cz2(){
int x,y;
read(x);read(y);
++x,++y;
while(top[x]!=top[y]){
if(dep[top[x]]<dep[top[y]]) swap(x,y);
L=New[top[x]],R=New[x];//错误笔记:写作L=New[x],R=New[y]
qj_change(1,1,n);
x=fa[top[x]];
}
if(dep[x]>dep[y]) swap(x,y);
L=New[x]+1,R=New[y];
if(L<=R) qj_change(1,1,n);
}
void cz3(){
int x,y,res=0;
read(x);read(y);
++x,++y;
while(top[x]!=top[y]){
if(dep[top[x]]<dep[top[y]]) swap(x,y);
L=New[top[x]],R=New[x];
res+=sum_query(1,1,n);
x=fa[top[x]];
}
if(dep[x]>dep[y]) swap(x,y);
L=New[x]+1,R=New[y];
if(L<=R) res+=sum_query(1,1,n);
printf("%d\n",res);
}
void cz4(){
int x,y,res=-INF;
read(x);read(y);
++x,++y;
while(top[x]!=top[y]){
if(dep[top[x]]<dep[top[y]]) swap(x,y);
L=New[top[x]],R=New[x];
res=max(res,max_query(1,1,n));
x=fa[top[x]];
}
if(dep[x]>dep[y]) swap(x,y);
L=New[x]+1,R=New[y];
if(L<=R) res=max(res,max_query(1,1,n));
printf("%d\n",res);
}
void cz5(){
int x,y,res=INF;
read(x);read(y);
++x,++y;
while(top[x]!=top[y]){
if(dep[top[x]]<dep[top[y]]) swap(x,y);
L=New[top[x]],R=New[x];
res=min(res,min_query(1,1,n));
x=fa[top[x]];
}
if(dep[x]>dep[y]) swap(x,y);
L=New[x]+1,R=New[y];
if(L<=R) res=min(res,min_query(1,1,n));
printf("%d\n",res);
}
int main(){
#ifndef ONLINE_JUDGE
freopen("travel.in","r",stdin);
#endif
read(n);
for(int i=1,x,y,z;i<n;i++){
read(x);read(y);read(z);
++x,++y;
add(x,y,z);add(y,x,z);
}
dfs1(1,0,1);dfs2(1,1);
preprocess();build(1,1,n);
read(T);
while(T--){
fr(op);
if(op==1) cz1();
else if(op==2) cz2();
else if(op==3) cz3();
else if(op==4) cz4();
else cz5();
}
return 0;
}
LG1505 [国家集训队]旅游 树链剖分的更多相关文章
- 洛谷 P1505 [国家集训队]旅游 树链剖分
目录 题面 题目链接 题目描述 输入输出格式 输入格式 输出格式 输入输出样例 输入样例: 输出样例: 说明 思路 AC代码 总结 题面 题目链接 P1505 [国家集训队]旅游 题目描述 Ray 乐 ...
- LUOGU P1505 [国家集训队]旅游 (树链剖分+线段树)
传送门 解题思路 快被调死的码农题,,,其实就是一个边权下放到点权的线段树+树剖. #include<iostream> #include<cstdio> #include&l ...
- BZOJ 2157: 旅游( 树链剖分 )
树链剖分.. 样例太大了根本没法调...顺便把数据生成器放上来 -------------------------------------------------------------------- ...
- BZOJ2157旅游——树链剖分+线段树
题目描述 Ray 乐忠于旅游,这次他来到了T 城.T 城是一个水上城市,一共有 N 个景点,有些景点之间会用一座桥连接.为了方便游客到达每个景点但又为了节约成本,T 城的任意两个景点之间有且只有一条路 ...
- BZOJ3999:[TJOI2015]旅游(树链剖分)
Description 为了提高智商,ZJY准备去往一个新世界去旅游.这个世界的城市布局像一棵树.每两座城市之间只有一条路径可 以互达.每座城市都有一种宝石,有一定的价格.ZJY为了赚取最高利益,她会 ...
- BZOJ2157: 旅游 树链剖分 线段树
http://www.lydsy.com/JudgeOnline/problem.php?id=2157 在对树中数据进行改动的时候需要很多pushdown(具体操作见代码),不然会wa,大概原因 ...
- BZOJ 2157: 旅游 (树链剖分+线段树)
树链剖分后线段树维护区间最大最小值与和. 支持单点修改与区间取反. 直接写个区间取反标记就行了.线段树板题.(200行6000B+ 1A警告) #include <cstdio> #inc ...
- 【BZOJ2157】旅游 树链剖分+线段树
[BZOJ2157]旅游 Description Ray 乐忠于旅游,这次他来到了T 城.T 城是一个水上城市,一共有 N 个景点,有些景点之间会用一座桥连接.为了方便游客到达每个景点但又为了节约成本 ...
- P1505 [国家集训队]旅游[树剖]
题目描述 Ray 乐忠于旅游,这次他来到了T 城.T 城是一个水上城市,一共有 N 个景点,有些景点之间会用一座桥连接.为了方便游客到达每个景点但又为了节约成本,T 城的任意两个景点之间有且只有一条路 ...
随机推荐
- acwing 25. 剪绳子
习题地址 https://www.acwing.com/problem/content/description/24/ 题目描述 给你一根长度为 nn 绳子,请把绳子剪成 mm 段(mm.nn 都是整 ...
- AcWing 836. 合并集合
地址 https://www.acwing.com/problem/content/838/ 一共有n个数,编号是1~n,最开始每个数各自在一个集合中. 现在要进行m个操作,操作共有两种: “M a ...
- 算法问题实战策略 CHILDRENDAY
地址 https://algospot.com/judge/problem/read/CHILDRENDAY 题解 ac代码 #include <iostream> #include &l ...
- 发布Jar包到中央仓库
参考流程 https://blog.csdn.net/qq_36838191/article/details/81027586 备份还原 gpg-keys https://blog.rathena.c ...
- ROS kinetic + Realsens D435i + ORK + LINEMOD 物体识别
1. ORK 网址:https://wg-perception.github.io/object_recognition_core/ ORK (Object Recognition Kitchen) ...
- PHP 构造函数
在PHP5以前的版本中,构造函数的名称必须与类名相同,这种方法在PHP5中仍然可以使用,但现在已经很少有人用了.PHP5以及之后的版本,构造函数用__construct()方法来声明,这样做的好处是可 ...
- 21个Java Collections面试问答
Java Collections框架是Java编程语言的核心API之一. 这是Java面试问题的重要主题之一.在这里,我列出了一些重要的Java集合面试问题和解答,以帮助您进行面试.这直接来自我14年 ...
- three.js实现世界3d地图
概况如下: 1.THREE.Shape绘制世界地图平面地图: 2.THREE.ExtrudeGeometry将绘制的平面沿着Z轴拉伸,实现3d效果: 效果图如下: 预览地址:three.js实现世界3 ...
- idea使用心得
简单的概括如下: IntelliJ系中的 Project 相当于Eclipse系中的 Workspace : IntelliJ系中的 Module 相当于Eclipse系中的 Project ...
- Android集成极光推送
要说学习极光推送,个人感觉官方文档就非常好啦,但是没法,人太懒啦,为了下次能够快速的将极光推送集成到项目中,故结合之前开发的项目和官方文档记录下简单的Android集成极光推送,在这之前,先上一张简单 ...