兰州大学第一届 飞马杯 ★★快乐苹果树★★ 树链剖分 + 懒标记 + 树状数组

传送门

文章目录

  • 题意:
  • 思路:

题意:

在这里插入图片描述

思路:

第一次听说树链剖分能在fa[top[i]]fa[top[i]]fa[top[i]]的地方加懒标记,学到了学到了。

首先不能被题目吓住,这个题目仔细剖析一下不难发现一些性质:
以下用seise_isei表示iii子树的大小。
(1)(1)(1)uuu选择了fff点,那么vvv选任一点都可以,也就是让fff为根的子树的每个点权值都加1sef1sefk\frac{1}{se_f}\frac{1}{se_f}ksef1sef1k
(2)(2)(2)uuu选择了非fff点的时候,其选择的点假设落在与fff直接相连的xxx点所在的子树中,那么vvv只能选非xxx的子树中的点,也就是让非xxx的子树的点都加上sexsef1sef−sexk\frac{se_x}{se_f}\frac{1}{se_f-se_x}ksefsexsefsex1k
以上两个操作显然可以先将fff所在的子树都加上1sef1sefk+∑sexsef1sef−sexk\frac{1}{se_f}\frac{1}{se_f}k+\sum \frac{se_x}{se_f}\frac{1}{se_f-se_x}ksef1sef1k+sefsexsefsex1k,现在就有个比较显然的做法就是遍历fff的出边给子树减去sexsef1sef−sexk\frac{se_x}{se_f}\frac{1}{se_f-se_x}ksefsexsefsex1kfff的全部子树都加上1sef1sefk+∑sexsef1sef−sexk\frac{1}{se_f}\frac{1}{se_f}k+\sum \frac{se_x}{se_f}\frac{1}{se_f-se_x}ksef1sef1k+sefsexsefsex1k即可。
但这样的话给你个菊花图就炸了,复杂度升到qnlognqnlognqnlogn,由于是对子树操作的,比较容易想到能不能加个懒标记呢?我们加了懒标记又如何在计算答案的时候能减去这个懒标记呢?解决了这个两个问题我们就可以顺利解决这个问题了。
树上倍增解决这个题不是很容易,我们考虑重链剖分。
首先我们加点的时候,思路就是按照上面哪个思路,即先给所有点都加上,再给子树减去某个值,所以我们lazylazylazy标记应该是标记了减去了多少。首先加操作显然可以写一个dfsdfsdfs序,让后用树状数组维护一下。考虑书剖的特殊性质,我们只需要在tag[f]tag[f]tag[f]加上kkk,让后将重儿子直接减去sexsef1sef−sexk\frac{se_x}{se_f}\frac{1}{se_f-se_x}ksefsexsefsex1k,因为重儿子在之后跳toptoptop的时候会直接跳过的,而轻儿子在跳到这条链的top[i]top[i]top[i]的时候,需要减去s2[top[i]]∗tag[fa[top[i]]]s2[top[i]]*tag[fa[top[i]]]s2[top[i]]tag[fa[top[i]]]s2s2s2是预处理的sexsef1sef−sex\frac{se_x}{se_f}\frac{1}{se_f-se_x}sefsexsefsex1
这样问题就顺利解决啦~

// Problem: ★★快乐苹果树★★
// Contest: NowCoder
// URL: https://ac.nowcoder.com/acm/contest/16520/K
// Memory Limit: 524288 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>
#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("---")
#define lowbit(x) ((x)&(-x))
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=200010,mod=998244353,INF=0x3f3f3f3f;
const double eps=1e-6;int n,m;
vector<int>v[N];
int se[N],son[N],top[N],fa[N],in[N],tot;
LL inv[N],s1[N],s2[N],s3[N];
LL tr[N],tag[N];void dfs1(int u,int f) {se[u]=1; fa[u]=f;for(auto x:v[u]) {if(x==f) continue;dfs1(x,u);se[u]+=se[x];if(se[son[u]]<se[x]) son[u]=x;}s1[u]=inv[se[u]]*inv[se[u]]%mod;for(auto x:v[u]) {if(x==f) continue;s2[x]=se[x]*inv[se[u]]%mod*inv[se[u]-se[x]]%mod;s3[u]+=s2[x]; s3[u]%=mod;}
}void dfs2(int u,int t) {top[u]=t; in[u]=++tot;if(son[u]) dfs2(son[u],t); for(auto x:v[u]) {if(x==fa[u]||x==son[u]) continue;dfs2(x,x);}
}void add(int x,LL c) {for(int i=x;i<=n;i+=lowbit(i)) tr[i]+=c,tr[i]%=mod,tr[i]+=mod,tr[i]%=mod;
}LL query(int x) {LL ans=0;for(int i=x;i;i-=lowbit(i)) ans+=tr[i],ans%=mod;return ans;
}int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);inv[1]=1;for(int i=2;i<N;i++)inv[i]=(mod-mod/i)*inv[mod%i]%mod;int _; scanf("%d",&_);while(_--) {scanf("%d%d",&n,&m); tot=0;for(int i=1;i<=n;i++) se[i]=son[i]=top[i]=fa[i]=s3[i]=s2[i]=s1[i]=0,v[i].clear();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);while(m--) {int op,x; scanf("%d%d",&op,&x);if(op==1) {int k; scanf("%d",&k);add(in[x],1ll*(s1[x]+s3[x])%mod*k%mod); add(in[x]+se[x],(((mod-s1[x]-s3[x])%mod+mod)%mod)*k%mod);if(son[x]) add(in[son[x]],(mod-s2[son[x]]*k%mod)%mod),add(in[son[x]]+se[son[x]],s2[son[x]]*k%mod);tag[x]+=k; tag[x]%=mod;} else {LL ans=query(in[x]);for(int i=top[x];fa[i];i=top[fa[i]]) {ans-=tag[fa[i]]*s2[i]%mod,ans+=mod,ans%=mod;}printf("%lld\n",ans%mod);}}for(int i=1;i<=n;i++) tr[i]=tag[i]=0;}return 0;
}
/**/

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

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

