【GDOI 2016 Day1】第四题 疯狂动物城

最新推荐文章于 2021-07-10 16:01:49 发布
最新推荐文章于 2021-07-10 16:01:49 发布
阅读量397 收藏
点赞数 2
分类专栏： 初级模板库及巧径
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/huangjingyuan107/article/details/77744762
版权

题目大意

给定一棵树，每个点x有一个权值a[x]，可以询问从x到y的路径上的S值。
具体的，

S=Σu∈x到y的路径Σu到y的距离i=1a[u]∗i
。
还有修改操作。每一次给点l,l+1,…,r的权值加上delta，这个算一个新的历史版本。
原来没修改的点的权值算历史版本0。
最后，还有一种操作，就是将权值回到历史版本x去。
操作中的x,y要疑惑一下上次输出的答案。

题解

题目就是要求
Σu∈x到lca(x,y)a[u]∗[1+...+(deep[u]+deep[y]−2∗deep[lca])]+Σu∈lca(x,y)到ya[u]∗[1+...+(deep[y]−deep[u]])
那么对于x到lca(x,y)，我们设对答案的贡献为S1。（除以2先不管）
设T=deep[y]-2*deep[lca]
S1=a[i]∗(deep[i]+T)∗(deep[i]+T+1)
那么对于lca(x,y)到y，我们设对答案的贡献为S2。（除以2先不管）
设T=deep[y]
S2=a[i]∗(−deep[i]+T)∗(−deep[i]+T+1)
那么我们只要在线段树上维护 a[i],a[i]∗deep[i],a[i]∗deep[i]2 就可以了。
为了方便修改，还要维护 deep[i],deep[i]2 。
然后进行树链剖分。
对了，还有回到第几个历史版本的操作，所以还要套主席树。

如何快捷地建立主席树？

首先，新的历史版本的节点维护的值是由旧历史版本的对应节点维护的值变化而来的。
我认为具体的就是先将旧历史版本的对应节点复制到新历史版本的节点去，然后根据题意修改新历史版本的节点维护的相应值。

树链剖分区间修改

