【ZJOI2018】历史【结论】【LCT思想】

题意:一棵nnn个点的树,每个点有权值aia_iai,初始时给定。维护:

  1. 单点权值加上一个正数。

  2. 询问每个点恰好执行aia_iai次 access 操作,任意安排顺序的条件下切换轻重链总次数的最大值。

n≤4×105n\leq 4\times 10^5n4×105

一道奇怪的题(?)

对于一个点uuu,考虑它子树内两个相邻 access的点,它们在uuu这里切换了一次当且仅当它们属于 不同的 uuu的儿子 为根 的 子树。(为了方便,我们把uuu单独这一个点看成uuu的子树。)

把来自同一个子树内的一次access看成一种颜色的小球,我们希望相邻的不同色小球尽量多。

结论 设总小球数为nnn,颜色最多的小球个数为mmm,那么最大贡献为min⁡{n−1,2(n−m)}\min\{n-1,2(n-m)\}min{n1,2(nm)}

证明

考虑把最多的颜色摆出来,在m−1m-1m1个空隙中放其他颜色。

  1. 假设最多的颜色只有一种

把所有其他颜色的球依次往空隙中放,放完一个颜色继续往下一个空隙放,空隙放完了再从第一个开始。

如果空隙都放完了,即n−m≥m−1n-m\geq m-1nmm1,即n−1≤2(n−m)n-1\leq2(n-m)n12(nm),所有相邻的小球颜色都不同,答案为n−1n-1n1

否则颜色最多的球一定会相邻。考虑开始时贡献为000,插入一个任意颜色的球最多让贡献+2+2+2,上面的构造方法可以达到这个最大值,所以最大贡献为2(n−m)2(n-m)2(nm)

  1. 最多的颜色有多种

显然这时m≤n2m\leq \frac{n}{2}m2n,把最多的颜色各取一个构成一组绑在一起,然后其他颜色用同样的方法放在组之间的空隙内,这样不会有同色相邻。

就可以O(n)O(n)O(n)处理了


然后考虑怎么用数据结构维护这玩意

上面已经推过了,虽然表达式看起来很奇怪,但实际上本质上是判断最大的子树是否超过自己的一半(这里的大小都是指aaa的和)

一点也不自然地想到实链剖分

具体而言,设sumusum_usumu表示uuu子树内权值和,对于uuu和它的儿子vvv,如果2sumv>sumu+12sum_v>sum_u+12sumv>sumu+1(先不管uuu本身什么的),那么vvvuuu的重儿子,显然重儿子最多有一个。如果不存在这样的儿子就没有重儿子。

这样做的好处是重链上都是取的min⁡{n−1,2(n−m)}\min\{n-1,2(n-m)\}min{n1,2(nm)}的右边这项,而把父亲和儿子加上同一个正数时左边会变大,右边不会变,而右边本来就比左边小,不仅仍然是重边,连这个点的答案都不会改变,直接跳过去就可以了。而到根的轻边是O(log⁡V)O(\log V)O(logV)的,直接讨论一下轻边的情况即可。

具体而言就是魔改一下LCT的access,跳轻边的时候需要判一下。具体实现的时候很诡异,建议直接看代码。

注意复杂度是通过轻边条数而非finger-search保证的,所以必须在修改前dfs一遍处理出重儿子

复杂度应该是O(nlog⁡nlog⁡V)O(n\log n\log V)O(nlognlogV)

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
#include <vector>
#define MAXN 400005
using namespace std;
typedef long long ll;
inline int read()
{int ans=0;char c=getchar();while (!isdigit(c)) c=getchar();while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();return ans;
}
int ch[MAXN][2],fa[MAXN];
ll val[MAXN],lz[MAXN];
inline void pushlzy(int x,ll v){val[x]+=v,lz[x]+=v;}
inline void pushdown(int x)
{if (lz[x]){if (ch[x][0]) pushlzy(ch[x][0],lz[x]);if (ch[x][1]) pushlzy(ch[x][1],lz[x]);lz[x]=0;}
}
inline int get(int x){return ch[fa[x]][1]==x;}
inline bool isroot(int x){return ch[fa[x]][0]!=x&&ch[fa[x]][1]!=x;}
inline void rotate(int x)
{int y=fa[x],z=fa[y];int l=get(x),r=l^1;int w=ch[x][r];if (!isroot(y)) ch[z][get(y)]=x;ch[x][r]=y,ch[y][l]=w;if (w) fa[w]=y;fa[y]=x,fa[x]=z;
}
int q[MAXN],tp;
inline void splay(int x)
{q[tp=1]=x;for (int i=x;!isroot(i);i=fa[i]) q[++tp]=fa[i];for (int i=tp;i>=1;i--) pushdown(q[i]);while (!isroot(x)){int y=fa[x];if (!isroot(y)){if (get(x)==get(y)) rotate(y);else rotate(x);}rotate(x);}
}
ll res[MAXN],a[MAXN],ans;
inline int findrt(int x){while (ch[x][0]) pushdown(x),x=ch[x][0];return x;}
inline void access(int x,int v)
{a[x]+=v;for (int y=0;x;y=x,x=fa[x]){splay(x);val[x]+=v;if (ch[x][0]) pushlzy(ch[x][0],v);ans-=res[x];int t=findrt(ch[x][1]);if (2*(val[x]-val[t])>=val[x]-1) ch[x][1]=0;t=findrt(y);if (2*(val[x]-val[t])<val[x]-1) ch[x][1]=y;t=findrt(ch[x][1]);res[x]=(t? 2*(val[x]-val[t]):val[x]-1);if (2*a[x]>val[x]+1) res[x]=2*(val[x]-a[x]);ans+=res[x];}
}
vector<int> e[MAXN];
void dfs(int u)
{val[u]=a[u];for (int i=0;i<(int)e[u].size();i++)if (e[u][i]!=fa[u]){fa[e[u][i]]=u;dfs(e[u][i]);val[u]+=val[e[u][i]];}for (int i=0;i<(int)e[u].size();i++)if (e[u][i]!=fa[u]&&2*val[e[u][i]]>val[u]+1)ch[u][1]=e[u][i];res[u]=(ch[u][1]? 2*(val[u]-val[ch[u][1]]):val[u]-1);if (2*a[u]>val[u]+1) res[u]=2*(val[u]-a[u]);ans+=res[u];
}
int main()
{int n,m;n=read(),m=read();for (int i=1;i<=n;i++) a[i]=read();for (int i=1;i<n;i++){int u,v;u=read(),v=read();e[u].push_back(v),e[v].push_back(u);}dfs(1);	printf("%lld\n",ans);while (m--){int x,v;x=read(),v=read();access(x,v);printf("%lld\n",ans);}return 0;
}

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

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