相关文章

【BZOJ3328】PYXFIB【矩阵快速幂】【单位根反演】【二项式定理】

传送门 题意&#xff1a; ∑i0⌊nk⌋(nik)Fik\sum_{i0}^{\lfloor\frac nk\rfloor}\binom n{ik}F_{ik}i0∑⌊kn​⌋​(ikn​)Fik​ FFF为斐波拉契数列 n≤1e18,k≤2e4,p≤1e9n\leq 1e18,k\leq2e4,p\leq1e9n≤1e18,k≤2e4,p≤1e9且为质数且模kkk余111 显然就是求 ∑i0n[k∣i](ni…

.NET分布式框架 | Orleans 知多少

引言公司物联网项目集成Orleans以支持高并发的分布式业务&#xff0c;对于Orleans也是第一次接触&#xff0c;本文就分享下个人对Orleans的理解。这里先抛出自己的观点&#xff1a;Orleans 是一个支持有状态云生应用/服务水平伸缩的基于Virtual Actor 模型的.NET分布式框架。下…

nowcoder 牛牛的最大兴趣组 质因子 + 思维

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 首先nnn很小的话可以暴力连边&#xff0c;让后染个色求一个颜色最多的即可。但是这个题显然不行&#xff0c;由于是三次方&#xff0c;所以考虑质因子入手。 首先很容易就能想到将所有的数…

ASP.NET Core on K8S学习初探(3)部署API到K8S

“ 终于可以部署ASP.NET Core到K8S中了...”在上一篇《基本概念快速一览》中&#xff0c;我们把基本的一些概念快速地简单地不求甚解地过了一下&#xff0c;本篇开始我们会将ASP.NET Core WebAPI部署到K8S&#xff0c;从而结束初探的旅程。01—准备一个WebAPI这里准备一个空的A…

K-D Tree学习笔记

引入 K-D Tree 是一种处理高维空间的数据结构。 支持O(nk−1k)O(n^{\frac {k-1}k})O(nkk−1​)查询给定超矩形内的点的信息&#xff0c; kkk 为维数。 可以用替罪羊树的思想实现动态插入。 但其实用的最多的是在错误的复杂度内查询奇奇怪怪的点信息。 构建 K-D Tree 是平…

Codeforces Round #633 B. Edge Weight Assignment 结论题 + dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 首先考虑最小值&#xff0c;如果从一个叶子结点出发到任意叶子的距离都为偶数&#xff0c;那么只需要一个值就可以满足条件。如果有奇数的&#xff0c;考虑111 ^ 222 ^ 303030&#xff0c;…

ASP.NET Core Web API中使用Swagger

本节导航Swagger介绍在ASP.NET CORE 中的使用swagger在软件开发中,管理和测试API是一件重要而富有挑战性的工作。在我之前的文章《研发团队,请管好你的API文档》 也专门阐述了通过文档管理工具,来保证API文档和代码的一致性,这样更加有助于团队的协作。以往我们总是通过第三方平…

【ZJOI2015】幻想乡战略游戏【点分树】【带权重心】

题意&#xff1a;nnn个点带边权的树&#xff0c;动态修改点权viv_ivi​&#xff0c;最小化 钦定一个点xxx 后 ∑idist(x,i)∗vi\sum\limits_{i} dist(x,i)*v_ii∑​dist(x,i)∗vi​的值。 n,q≤105n,q \leq10^5n,q≤105&#xff0c;度数不超过202020 限制度数的树上的一些诡异…

