XXI Open Cup. Grand Prix of Korea I. Query On A Tree 17 树剖 + 二分 + 树带权重心

传送门

文章目录

  • 题意:
  • 思路:

题意:

给你一棵树,每棵树初始权值都为000,现在给你两个操作:

(1)(1)(1)uuu的子树权值全部加111

(2)(2)(2)(u,v)(u,v)(u,v)路径上的点权值都加111

每次输出一个点xxx,满足∑y=1Na[y]∗dis(x,y)\sum_{y=1}^Na[y]*dis(x,y)y=1Na[y]dis(x,y)最小,如果多个点相同,输出深度最小的点。

思路:

观察这个式子,发现就是一个带权重心的问题,如果不在重心显然向重心走更优,需要支持修改。

这个题有一个结论:深度最小的带权重心的子树带权和一定>>>总和的一半。

用反证法来证明,假设这个不满足,即非这个子树的权值和,一定≥\ge权值和的一半,显然向其父亲节点走之后,会使得减少至少一半的权值,加上不到一半的权值,会更优。

考虑如何找到这样的一个点,显然我们可与从头dfsdfsdfs每个点,判断是否合法,这个复杂度O(n)O(n)O(n)

但是根据题面,我们显然要打一个树剖,这就跟dfsdfsdfs序扯上关系了,所以考虑能否将上面的问题转换到dfsdfsdfs上,用线段树来维护呢?

假设当前权值总和为sumsumsum,考虑在线段树上二分一个前缀≤sum2+1\le\frac{sum}{2}+12sum+1的最大位置,将dfsdfsdfs映射到点,假设为pospospos,此时满足了什么呢?满足pospospos子树的权值总和≤sum2+1\le \frac{sum}{2}+12sum+1,让后我们再向上倍增找到答案即可。

写完就过有被爽到。。

// Problem: I. Query On A Tree 17
// Contest: Codeforces - XXI Open Cup. Grand Prix of Korea
// URL: https://codeforces.com/gym/102759/problem/I
// Memory Limit: 1024 MB
// Time Limit: 2000 ms
// 
// Powered by CP Editor (https://cpeditor.org)//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#include<random>
#include<cassert>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid ((tr[u].l+tr[u].r)>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;const int N=1000010,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;int n;
int dfn[N],se[N],fa[N][20],top[N],son[N],tot,depth[N];
int inv[N];
LL sum;
vector<int>v[N];
struct Node {int l,r;LL sum,lazy;
}tr[N<<2];void dfs1(int u,int f) {fa[u][0]=f; depth[u]=depth[f]+1; se[u]=1;for(int i=1;i<=18;i++) fa[u][i]=fa[fa[u][i-1]][i-1];for(auto x:v[u]) {if(x==f) continue;dfs1(x,u);se[u]+=se[x];if(se[x]>se[son[u]]) son[u]=x;}
}void dfs2(int u,int t) {top[u]=t; dfn[u]=++tot;inv[tot]=u;if(son[u]) dfs2(son[u],t);for(auto x:v[u]) {if(x==fa[u][0]||x==son[u]) continue;dfs2(x,x);}
}void pushup(int u) {tr[u].sum=tr[L].sum+tr[R].sum;
}void pushdown(int u) {LL lazy=tr[u].lazy; tr[u].lazy=0;tr[L].sum+=Len(L)*lazy; tr[L].lazy+=lazy;tr[R].sum+=Len(R)*lazy; tr[R].lazy+=lazy;
}void build(int u,int l,int r) {tr[u]={l,r};if(l==r) return;build(L,l,Mid); build(R,Mid+1,r);
}void change(int u,int l,int r) {if(tr[u].l>=l&&tr[u].r<=r) {tr[u].sum+=Len(u);tr[u].lazy++;return;}pushdown(u);if(l<=Mid) change(L,l,r);if(r>Mid) change(R,l,r);pushup(u);
}LL query_sum(int u,int l,int r) {if(tr[u].l>=l&&tr[u].r<=r) return tr[u].sum;pushdown(u);LL ans=0;if(l<=Mid) ans+=query_sum(L,l,r);if(r>Mid) ans+=query_sum(R,l,r);return ans;
}int query_k(int u,LL sum) {if(tr[u].l==tr[u].r) return inv[tr[u].l];pushdown(u);if(sum<=tr[L].sum) return query_k(L,sum);else return query_k(R,sum-tr[L].sum);
}void update(int x,int y) {while(top[x]!=top[y]) {if(depth[top[x]]<depth[top[y]]) swap(x,y);change(1,dfn[top[x]],dfn[x]);x=fa[top[x]][0]; }	if(depth[x]>depth[y]) swap(x,y);change(1,dfn[x],dfn[y]); 
}int solve() {LL now=tr[1].sum/2+1;int pos=query_k(1,now);if(query_sum(1,dfn[pos],dfn[pos]+se[pos]-1)>=now) return pos;for(int i=18;i>=0;i--) if(fa[pos][i]&&query_sum(1,dfn[fa[pos][i]],dfn[fa[pos][i]]+se[fa[pos][i]]-1)<now) pos=fa[pos][i];return fa[pos][0];
}int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);	scanf("%d",&n);for(int i=1;i<=n-1;i++) {int a,b; scanf("%d%d",&a,&b);v[a].pb(b); v[b].pb(a);}dfs1(1,0); dfs2(1,1);build(1,1,n);int q; scanf("%d",&q);while(q--) {int op,l,r;scanf("%d%d",&op,&l);if(op==1) change(1,dfn[l],dfn[l]+se[l]-1);else scanf("%d",&r),update(l,r);printf("%d\n",solve());} return 0;
}
/**/

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/314631.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

