P6845 [CEOI2019] Dynamic Diameter
P6845 [CEOI2019] Dynamic Diameter
题意
一颗带权树,每次更改一条边的权,每次修改后求出最大直径。强制在线。
思路
\(O(n\log^2n)\) 的暴力做法。
根据经典结论,对于两个点集的树上最大直径(权值非负),并集点集的树上最大直径的端点一定是原四个端点中的两个。
那么我们就可以用线段树维护点集,合并时 \(O(\log n)\) 查询两点间距离合并就可以做到 \(O(n\log^2n)\) 的复杂度了。
思考如何支持在线修改边权。两点间距离为 \(dis(x)+dis(y)-2*dis(lca(x,y))\) 那么先树剖发现修改一条边的权实际上是将子树内 \(dis\) 增加,那么子树就是区间加。我们可以用树状数组维护修改。
考虑哪些部分的线段树上的点集最大直径被改了。发现子树内和子树外的最大直径不会变化,那么我们只需要更改子树内与外之间合并的最大直径就行了。我们按照 \(dfn\) 用线段树进行维护,那么我们每次只需要将子树区间重新 pushup
一下就可以了。
代码
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cctype>
#include<cstdlib>
#include<utility>
using namespace std;
inline long long read(){
long long x=0;int w=0;char c=getchar();
while(!isdigit(c)) w|=c=='-',c=getchar();
while(isdigit(c)) x=x*10+(c^48),c=getchar();
return w?-x:x;
}
namespace star
{
const int maxn=1e5+10;
int n,q;
int ecnt=1,head[maxn],to[maxn<<1],nxt[maxn<<1],go[maxn];
long long W,w[maxn],dis[maxn];
inline void addedge(){
int a=read(),b=read();
to[++ecnt]=b,nxt[ecnt]=head[a],head[a]=ecnt;
to[++ecnt]=a,nxt[ecnt]=head[b],head[b]=ecnt;
w[ecnt>>1]=read();
}
int dep[maxn],fa[maxn],top[maxn],siz[maxn],son[maxn],dfn[maxn],id[maxn];
void dfs1(int x,int f){
fa[x]=f,dep[x]=dep[f]+1,siz[x]=1;
for(int u,i=head[x];i;i=nxt[i]) if((u=to[i])!=f){
dis[u]=dis[x]+w[i>>1],dfs1(u,x),go[i>>1]=u;
if(siz[u]>siz[son[x]]) son[x]=u;
siz[x]+=siz[u];
}
}
void dfs2(int x,int topf){
top[x]=topf,dfn[x]=++dfn[0],id[dfn[0]]=x;
if(!son[x]) return;
dfs2(son[x],topf);
for(int u,i=head[x];i;i=nxt[i]) if((u=to[i])!=fa[x] and u!=son[x]) dfs2(u,u);
}
inline int LCA(int x,int y){
while(top[x]!=top[y]) dep[top[x]]>dep[top[y]]?(x=fa[top[x]]):(y=fa[top[y]]);
return dep[x]<dep[y]?x:y;
}
#define ls (ro<<1)
#define rs (ro<<1|1)
#define mid ((l+r)>>1)
namespace S{
long long c[maxn];
inline void upd(int x,long long k){for(;x<=n;x+=x&-x) c[x]+=k;}
inline long long query(int x){long long ans=dis[id[x]];for(;x;x-=x&-x) ans+=c[x];return ans;}
inline void update(int x,int y,long long w){upd(x,w),upd(y+1,-w);}
}
inline long long Dis(pair<int,int> a){return S::query(dfn[a.first])+S::query(dfn[a.second])-2*S::query(dfn[LCA(a.first,a.second)]);}
namespace T{
pair<int,int> e[maxn<<2];
inline pair<int,int> merge(const pair<int,int>& a,const pair<int,int>& b){
pair<int,int> p[6]={a,b,make_pair(a.first,b.first),make_pair(a.first,b.second),make_pair(a.second,b.first),make_pair(a.second,b.second)};
long long dis[6];
for(int i=0;i<6;i++) dis[i]=Dis(p[i]);
return p[max_element(dis,dis+6)-dis];
}
void build(int ro=1,int l=1,int r=n){
if(l==r) return e[ro]=make_pair(id[l],id[l]),void();
build(ls,l,mid),build(rs,mid+1,r);
e[ro]=merge(e[ls],e[rs]);
}
void update(int x,int y,int ro=1,int l=1,int r=n){
if(x==l and y==r) return;
if(y<=mid) update(x,y,ls,l,mid);
else if(x>mid) update(x,y,rs,mid+1,r);
else update(x,mid,ls,l,mid),update(mid+1,y,rs,mid+1,r);
e[ro]=merge(e[ls],e[rs]);
}
}
#undef ls
#undef rs
#undef mid
inline void work(){
n=read(),q=read(),W=read();
for(int i=1;i<n;i++) addedge();
dfs1(1,0),dfs2(1,1);
T::build();
long long ans=0;
while(q--){
int e=(read()+ans)%(n-1)+1;
long long v=(read()+ans)%W;
S::update(dfn[go[e]],dfn[go[e]]+siz[go[e]]-1,v-w[e]),w[e]=v;
T::update(dfn[go[e]],dfn[go[e]]+siz[go[e]]-1);
printf("%lld\n",ans=Dis(T::e[1]));
}
}
}
signed main(){
star::work();
return 0;
}
P6845 [CEOI2019] Dynamic Diameter的更多相关文章
- CEOI2019 / CodeForces 1192B. Dynamic Diameter
题目简述:给定一棵$N \leq 10^5$个节点的树,边上带权,维护以下两个操作: 1. 修改一条边的边权: 2. 询问当前树的直径长度. 解1:code 注意到树的直径有以下性质: 定理:令$\t ...
- NOIP前做题记录
鉴于某些原因(主要是懒)就不一题一题写了,代码直接去\(OJ\)上看吧 CodeChef Making Change 传送门 完全没看懂题解在讲什么(一定是因为题解公式打崩的原因才不是曲明英语太差呢- ...
- Diameter协议摘要
---------选择同学整理文档 1. 协议概述 Diameter协议主要为应用程序提供认证.鉴权.计费框架,即AAA,并支持本地AAA和漫游场景下的AAA. 1.1. 特点介绍 以前的AAA ...
- var和dynamic的区别
1.var 1.均是声明动态类型的变量. 2.在编译阶段已经确定类型,在初始化的时候必须提供初始化的值. 3.无法作为方法参数类型,也无法作为返回值类型. 2.dynamic 1.均是声明动态类型的变 ...
- 遍历dynamic的方式
一.遍历ExpandoObject /// <summary> /// 遍历ExpandoObject /// </summary> [TestMethod] public v ...
- C# dynamic 动态创建 json
1. 如何通过C# 的dynamic 创建如下json 对象? { "query": { "match": [{ "name": " ...
- BZOJ 1901: Zju2112 Dynamic Rankings[带修改的主席树]【学习笔记】
1901: Zju2112 Dynamic Rankings Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 7143 Solved: 2968[Su ...
- 当类型为dynamic的视图模型遭遇匿名对象
当年在ASP.NET MVC 1.0时代我提到,在开发时最好将视图的Model定制为强类型的,这样可以充分利用静态检查功能进行排错.不过有人指出,这么做虽然易于静态检查,但是定义强类型的Model类型 ...
- 动态规划 Dynamic Programming
March 26, 2013 作者:Hawstein 出处:http://hawstein.com/posts/dp-novice-to-advanced.html 声明:本文采用以下协议进行授权: ...
随机推荐
- JUC 并发编程--02,生产者和消费者 synchronized的写法 , juc的写法. Condition的用法
synchronized的写法 class PCdemo{ public static void main(String[] args) { //多个线程操作同一资源 Data data = new ...
- mysql-8.0.16-winx64.zip安装教程详解
摘自:https://zhuanlan.zhihu.com/p/48531203 mysql-8.0.16-winx64.zip安装教程详解下载地址:https://dev.mysql.com/get ...
- Java基础_循环嵌套_打印乘法口诀、菱形,各种图形,计算二元一次和三元一次方程组_7
循环嵌套 打印乘法口诀 for(int j=1;j<=9;j++){ for(int i=1;i<=j;i++){ System.out.print(i+"*"+j+& ...
- 二、部署监控服务器-Zabbix Server
二.部署监控服务器-Zabbix Server 1)源码安装Zabbix Server 多数源码包都是需要依赖包的,zabbix也- 样,源码编译前需要先安装相关依赖包. [root@zabbixse ...
- 2.docker下centos镜像
1.下载并运行 # 交互模式下载并运行centos容器 $ docker run -it centos:latest /bin/bash 1.1 配置centos的环境别名 $ vi /etc/bas ...
- 深度解密:Java与线程的关系
并发不一定要依赖多线程(如PHP的多进程并发),但在Java中谈论并发,大多数都与线程脱不开关系. 线程的实现 线程是CPU调度的基本单位,Thread类与大部分的Java API有显著的差别,它的所 ...
- 题解 P3940 分组
有些梦想虽然遥不可及,但不是不可能实现.只要我足够的强. 前言 调了挺长时间的,并查集合并的时候需要 find 一下,不然会炸内存.... 解题思路 参考了题解区一篇思路非常好的题解,在这里讲一下自己 ...
- 二、RabbitMQ 进阶特性及使用场景 [.NET]
前言 经过上一篇的介绍,相信大家对RabbitMQ 的各种概念有了一定的了解,及如何使用RabbitMQ.Client 去发送和消费消息. 特性及使用场景 1. TTL 过期时间 TTL可以用来指定q ...
- 对volatile的理解--从JMM以及单例模式剖析
请谈谈你对volatile的理解 1.volitale是Java虚拟机提供的一种轻量级的同步机制 三大特性1.1保证可见性 1.2不保证原子性 1.3禁止指令重排 首先保证可见性 1.1 可见性 概念 ...
- 垃圾处理器-CMS
一.简介 CMS垃圾收集器是一款用于老年代的,使用复制-清除-整理算法的垃圾收集器. 二.GC阶段 1.初始化标记(STW) 暂停应用程序线程,遍历 GC ROOTS 直接可达的对象并将其压入标记栈( ...