P4559 [JSOI2018]列队 主席树

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你nnn个学生以及其位置&#xff0c;mmm个询问&#xff0c;每次询问[l,r][l,r][l,r]的学生跑到[k,kr−l][k,kr-l][k,kr−l]的位置的最小总的移动距离。当然不同的学生不能站在同一位置。 n,m≤5e5,1≤ai,k≤…

C# 8.x 先睹为快

尽管 C# 8.0 还需要几个月的时间才能发布&#xff0c;但是 C# 8.x 的特性计划已经开始了。该计划中的特性有一些是新提出的&#xff0c;而另外一些则是之前在 C# 8 中考虑过的。和往常一样&#xff0c;这个特性列表仍可能会发生变化。如之前的一个报告所讨论的那样&#xff0c;…

【WC2014】紫荆花之恋【替罪羊思想】【动态点分树】【替罪羊树】

传送门 题意&#xff1a;有一棵开始时没有结点的树&#xff0c;nnn次询问&#xff0c;每次新加一点并给定父结点、到父亲的距离、参数rir_iri​,并询问满足dist(u,v)≤rurvdist(u,v)\leq r_ur_vdist(u,v)≤ru​rv​的点对(u,v)(u,v)(u,v)的对数。 n≤105n\leq 10^5n≤105&…

NuGet是什么?理解与使用

如果你了解python&#xff0c;那么它类似pip。如果你了解nodejs&#xff0c;那么它类似npm。如果你了解ruby&#xff0c;那么它类似gem。对&#xff0c;它就是一个包&#xff08;package&#xff09;管理平台&#xff0c;确切的说是 .net平台的包管理工具&#xff0c;它提供了一…

P2468 [SDOI2010]粟粟的书架 主席树 + 二分 + 二维前缀和

传送门 文章目录题意&#xff1a;题意&#xff1a;题意&#xff1a; 给你一个n∗mn*mn∗m的矩阵&#xff0c;每次询问一个矩形&#xff0c;左上角是(x1,y1)(x_1,y_1)(x1​,y1​)&#xff0c;右上角是(x2,y2)(x_2,y_2)(x2​,y2​)&#xff0c;和一个hhh&#xff0c;问是否能选出…

带花树学习笔记

引入 带花树解决的是一般图的最大匹配问题。 学习此算法前你需要掌握二分图匹配的匈牙利算法&#xff0c;至少需要知道它的思想。 理论 众所周知&#xff0c;二分图匹配的思想是不断地找增广路。 严谨地讲&#xff0c;是找到了一条简单路径&#xff0c;它的起点和终点都未…

你竟然没用 Nuget 构建项目?

想写这篇文章缘起于公众号文章里的一条留言&#xff1a;对于一个现代化的开发平台&#xff0c;建立一种让开发者创建&#xff0c;分享与使用可复用代码的机制是十分必要的。这种“可复用代码”被打包后的文件通常被称作“包”&#xff08;package&#xff09;&#xff0c;对于.…

Codeforces Round #339 (Div. 1) D. Kingdom and its Cities 虚树 + dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; n≤1e5,∑k≤1e5,q≤1e5n\le1e5,\sum k\le1e5,q\le1e5n≤1e5,∑k≤1e5,q≤1e5。 思路&#xff1a; 经过分析&#xff0c;外敌占领的城市只有可能是两点的lcalcalca或者两点之间任意一点&#xff0c;是哪个点…

小白开学Asp.Net Core《四》 —— 使用AspectCore-Framework

小白开学Asp.Net Core《四》—— 使用AspectCore-Framework一、AspectCore-Framework说AspectCore-Framework不得不先谈谈的AOP&#xff0c;AOP&#xff1a;在软件业&#xff0c;AOP为Aspect Oriented Programming的缩写&#xff0c;意为&#xff1a;面向切面编程&#xff0c;通…

【WC2016】挑战NPC 【带花树】【建图】

传送门 题意&#xff1a;有nnn个球和mmm个筐&#xff0c;每个筐最多放333个球&#xff0c;每个球只能放入特定的一些筐中&#xff0c;在题中给出。构造一种放球的方案使得nnn个球都被放在某个筐中且 球的个数不超过111 的筐的数量尽量大。 m≤100,n≤3mm\leq 100,n\leq 3mm≤1…

P4103 [HEOI2014]大工程 虚树 + dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 简化一下题意就是求树上给定点集中每两个点之间的距离之和&#xff0c;相距最远的点&#xff0c;相距最近的点。 对于距离和我们统计一下边的贡献就好啦&#xff0c;边两头的sizesizesize乘…

【翻译】Keras.NET简介 - 高级神经网络API in C#

用C#&#xff0c;远离996Keras.NETKeras.NET是一个高级神经网络API&#xff0c;它使用C#编写&#xff0c;并带有Python绑定&#xff0c;可以在Tensorflow、CNTK或Theano上运行。其关注点是实现快速实验。因为做好研究的关键是&#xff1a;能在尽可能短的时间内从一个想法发展出…