使用Kubeadm创建k8s集群之节点部署(三十二)

前言 由于上次忘开申明原创&#xff0c;特再发一次。本篇部署教程将讲述k8s集群的节点&#xff08;master和工作节点&#xff09;部署&#xff0c;请先按照上一篇教程完成节点的准备。本篇教程中的操作全部使用脚本完成&#xff0c;并且对于某些情况&#xff08;比如镜…

Trie:hdu 4825、1251、1247、Poj 3764

hdu 4825链接 题目意思很简单&#xff0c;就是要求最大异或值的数。 我们可以从二进制的最高位开始选择&#xff0c;不断的排除一些数。我们先假设存在某些数字的二进制数是与当前查找的数不一样的&#xff0c;我们进入这一部分数进行查找&#xff0c;以此重复&#xff0c;不断…

Codeforces Round #675 (Div. 2) F. Boring Queries 区间lcm + 主席树

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个长度为nnn的序列aaa&#xff0c;qqq个询问&#xff0c;每次询问[l,r][l,r][l,r]内的lcmlcmlcm是多少&#xff0c;对1e971e971e97取模。 n≤1e5,a≤2e5,q≤1e5n\le1e5,a\le2e5,q\le1e5n≤1e5,a≤2e5,…

ASP.NET Core on K8S深入学习(2)部署过程解析与部署Dashboard

上一篇《K8S集群部署》中搭建好了一个最小化的K8S集群&#xff0c;这一篇我们来部署一个ASP.NET Core WebAPI项目来介绍一下整个部署过程的运行机制&#xff0c;然后部署一下Dashboard&#xff0c;完成可视化管理。本篇已加入了《.NET Core on K8S学习实践系列文章索引》&#…

字符Hash初步

兔子与兔子 很久很久以前&#xff0c;森林里住着一群兔子。 有一天&#xff0c;兔子们想要研究自己的 DNA 序列。 我们首先选取一个好长好长的 DNA 序列&#xff08;小兔子是外星生物&#xff0c;DNA 序列可能包含 26 个小写英文字母&#xff09;。 然后我们每次选择两个区间&…

02 | 健康之路 kubernetes(k8s) 实践之路 : 生产可用环境及验证

上一篇《 01 | 健康之路 kubernetes(k8s) 实践之路 : 开篇及概况 》我们介绍了我们的大体情况&#xff0c;也算迈出了第一步。今天我们主要介绍下我们生产可用的集群架设方案。涉及了整体拓补图&#xff0c;和我们采用的硬件配置&#xff0c;目前存在的问题等内容。遵循上一篇提…

NWERC 2018 C. Circuit Board Design 树 + 构造

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一颗nnn个点的树&#xff0c;让你在二维平面中构造一颗树&#xff0c;保证相邻点的距离正好为111&#xff0c;并且线段不能有相交&#xff0c;坐标绝对值≤3e3\le3e3≤3e3。 n≤1e3n\le1e3n≤1e3 思路&…

[开源] .NETCore websocket 即时通讯组件---ImCore

ImCore 是一款 .NETCore 下利用 WebSocket 实现的简易、高性能、集群即时通讯组件&#xff0c;支持点对点通讯、群聊通讯、上线下线事件消息等众多实用性功能。开源地址&#xff1a;https://github.com/2881099/im &#xff0c;求 star~~dotnet add package ImCoreIM服务端publ…

算法竞赛进阶指南——后缀数组