相关文章

2020 ICPC 上海 Sum of Log 数位dp + 状态剪枝

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 观察可以发现&#xff0c;由于i&j0i \And j0i&j0&#xff0c;所以log2(ij)log_2(ij)log2​(ij)表示的应该是ijijij的111的最高位。 一个显然的dpdpdp状态就出来了f[pos][flag1][fl…

从MVP到微软产品经理的几点心得

昨天&#xff0c;我应邀在微软2019年暑期实习生的入职培训活动上做了一个简短的分享 —— 以我自身的职业发展过程为例&#xff0c;给大家讲了一讲我是怎样从一个技术爱好者&#xff0c;到获得第一个MVP&#xff08;微软最有价值专家&#xff09;的荣誉奖项&#xff0c;以及后来…

【UOJ188】 Sanrd【类min_25筛】

题意&#xff1a;设f(i)f(i)f(i)表示iii的不严格次大质因子&#xff08;没有为000&#xff09;&#xff0c;求∑ilrf(i)\sum_{il}^rf(i)∑ilr​f(i) l≤r≤1011l\leq r\leq10^{11}l≤r≤1011 这种和质因数有关的奇奇怪怪的函数的前缀和可以试试魔改min_25筛 设 S(n,j)∑i2n[m…

小白开学Asp.Net Core 《六》 —— 探究.Net Core 跨平台的奥秘

1、写这篇文章的初衷有好多朋友反馈看不懂我写的开源的一个练手项目&#xff08;GitHub&#xff1a;https://github.com/AjuPrince/Aju.Carefree&#xff09;也有好多人都希望我能写一些简单的入门的文章&#xff0c;记得前几天在群里有人问为什么 .Net Core 能跨平台&#xff…

Codeforces Round #636 (Div. 3) D. Constant Palindrome Sum 思维 + 差分

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 首先有一个显然的性质就是每组操作最多不会超过两次。 很容易想到一个很暴力的思路&#xff0c;就是枚举x∈[1,2∗k]x \in [1,2*k]x∈[1,2∗k]&#xff0c;让后判断一下每组需要操作几次取…

C#8.0: 在 LINQ 中支持异步的 IAsyncEnumerableT接口

C# 8.0中&#xff0c;提供了一种新的IAsyncEnumerable<T>接口&#xff0c;在对集合进行迭代时&#xff0c;支持异步操作。比如在读取文本中的多行字符串时&#xff0c;如果读取每行字符串的时候使用同步方法&#xff0c;那么会导致线程堵塞。IAsyncEnumerable<T>可…

SCOI2020游记

考试期间发现只有T1看起来可做&#xff0c;脑子发麻硬刚5小时&#xff0c;最后还是只交了40分暴力 考完后用一节数学课口胡出了T2 100pts和T3 80pts&#xff0c;并且发现只有T1不会做 退役了&#xff0c;告辞

Codeforces Round #636 (Div. 3) E. Weights Distributing 思维 + bfs

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; n≤2e5,m≤2e5n\le2e5,m\le2e5n≤2e5,m≤2e5 思路&#xff1a; 怎么感觉每场div3div3div3都有一个巧妙的图论题。 首先如果只有两个点的话&#xff0c;肯定是一次bfsbfsbfs之后取前disdisdis小的www作为边权…

解读 Microsoft.NET.Sdk 的源码,你能定制各种奇怪而富有创意的编译过程

