九条可怜是一个热爱阅读的女孩子。

这段时间,她看了一本非常有趣的小说,这本小说的架空世界引起了她的兴趣。

这个世界有 \(n\) 个城市,这 \(n\) 个城市被恰好 \(n-1\) 条双向道路联通,即任意两个城市都可以互相到达。同时城市 \(1\) 坐落在世界的中心,占领了这个城市就称霸了这个世界。

在最开始,这 \(n\) 个城市都不在任何国家的控制之下,但是随着社会的发展,一些城市会崛起形成国家并夺取世界的霸权。为了方便,我们标记第 \(i\) 个城市崛起产生的国家为第 \(i\) 个国家。在第 \(i\) 个城市崛起的过程中,第 \(i\) 个国家会取得城市 \(i\) 到城市 \(1\) 路径上所有城市的控制权。

新的城市的崛起往往意味着战争与死亡,若第 \(i\) 个国家在崛起中,需要取得一个原本被国家 \(j\)(\(j\ne i\))控制的城市的控制权,那么国家 \(i\) 就必须向国家 \(j\) 宣战并进行战争。

现在,可怜知道了,在历史上,第 \(i\) 个城市一共崛起了 \(a_i\) 次。但是这些事件发生的相对顺序已经无从考究了,唯一的信息是,在一个城市崛起称霸世界之前,新的城市是不会崛起的。

战争对人民来说是灾难性的。可怜定义一次崛起的灾难度为崛起的过程中会和多少不同的国家进行战争(和同一个国家进行多次战争只会被计入一次)。可怜想要知道,在所有可能的崛起顺序中,灾难度之和最大是多少。

同时,在考古学家的努力下,越来越多的历史资料被发掘了出来,根据这些新的资料,可怜会对 \(a_i\) 进行一些修正。具体来说,可怜会对 \(a_i\) 进行一些操作,每次会将 \(a_x\) 加上 \(w\)。她希望在每次修改之后,都能计算得到最大的灾难度。

然而可怜对复杂的计算并不感兴趣,因此她想让你来帮她计算一下这些数值。