后缀数组 后缀数组 (SA) 是一种重要的数据结构&#xff0c;通常使用倍增或者DC3算法实现&#xff0c;这超出了我们的讨论范围。 在本题中&#xff0c;我们希望使用快排、Hash与二分实现一个简单的O(nlog2n)的后缀数组求法。 详细地说&#xff0c;给定一个长度为 n 的字符串S&a…

NWERC 2018 A. Access Points 二维转一维 + 单调栈

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你平面上nnn个点&#xff0c;你需要对于1−n1-n1−n依次选择nnn个点作为每个点的终点&#xff0c;满足选择的点i<j,xi≤xj,yi≤yji<j,x_i\le x_j,y_i\le y_ji<j,xi​≤xj​,yi​≤yj​&#xff0…

【译】在 Linux 上不安装 Mono 构建 .NET Framework 类库

在这篇文章中&#xff0c;我展示了如何在 Linux 上构建针对 .NET Framework 版本的.NET项目&#xff0c;而不使用 Mono。通用使用微软新发布的 Mocrosoft.NETFramework.ReferenceAssemblies NuGet 包&#xff0c;您将不需要安装除 .NET Core SDK 之外的任何其他软件包&#xff…

hdu 5023 线段树染色问题

题目链接 A Corrupt Mayor’s Performance Art Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 100000/100000 K (Java/Others) Total Submission(s): 4094 Accepted Submission(s): 1418 Problem Description Corrupt governors always find ways to get dirty money…

Codeforces Round #740 (Div. 2) D2. Up the Strip dp + 分块优化 + 逆向思维

传送门 文章目录题意&#xff1a;思路题意&#xff1a; 有nnn个细胞&#xff0c;你初始在第nnn细胞上&#xff0c;假设你当前在xxx处&#xff0c;你每次可以进行如下两个操作&#xff1a; (1)(1)(1)选择[1,x−1][1,x-1][1,x−1]内一个数yyy&#xff0c;跳到第x−yx-yx−y个细胞…

poj 2528 线段树离散化+染色

题目链接 Mayor’s posters Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 92628 Accepted: 26452 Description The citizens of Bytetown, AB, could not stand that the candidates in the mayoral election campaign have been placing their electoral pos…

Codeforces Round #740 (Div. 2) E. Bottom-Tier Reversals 构造

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个长度为奇数nnn的排列aaa&#xff0c;每次可以选择长度为奇数的前缀&#xff0c;并将[1,len][1,len][1,len]翻转&#xff0c;你需要用不超过5n2\frac{5n}{2}25n​次操作将其变成有序的&#xff0c;输…

DevOps书单:调研了101名专家,推荐这39本必读书籍

任何一个领域都遵循从新人到熟手&#xff0c;从熟手到专家的路径。在成长过程中&#xff0c;DevOps人经常会陷入没人带&#xff0c;没人管&#xff0c;找不到职业方向的迷茫。DevOps是在商业演进与企业协作的进化过程中诞生的一个全新职业&#xff0c;被很多人看成是一个“全栈…

模板:Prime最小生成树堆优化 + Dijkstra单源最短路堆优化

Dijkstra 单源最短路堆优化 #include<bits/stdc.h> using namespace std; typedef pair<int, int> PII; const int N 2e5 10; int head[N], to[N], value[N], nex[N], cnt 1; int n, m, rt, st, dis[N], visit[N]; struct cmp {bool operator()(PII a, PII b) …

Gym - 102001K Boomerangs 构造 + 三元环

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一张nnn个点mmm条边的简单图&#xff0c;让你找出尽可能多的三元环&#xff0c;要求每个三元环都不能共边&#xff0c;输出三元环数量和具体是那个。 n,m≤1e5n,m\le1e5n,m≤1e5 思路&#xff1a; 其实…

面对人性,有的选择向左,有的向右

这里是Z哥的个人公众号每周五11&#xff1a;45 按时送达有时也会有感而发&#xff0c;来加个餐&#xff5e;我的第「84」篇原创敬上前天早上&#xff0c;在36kr看到一篇文章《一个负能量的人&#xff0c;可以轻易搞垮周围人的生活》&#xff0c;讲述了可能我们每个人身边都存在…

10分钟了解分布式CAP、BASE理论

CAP理论2000年7月&#xff0c;Eric Brewer教授提出CAP猜想&#xff1b;2年后&#xff0c;Seth Gilbert和Nancy Lynch从理论上证明了CAP&#xff1b;之后&#xff0c;CAP理论正式成为分布式计算领域的公认定理。CAP定律说的是在一个分布式计算机系统中&#xff0c;一致性&#x…