在 csproj 中&#xff0c;Project 中的 Sdk 属性是 MSBuild 15.0 开始支持的&#xff0c;也就是 Visual Studio 2017 开始支持。有了 Sdk 属性的存在&#xff0c;MSBuild 编译过程能够扩展得非常强大&#xff0c;而不止是过去 Import 的一个 props 和 targets 文件。本文将介绍…

【省选2020A卷】作业题【矩阵树】【扩域】【莫比乌斯反演】

传送门 为什么世界上会有这么傻的题啊……我佛了 很显然就是矩阵树强行缝合莫反 设f(n)f(n)f(n)表示所有边权gcd⁡\gcdgcd为nnn的生成树权值和&#xff0c;g(d)g(d)g(d)表示所有边权都是 ddd的倍数的生成树权值和 g(d)∑d∣nf(n)g(d)\sum_{d\mid n}f(n)g(d)d∣n∑​f(n) f(d…

Libre OJ 「BalticOI 2013」非回文数 数位dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 问[l,r][l,r][l,r]内有多少个数是非回文数&#xff0c;即数字中不存在连续几个数为回文数。 l,r≤1e18l,r\le1e18l,r≤1e18 思路&#xff1a; 这么大的范围很明显数位dpdpdp了&#xff0c;容易知道当一个数…

《ASP.NET Core 微服务实战》译者序

最近&#xff0c;我将《ASP.NET Core 微服务实战》一书由英文翻译为中文。这本书是由清华大学出版社引进的&#xff0c;目前还处于最后的排版校对过程中&#xff0c;现将该书的译者序发表于此。以下为译者译全文&#xff1a;“微服务”的概念在 2014 年正式提出之后&#xff0c…

【杭电多校2020】Distinct Sub-palindromes【结论】【构造】

题意&#xff1a;长度为nnn的 本质不同的回文子串个数最少 的小写字母组成的字符串 的个数。 n≤109n\leq 10^9n≤109 回文子串个数至少为不同字母的数量 显然 最少回文子串个数为min⁡(n,3)\min(n,3)min(n,3) n≤2n\leq2n≤2显然 n≥3n\geq3n≥3时 构造abcabcabc...\textt…

.NET 云原生技术使用调查

.NET Core 3.0 即将发布RC&#xff0c;.NET Core 3.0默认更好的支持Docker资源限制&#xff0c;官方团队也在努力让.NET Core成为真正的容器运行时&#xff0c;使其在低内存环境中具有容器感知功能并高效运行。GC堆限制.NET Core减少了CoreCLR默认使用的内存&#xff0c;如G0代…

【杭电多校2020】Minimum Index【Lyndon Word】

题意&#xff1a;给定字符串&#xff0c;求所有前缀的最小后缀。 n≤2107n\leq 2\times10^7n≤2107 最小后缀就是Lyndon分解的最后一段。而Duval本质上是可以重复修改的增量算法&#xff0c;所以是可以做的。 记ansians_iansi​为前缀iii的最小后缀。设维护未确定的循环节的指…

FZU - 2042 The Mad Mathematician 数位dp + 算贡献

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 求sumsumsum。 a,b,c,d,e≤1e18a,b,c,d,e\le1e18a,b,c,d,e≤1e18 思路&#xff1a; 这是一篇无从考究的题解&#xff0c;因为fzu现在进不去。 看到这种题直接考虑数位dpdpdp&#xff0c;对于[A,B],[C,D][…

「PowerBI」分析服务多维数据结构重回关系数据库的一大绝招

在过往Excel催化剂定位的轻量级ExcelPowerBIDesktop的解决方案中&#xff0c;已经做过了近乎完美的PowerBIDesktop数据模型数据导出到Excel工作表的应用&#xff0c;这也是个人版数据应用的重要的刚需功能&#xff0c;由Excel催化剂免费奉献给PowerBI社区作完结。此篇再次带领大…

【杭电多校2020】Fibonacci Sum【斐波拉契通项】【推式子】

题意&#xff1a;设FiF_iFi​为斐波拉契数列&#xff0c;求 ∑i0N(FiC)k\sum_{i0}^N(F_{iC})^ki0∑N​(FiC​)k 模109910^991099 N,C≤1018,k≤105N,C\leq10^{18},k\leq10^5N,C≤1018,k≤105 把斐波拉契暴力拆开 FiC15[(152)iC−(1−52)iC]F_{iC}\frac 1 {\sqrt{5}}[(\frac {…

Codeforces Round #257 (Div. 1) D. Jzzhu and Numbers 高维前缀和 + 容斥

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 完全想不到容斥啊&#xff0c;看了半天也没看懂渍渍渍。 定义f[i]f[i]f[i]表示iii的超集个数&#xff0c;那么选择的方案就是2f[i]−12^{f[i]}-12f[i]−1了&#xff0c;因为不能一个不选所…

花5分钟时间来了解一下高性能网关Kong会有意外收获

前言前几天开源发布了 Kong.Net 项目&#xff0c;收到了大量园友的反馈&#xff0c;开源当天就突破了 100 个star &#xff0c;可喜可贺&#xff0c;但是从侧面也说明&#xff0c;我们 .NetCore 阵营真的非常需要拥抱开源&#xff0c;应该敞开心扉&#xff0c;集众家之长&#…