对题面的一些补充:

  • 同一个城市多次崛起形成的国家是同一个国家,这意味着同一个城市连续崛起两次是不会和任何国家开战的:因为这些城市原来就在它的控制之下。
  • 在历史的演变过程中,第 $i$ 个国家可能会有一段时间没有任何城市的控制权。但是这并不意味着第 $i$ 个国家灭亡了,在城市 $i$ 崛起的时候,第 $i$ 个国家仍然会取得 $1$ 到 $i$ 路径上的城市的控制权。
  • 输入格式

    第一行输入两个整数 \(n,m\) 表示城市个数和操作个数。

    第二行输入 \(n\) 个整数表示 \(a_i\) 的初始值。

    接下来 \(n-1\) 行,每行输入两个整数 \(u_i,v_i\)(\(1\le u_i,v_i\le n\))描述了一条道路。

    接下来 \(m\) 行每行输入两个整数 \(x_i,w_i\) 表示将 \(a_{x_i}\) 加上 \(w_i\)。

    输出格式

    输出共 \(m+1\) 行,第一行表示初始的 \(a_i\) 的答案,接下来 \(m\) 行每行表示这次修正后的答案。

    样例一

    input

    5 3

    1 1 1 1 1

    1 2

    1 3

    2 4

    2 5

    2 1

    3 1

    4 1

    output

    6

    7

    9

    10

    explanation

    在修正开始之前,如果按照所在城市 \(4,1,5,3,2\) 的顺序崛起,那么依次会和 \(0,1,2,1,2\) 个国家进行战争。这时一共会产生 \(6\) 对敌对关系。可以证明这是所有崛起顺序中的最大值。

    样例二

    见样例数据下载。

    限制与约定

    测试点 $n$ $m$ 其他约定
    1 $\le 10$ $=0$ $a_i=1$
    2 $\le 150000$ $\le 150000$ 第 $i$ 条道路连接 $i$ 和 $i + 1$
    3
    4 $=0$
    5
    6 $\le 150000$
    7
    8
    9 $\le 4\times 10^5$ $\le 4\times 10^5$
    10

    对于 100% 的数据,\(1\le a_i,w_i\le 10^7\),\(1\le x_i\le n\)。

    时间限制:\(2\texttt{s}\)

    空间限制:\(512\texttt{MB}\)

    题解

    首先考虑静态询问答案,我们令 \(sum_i\) 表示 \(i\) 的子树的 \(a_v\) 之和,然后想要求这个点的儿子的子树中对当前点能产生的最大贡献

    正常情况下是:\(sum_i-1\),就是每次从不同儿子的子树中交错上来

    但是如果存在一个儿子 \(v\) 的 \(sum_v\) 大于 \(sum_i\) 的一半,那么我们就做不到每次交错地选,最后一定会有剩下的从同一子树上来,这样是无法产生贡献的(与同一个国家的战争只算一次贡献,这样的贡献在儿子处已经计算完了)。这样的情况下,答案是 \(2(sum_i-\max\{\max\{sum_v\},a_i\})\)

    这样我们就得到了一个 \(O(n)\) 的静态算法

    考虑动态的时候,我们发现修改一个点的 \(a_i\) 改变的就是这个点到根路径上的 \(sum_f\) 值,我们要对这些点进行修改。这不就想LCT的access吗,我们只要把重儿子的定义修改为存在上述第二种情况,然后每次修改的时候就可以重新计算贡献了

    #include<bits/stdc++.h>
    #define ui unsigned int
    #define ll long long
    #define db double
    #define ld long double
    #define ull unsigned long long
    #define ft first
    #define sd second
    #define pb(a) push_back(a)
    #define PII std::pair<int,int>
    #define PLL std::pair<ll,ll>
    #define mp(a,b) std::make_pair(a,b)
    #define ITR(a,b) for(auto a:b)
    #define REP(a,b,c) for(register int a=(b),a##end=(c);a<=a##end;++a)
    #define DEP(a,b,c) for(register int a=(b),a##end=(c);a>=a##end;--a)
    const int MAXN=400000+10;
    int n,m,hson[MAXN],e,beg[MAXN],nex[MAXN<<1],to[MAXN<<1];
    ll size[MAXN],a[MAXN],ans,val[MAXN];
    template<typename T> inline void read(T &x)
    {
    T data=0,w=1;
    char ch=0;
    while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
    if(ch=='-')w=-1,ch=getchar();
    while(ch>='0'&&ch<='9')data=((T)data<<3)+((T)data<<1)+(ch^'0'),ch=getchar();
    x=data*w;
    }
    template<typename T> inline void write(T x,char ch='\0')
    {
    if(x<0)putchar('-'),x=-x;
    if(x>9)write(x/10);
    putchar(x%10+'0');
    if(ch!='\0')putchar(ch);
    }
    template<typename T> inline bool chkmin(T &x,T y){return y<x?(x=y,true):false;}
    template<typename T> inline bool chkmax(T &x,T y){return y>x?(x=y,true):false;}
    template<typename T> inline T min(T x,T y){return x<y?x:y;}
    template<typename T> inline T max(T x,T y){return x>y?x:y;}
    inline void calc(int x)
    {
    ans-=val[x];
    val[x]=hson[x]?2*(size[x]-size[hson[x]]):size[x]-1;
    if(a[x]*2>size[x]+1)val[x]=2*(size[x]-a[x]);
    ans+=val[x];
    }
    #define lc(x) ch[(x)][0]
    #define rc(x) ch[(x)][1]
    struct LinkCut_Tree{
    int fa[MAXN],ch[MAXN][2],stk[MAXN],cnt;
    ll add[MAXN];
    inline bool nroot(int x)
    {
    return lc(fa[x])==x||rc(fa[x])==x;
    }
    inline void plus(int x,ll v)
    {
    size[x]+=v;add[x]+=v;
    }
    inline void pushdown(int x)
    {
    if(add[x])
    {
    if(lc(x))plus(lc(x),add[x]);
    if(rc(x))plus(rc(x),add[x]);
    add[x]=0;
    }
    }
    inline void rotate(int x)
    {
    int f=fa[x],p=fa[f],c=(rc(f)==x);
    if(nroot(f))ch[p][rc(p)==f]=x;
    fa[ch[f][c]=ch[x][c^1]]=f;
    fa[ch[x][c^1]=f]=x;
    fa[x]=p;
    }
    inline void splay(int x)
    {
    cnt=0;
    stk[++cnt]=x;
    for(register int i=x;nroot(i);i=fa[i])stk[++cnt]=fa[i];
    while(cnt)pushdown(stk[cnt--]);
    for(register int y=fa[x];nroot(x);rotate(x),y=fa[x])
    if(nroot(y))rotate((lc(y)==x)==(lc(fa[y])==y)?y:x);
    }
    inline int findroot(int x)
    {
    while(lc(x))pushdown(x),x=lc(x);
    return x;
    }
    inline void access(int x,int v)
    {
    for(register int y=0;x;x=fa[y=x])
    {
    splay(x);size[x]+=v;
    if(lc(x))plus(lc(x),v);
    if(hson[x])
    {
    cnt=0;
    stk[++cnt]=hson[x];
    for(register int i=hson[x];nroot(i);i=fa[i])stk[++cnt]=fa[i];
    while(cnt)pushdown(stk[cnt--]);
    if(size[hson[x]]*2<=size[x]+1)hson[x]=rc(x)=0;
    }
    int nrt=findroot(y);
    if(size[nrt]*2>size[x]+1)hson[x]=nrt,rc(x)=y;
    calc(x);
    }
    }
    };
    LinkCut_Tree T;
    #undef lc
    #undef rc
    inline void insert(int x,int y)
    {
    to[++e]=y;
    nex[e]=beg[x];
    beg[x]=e;
    }
    inline void dfs(int x,int f)
    {
    ll res=0;
    size[x]=a[x];T.fa[x]=f;
    for(register int i=beg[x];i;i=nex[i])
    if(to[i]==f)continue;
    else
    {
    dfs(to[i],x);
    size[x]+=size[to[i]];
    if(chkmax(res,size[to[i]]))hson[x]=to[i];
    }
    if(size[hson[x]]*2<=size[x])hson[x]=0;
    T.ch[x][1]=hson[x];
    calc(x);
    }
    int main()
    {
    read(n);read(m);
    REP(i,1,n)read(a[i]);
    REP(i,1,n-1)
    {
    int u,v;read(u);read(v);
    insert(u,v);insert(v,u);
    }
    dfs(1,0);
    printf("%lld\n",ans);
    while(m--)
    {
    int x,w;read(x);read(w);
    a[x]+=w;T.access(x,w);
    printf("%lld\n",ans);
    }
    return 0;
    }

    【刷题】UOJ #374 【ZJOI2018】历史的更多相关文章

    1. ZJOI2019一轮停课刷题记录

      Preface 菜鸡HL终于狗来了他的省选停课,这次的时间很长,暂定停到一试结束,不过有机会二试的话还是可以搞到4月了 这段时间的学习就变得量大而且杂了,一般以刷薄弱的知识点和补一些新的奇怪技巧为主. ...

    2. 【BZOJ5212】[ZJOI2018]历史(Link-Cut Tree)

      [BZOJ5212][ZJOI2018]历史(Link-Cut Tree) 题面 洛谷 BZOJ 题解 显然实际上就是给定了一棵树和每个点被\(access\)的次数,求解轻重链切换的最大次数. 先考 ...

    3. 洛谷P1926 小书童—刷题大军【01背包】

      题目链接:https://www.luogu.org/problemnew/show/P1926 题目背景 数学是火,点亮物理的灯:物理是灯,照亮化学的路:化学是路,通向生物的坑:生物是坑,埋葬学理的 ...

    4. 好像leeceode题目我的博客太长了,需要重新建立一个. leecode刷题第二个

      376. Wiggle Subsequence               自己没想出来,看了别人的分析. 主要是要分析出升序降序只跟临近的2个决定.虽然直觉上不是这样. 455. 分发饼干     ...

    5. OI刷题录——hahalidaxin

      16-3-25  —— bzoj 2049 [Sdoi2008]Cave 洞穴勘测:LCT入门 bzoj 2002 [Hnoi2010]Bounce 弹飞绵羊:LCT Tsinsen A1303. t ...

    6. 牛客网Java刷题知识点之Map的两种取值方式keySet和entrySet、HashMap 、Hashtable、TreeMap、LinkedHashMap、ConcurrentHashMap 、WeakHashMap

      不多说,直接上干货! 这篇我是从整体出发去写的. 牛客网Java刷题知识点之Java 集合框架的构成.集合框架中的迭代器Iterator.集合框架中的集合接口Collection(List和Set). ...

    7. LeetCode刷题感想

      断断续续用了半年的时间把LeetCode刷完了,之前复习了数据结构与算法.将刷题与复习数据结构结合起来会更有效果.总之不是为了刷题而刷题,而是为了巩固和补充一部分知识. LeetCode真的是一个很好 ...

    8. 洛谷 P1926 小书童——刷题大军

      题目背景 数学是火,点亮物理的灯:物理是灯,照亮化学的路:化学是路,通向生物的坑:生物是坑,埋葬学理的人. 文言是火,点亮历史宫灯:历史是灯,照亮社会之路:社会是路,通向哲学大坑:哲学是坑,埋葬文科生 ...

    9. AC日记——小书童——刷题大军 洛谷 P1926

      题目背景 数学是火,点亮物理的灯:物理是灯,照亮化学的路:化学是路,通向生物的坑:生物是坑,埋葬学理的人. 文言是火,点亮历史宫灯:历史是灯,照亮社会之路:社会是路,通向哲学大坑:哲学是坑,埋葬文科生 ...

    10. 洛谷—— P1926 小书童——刷题大军

      https://www.luogu.org/problem/show?pid=1926#sub 题目背景 数学是火,点亮物理的灯:物理是灯,照亮化学的路:化学是路,通向生物的坑:生物是坑,埋葬学理的人 ...

    随机推荐

    1. React-redux-saga

      新建sagas.js import { takeEvery , put} from 'redux-saga/effects' import axios from 'axios'; import { G ...

    2. [转]The superclass "javax.servlet.http.HttpServlet" was not found on the Java Build Path

      完整错误信息: THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS"AS IS" AND ANY ...

    3. Scala学习(六)---Scala对象

      Scala中的对象 摘要: 在本篇中,你将会学到何时使用Scala的object语法结构.在你需要某个类的单个实例时,或者想为其他值或函数找一个可以挂靠的地方时,你就会用到它.本篇的要点包括: 1. ...

    4. 【nodejs】让nodejs像后端mvc框架(asp.net mvc)一样处理请求--参数自动映射篇(6/8)

      文章目录 前情概要 路由.action的扫描.发现.注册搞定之后,后来我发现在我们的action里面获取参数往往都是通过request对象来一个一个获取.同样的一行代码我们不厌其烦的重复写了无数次.遂 ...

    5. apache工作模式总结及网站访问缓慢处理记录

      apache目前主要有两种模式:prefork模式和worker模式:1)prefork模式(默认模式)prefork是Unix平台上的默认(缺省)MPM,使用多个子进程,每个子进程只有一个线程.每个 ...

    6. <a>标签中href="javascript:;"** 为什么 style不用src**

      &src/href <!--href 用于标示资源和文档关系,src 用于替换标签内容--> <img src="xxx.jpg"/> <sc ...

    7. 第三次作业 (一)----------------------Visual Studio 2015的安装及单元测试

      这是第三周的第一个作业,Visual Studio 2015的安装及单元测试. 我的电脑之前安装过Visual Studio 2015,但是在安装过程中我从来没有留意过各种注意事项,所集正好借此作业的 ...

    8. pair work结对编程(张艺 杨伊)

      一.结对编程人员: 张艺(学号后三位:185) 杨伊(学号后三位:151) 二.这是我们工作的样子:(图片) 三.结对编程优缺点:  优点:  1.结对编程时间紧密,在一定程度上可以督促双方学习,提高 ...

    9. 个人博客作业Week7(阅读文章,心得体会)

      Alpha阶段结束了,内心可以说是五味杂陈.不是说我们的产品拿不上台面那般差劲,复杂的心绪主要来源于和别的队的比较,别的队才刚刚发布没多久访问量和注册量就破百了,并且还发起了找bug送红包的活动.可能 ...

    10. 软件工程实践作业2 --梭哈游戏(java) 实践报告

      一,题目简介: 1.创建一副扑克牌 7------k 加入到集合对象中2.对扑克牌洗牌3.定义参与游戏的玩家的人,通过键盘输入,限定人数2-54.人数符合要求继续执行,不符合退出5.对玩家发牌,每个人 ...