步骤：
①比较deep[top[x]]和deep[top[y]]。
②如果deep[top[x]]>=deep[top[y]]就对路径x~top[x]进行操作，然后跳x，否则就对路径y~top[y]进行操作，然后跳y。转①
③如果top[x]=top[y]，终止循环，并对路径x~y进行操作。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 100010
#define mo 20160501
#define LL long long
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fd(i,a,b) for(i=a;i>=b;i--)
#define eg(i,x) for(i=head[x];i;i=edge[i].next)
using namespace std;
struct note{int to,next;};note edge[N*2];
struct sgm{
    int ls,rs;
    LL a,ad,ad2,d,d2,lazy;
};sgm tr[N*70];
LL i,j,k,l,n,m,T,temp1,temp2;
LL x,y,z,op,lastans;
int u,v,sgm_tot,ed;
int root[N];
LL head[N],val[N];
LL siz[N],dep[N],top[N],hv[N];
LL dfn[N],tar[N],fa[N][17];
LL tot,XXX,tt,ans;
void lb(int x,int y){edge[++tot].to=y;edge[tot].next=head[x];head[x]=tot;}
void dfs1(int x){
    siz[x]=1;
    for(int i=head[x];i;i=edge[i].next)
        if (!dep[edge[i].to]){
            dep[edge[i].to]=dep[x]+1; 
            fa[edge[i].to][0]=x;
            dfs1(edge[i].to);
            siz[x]+=siz[edge[i].to];
            if (!hv[x]||siz[edge[i].to]>siz[hv[x]])hv[x]=edge[i].to;
        }
}
void dfs2(int x){
    dfn[x]=++XXX;tar[XXX]=x;
    if (hv[x]) top[hv[x]]=top[x],dfs2(hv[x]);
    for(int i=head[x];i;i=edge[i].next){
        if (edge[i].to!=fa[x][0] && edge[i].to!=hv[x]){
            top[edge[i].to]=edge[i].to;
            dfs2(edge[i].to);
        }
    }
}
LL getlca(LL x,LL y){
    if (dep[x]<dep[y]) swap(x,y);
    int i;
    fd(i,16,0) if (dep[fa[x][i]]>=dep[y]) x=fa[x][i];
    if (x==y) return x;
    fd(i,16,0) if (fa[x][i]!=fa[y][i]) x=fa[x][i],y=fa[y][i];
    return fa[x][0];
}
void update(int ps){
    int L=tr[ps].ls,R=tr[ps].rs;
    tr[ps].a=(tr[L].a+tr[R].a)%mo;
    tr[ps].ad=(tr[L].ad+tr[R].ad)%mo;
    tr[ps].ad2=(tr[L].ad2+tr[R].ad2)%mo;
    tr[ps].d=(tr[L].d+tr[R].d)%mo;
    tr[ps].d2=(tr[L].d2+tr[R].d2)%mo;
}
void down(int ps,int l,int r){
    if (!tr[ps].lazy) return;
    if (l==r){
        tr[ps].lazy=0;
        return;
    }
    int wz=(l+r)/2,L=tr[ps].ls,R=tr[ps].rs,dt=tr[ps].lazy;
    tr[++sgm_tot]=tr[L];
    tr[sgm_tot].a=(tr[sgm_tot].a+dt*(wz-l+1))%mo;
    tr[sgm_tot].ad=(tr[sgm_tot].ad+dt*tr[sgm_tot].d)%mo;
    tr[sgm_tot].ad2=(tr[sgm_tot].ad2+dt*tr[sgm_tot].d2)%mo;
    tr[sgm_tot].lazy=(tr[sgm_tot].lazy+dt)%mo;
    tr[ps].ls=sgm_tot;
    tr[++sgm_tot]=tr[R];
    tr[sgm_tot].a=(tr[sgm_tot].a+dt*(r-wz))%mo;
    tr[sgm_tot].ad=(tr[sgm_tot].ad+dt*tr[sgm_tot].d)%mo;
    tr[sgm_tot].ad2=(tr[sgm_tot].ad2+dt*tr[sgm_tot].d2)%mo;
    tr[sgm_tot].lazy=(tr[sgm_tot].lazy+dt)%mo;
    tr[ps].rs=sgm_tot;
    tr[ps].lazy=0;
}
void build(int ps,int l,int r){
    if (l==r){
        tr[ps].a=val[tar[l]];
        tr[ps].ad=(tr[ps].a*dep[tar[l]])%mo;
        tr[ps].ad2=(tr[ps].ad*dep[tar[l]])%mo;
        tr[ps].d=dep[tar[l]];
        tr[ps].d2=dep[tar[l]]*dep[tar[l]];
        return;
    }
    int wz=(l+r)/2;
    tr[ps].ls=++sgm_tot;build(tr[ps].ls,l,wz);
    tr[ps].rs=++sgm_tot;build(tr[ps].rs,wz+1,r);
    update(ps);
}
void insert(int ps,int l,int r,int x,int y,int z){
    if (l==x && r==y){
        tr[ps].a=(tr[ps].a+(r-l+1)*z)%mo;
        tr[ps].ad=(tr[ps].ad+z*tr[ps].d)%mo;
        tr[ps].ad2=(tr[ps].ad2+z*tr[ps].d2)%mo;
        tr[ps].lazy=(tr[ps].lazy+z)%mo; 
        return;
    }
    int wz=(l+r)/2;
    int pp=0;
    if (tr[ps].lazy) pp=1; 
    down(ps,l,r);
    if (y<=wz){
        if (!pp) tr[++sgm_tot]=tr[tr[ps].ls],tr[ps].ls=sgm_tot;
        insert(tr[ps].ls,l,wz,x,y,z);
    } else
    if (x>wz){
        if (!pp) tr[++sgm_tot]=tr[tr[ps].rs],tr[ps].rs=sgm_tot;
        insert(tr[ps].rs,wz+1,r,x,y,z);
    } else{
        if (!pp) tr[++sgm_tot]=tr[tr[ps].ls],tr[ps].ls=sgm_tot;
        insert(tr[ps].ls,l,wz,x,wz,z);
        if (!pp) tr[++sgm_tot]=tr[tr[ps].rs],tr[ps].rs=sgm_tot;
        insert(tr[ps].rs,wz+1,r,wz+1,y,z);
    }
    update(ps);
}
void query(int ps,int l,int r,int x,int y,int z){
    if (l==x && r==y){
        ans=(ans+tr[ps].ad2+(tr[ps].a*(T*T+T)%mo)%mo)%mo;
        ans=(ans+z*tr[ps].ad*(2*T+1)%mo+mo)%mo;
        return;
    }
    down(ps,l,r);
    int wz=(l+r)/2;
    if (y<=wz)query(tr[ps].ls,l,wz,x,y,z);else
        if (x>wz)query(tr[ps].rs,wz+1,r,x,y,z);else{
            query(tr[ps].ls,l,wz,x,wz,z);
            query(tr[ps].rs,wz+1,r,wz+1,y,z);
        }
}
void oper(int x,int y,LL z){
    LL lca=getlca(x,y);
    temp1=dep[y]-2*dep[lca];
    temp2=dep[y];
    bool pp=0;
    LL tq=1; 
    while (1){
        LL tx=top[x],ty=top[y];
        if (dep[tx]<dep[ty]){
            swap(x,y);
            swap(tx,ty);
            tq=-tq;
        }
        if (tx==ty) break;
        if (tq==-1) T=temp2;else T=temp1;
        if(z==-1)query(root[ed],1,n,dfn[tx],dfn[x],tq);
            else insert(root[ed],1,n,dfn[tx],dfn[x],z);
        x=fa[tx][0];
    }
    if (dfn[x]>dfn[y]) swap(x,y),tq=-tq;
    if (tq==-1) T=temp1;else T=temp2;
    if (z==-1) query(root[ed],1,n,dfn[x],dfn[y],-tq);
          else insert(root[ed],1,n,dfn[x],dfn[y],z);
}
int main(){
    freopen("zootopia.in","r",stdin);
    freopen("zootopia.out","w",stdout);
    scanf("%lld%lld",&n,&m);
    fo(i,1,n-1){
        scanf("%lld%lld",&u,&v);
        lb(u,v);lb(v,u);
    }
    dep[1]=1;dfs1(1);
    top[1]=1;dfs2(1);
    fo(j,1,16)fo(i,1,n)fa[i][j]=fa[fa[i][j-1]][j-1];
    lastans=0;
    fo(i,1,n) scanf("%lld",&val[i]);
    root[0]=sgm_tot=1;
    build(root[0],1,n);
    ed=0;
    fo(i,1,m){
        scanf("%lld",&op);
        if (op==1){
            scanf("%lld%lld%lld",&x,&y,&z);
            x^=lastans;y^=lastans;
            root[++tt]=++sgm_tot;
            tr[sgm_tot]=tr[root[ed]];
            ed=tt;
            oper(x,y,z);
        } else
        if (op==2){
            scanf("%lld%lld",&x,&y);
            x^=lastans;y^=lastans;
            ans=0;
            oper(x,y,-1);
            ans=(ans*10080251)%mo;
            printf("%lld\n",ans);
            lastans=ans;
        } else{
            scanf("%lld",&x);
            x^=lastans;
            ed=x;
        }
    }
    return 0;
}
Luckfort
关注
专栏目录
SSL集训 2021.07.16 提高B组 GDOI2017 day2 小学生语文题【DP】【细节思维题】【难】
JA_yichao欢迎你 cnblogs账号：https://www.cnblogs.com/mayichao/
07-17 153
SSL集训
【GDOI2016模拟4.5】刺客 题解
CMC_YXY的博客
03-22 285
【GDOI2016模拟4.5】刺客 Description SillyHook是一个著名的刺客，虽然他比较Silly，但他精通打洞，善于深入敌方内部刺杀敌人。 现在SillyHook接到了若干个消灭敌人的任务，每次任务中他都会装备一把耐久度为m 的钩刃，并打洞潜入敌方内部。有n 个待消灭的敌人SillyHook消灭第i 个敌人需要消耗钩刃Ai 点耐久度（如果有其他武器则可以选择用其他武器消灭这个敌人，且不消耗钩刃耐久度），之后便得到他的武器，可以用来消灭任意Bi 个敌人。但SillyHook太Silly了
[GDOI2016][树链剖分+主席树]疯狂动物城
EZ_TJY的博客
07-19 1613
题面： Description Nick 是只在动物城以坑蒙拐骗为生的狐狸，儿时受到偏见的伤害，放弃了自己的理想。他被兔子 Judy 设下圈套，被迫与她合作查案，而卷入意想不到的阴谋，历尽艰险后成为搭档。他们识破了绵羊副市长 Bellwether 的计划，发现是 Bellwether 陷害食肉动物，用毒药让食肉动物发狂。Bellwether 被抓到了监狱里面， Nick 和 Judy 过上了一...
【GDOI2016】疯狂动物城（树链剖分+可持久化线段树）
一个蒟蒻的博客
08-14 542
码农题…… 调了我三个晚上…… 看来我的代码能力还是太弱了…… 首先我们不难发现在u到v这条链的答案为∑ni=1(n−i)(n−i+1)ai2∑i=1n(n−i)(n−i+1)ai2\sum_{i=1}^n\frac{(n-i)(n-i+1)a_i}{2} 然后把它拆开可以得到答案为∑ni=1(n+1)nai−(2n+1)iai+i2ai∑i=1n(n+1)nai−(2n+1)iai+i2...
【GDOI2016】疯狂动物城 题解
Four's
08-18 1659
【题目大意】      ~~~~~~n个节点的一棵树，有三种操作。       ~~~~~~1：将x到y的路径上的所有点的点权加上delta       ~~~~~~2：询问x到y的答案。答案的计算为：对于路径上的点i，设它到y的距离为s，则i的贡献为1加到s。       ~~~~~~3：将这棵树恢复到第x次1操作之后的版本。       ~~~~~~操作数为m，强制在线。
【GDOI 2016】第四题 疯狂动物城
Facico的博客
08-06 1672
DescriptionSolution题目好长，一看就烦操作1：1,x,y,z把x到y的点权值加z。 操作2：2,x,y把x到y的点一次放进一个数组a里面，统计答案ans=∑ni=1a[i]∗∑n−ij=1ans=\sum_{i=1}^na[i]*\sum_{j=1}^{n-i}。 操作3：3,x所有点的状态回到第x次修改后的。码农题啊！！！害人啊！！！不过知道是码农题，在比赛的最后30分钟，竟
【xsy2440】【GDOI2016】疯狂动物城
weixin_30565199的博客
06-25 123
感受一下这恐怖的题目长度~~~ 其实题意很裸，但是作为GDOI的一道防AK题，自然有他优秀的地方。 简化题意：给出一棵树，要求支持三个操作： 1.修改点值 2.询问点$x$到$y$之间的一些东东 3.回到某个版本之前 可持久化，强制在线。 考虑拆一下询问 如果把$x$到$y$路径上的所有点值放到一个数组$a$里，那么询问式子就是： $\sum\limits_{i=1}...
GDOI2016第一次终极模拟总结
PhilipsWeng的专栏
04-23 1031
感觉发挥的还算平稳，但每天都有分数没有搞到，不能做到分分必拿，还是比较亏的。Day1一开始看题，第一题一看题就开始推公式，发现要维护个凸壳，还要维护前后凸壳，看起来非常地奇怪，再看第二题，看数据范围都知道要矩阵乘法，推了一下状态发现比较简单，然后看第三题，画风比较奇怪，应该要求一个欧几里德最小生成树，然后模拟一下，但不会求，最后一题不怎么懂，只会没有障碍的20分做法。看完题想了一会儿大概过了1个钟，
【GDOI2018 Day1】密码锁 题解
CMC_YXY的博客
07-10 278
【GDOI2018 Day1】密码锁 Description Input Output Sample Input Sample Input1 4 3 1 2 1 0 Sample Input2 11 8 1 2 3 4 5 0 5 4 3 2 1 Sample Input3 20 100 30 91 15 72 61 41 10 37 98 41 94 80 26 96 10 88 59 5 84 14 Sample Output Sample Output1 2 Sample Output2 8 S
【GDOI 2016 Day1】疯狂动物城
chen1352
06-01 1001
【GDOI 2016 Day1】疯狂动物城
GDOI 2016 Day1 T4 疯狂动物城
时光真疯狂, 我一路执迷于匆忙.
05-17 1045
Description给出一个N个节点的数，和M次操作。每次操作的类型如下： 1,x,y,z，将x到y的路径上的ai加上z 2,x,y，询问x到y的路径上，ai*(1+2+..+n-i)的和 3,x，将所有的a变更回第x次修改之后的状态。 强制在线。 N,M<=10^5.Solution码农题（hehe） 愉快的农民生活 先考虑一下序列上的问题。 Ans=∑i=1nai∗(n−i+1
【GDOI 2016 Day1】第四题 疯狂动物城 题解+代码
Jacky35
05-23 914
我恨码农题！！！题目大意一棵树，三个操作 1、从X到Y的路径上每个点的权值加Z 2、查询X到Y的一些神奇东西 3、将树回复到第X次操作后的状态输入第一行n,m。n个点，m个操作 接下来n-1行每行两个数表示X连到Y 接下来一行n个数：每个点的初始权值 接下来M行对应上面的三种操作（其中X和Y要亦或上次的答案）题解deep[i]表示深度，a[i]为i的权值 其中有40% 的数据是一条链，
csust2019集训队选拔赛题解
JCC
03-06 1754
第一次参加比赛 崩的实在是太惨了.. 只A了一道板子题(已经被自己菜哭了 没想到大家都掉线了 一起崩 这也给了我一个血的教训 以后比赛一定要多开几题 简直是在死撑着坐满五个小时 当时其实老早就想溜了(以后比赛再也不吃香蕉了.. 肚子快疼炸了 接下来我就按我自己的顺序来写题解(我怕我先写难的写不下去就咕了 //耗时一星期终于竣工了emmmm 哭泣的阿木木 众所周知，阿木木是英雄联盟里面最...
模拟题【字符串、枚举、贪心】字符串
ZYzyZzzz的博客
09-27 212
字符串 【题目描述】 小林与亮亮正在做一个游戏。小林随意地写出一个字符串，字符串只由大写 字母组成，然后指定一个非负整数mmm，亮亮可以进行至多mmm次操作，每次操作 为交换相邻两个字符。亮亮的目标是使得操作后的字符串出现最长相同的字符的 长度最大。你能帮亮亮计算一下这个最大长度是多少吗？ 【输入格式】 第一行一个字符串SSS。 第二行一个整数mmm。 【输出格式】 只有一个整数，表示所求的最大长度。 【样例输入】 ABCCDCDDC 4 【样例输出】 4 【数据规模】 对于20%20\%20%的数据，
GDOI2016总结
WorldWide_D的博客
05-03 1076
Day0第三次参加GDOI，在初中已经是全勤了。前两次一次二等，一次一等。这一次考前挺有信心的。我复习了一些模板，大致总结了前面模拟犯过的错误。 比赛来到了四会，一个比较远的地方，住在城郊的一个宾馆，没有什么地方可以走走。。。Day1在第一场我就感觉到题目相比前几次加大了难度，结果（0+20+0+20）T1：中学生数学题。出题人本来想弄个很水的二次函数题的，结果一个向下取整让我在草稿纸上推了很久，
#191-[LCA]祖孙询问
2017gdgzoi999的博客
12-15 885
Description 已知一棵n个节点的有根树。有m个询问。每个询问给出了一对节点的编号x和y，询问x与y的祖孙关系。 Input 输入第一行包括一个整数n表示节点个数。 接下来n行每行一对整数对a和b表示a和b之间有连边。如果b是-1，那么a就是树的根。 第n+2行是一个整数m表示询问个数。 接下来m行，每行两个正整数x和y。 Output 对于每一个询问，输出1:如果x是y的...
GDKOI2021：信息学冬令营普及组真题解析
"2021年广东省重点中学信息学冬令营(GDKOI2021)的普及组第一试试题，包含了四道题目，分别是地图、灌水、配对和旅行。每道题都有特定的输入输出文件名、时间空间限制以及满分分数。竞赛环境采用NOI系列活动的标准...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值