链剖——BZOJ1036/Luogu2590 [ZJOI2008]树的统计
来源:互联网 发布:车载导航电子狗软件 编辑:程序博客网 时间:2024/05/18 03:30
题面:BZOJ1036 Luogu2590
树剖裸题,比luogu某模板题不知道水到哪里去了。。。
今天刚入门,看到一篇很好的树剖入门教程,写得很不错:传送门
对于初学者来说讲得很详细
最后来说说这道题
单点修改直接暴力找进线段树修改
求链上最大以及求和同链上修改,只不过不必更新权值
打得谨慎一点就行了
献上我的丑陋程序+注释
#include<cstdio>#include<cstring>#include<algorithm>#include<cmath>#include<iostream>#include<cstdlib>#include<string>#include<ctime>#include<queue>#include<climits>using namespace std;char c[10];int n,nedge=0,v[100001],p[200001],nex[200001],head[200001];int fa[100001],deep[100001],s[100001],son[100001],top[100001];int sx[100001],xs[100001],ne=0;int lt[400001],rt[400001],t[400001],mt[400001];inline void addedge(int a,int b){p[++nedge]=b;nex[nedge]=head[a];head[a]=nedge;}inline void dfs(int x,int Fa,int dep){//第一次dfs求父亲,子树数目,重儿子,深度 fa[x]=Fa;deep[x]=dep;s[x]=1; for(int k=head[x];k;k=nex[k]){ if(p[k]==Fa)continue; dfs(p[k],x,dep+1);s[x]+=s[p[k]]; if(!son[x]||s[p[k]]>s[son[x]])son[x]=p[k]; }}inline void dfss(int x,int bh){//第二次dfs求在线段树内位置,重链最顶端位置 top[x]=bh;sx[x]=++ne;xs[sx[x]]=x; if(!son[x])return; dfss(son[x],bh); for(int k=head[x];k;k=nex[k])if(p[k]!=son[x]&&p[k]!=fa[x])dfss(p[k],p[k]);}inline void build(int l,int r,int nod){//建线段树 lt[nod]=l;rt[nod]=r; if(l==r){t[nod]=mt[nod]=v[xs[l]];return;} int mid=l+r>>1; build(l,mid,nod*2);build(mid+1,r,nod*2+1); t[nod]=t[nod*2]+t[nod*2+1];mt[nod]=max(mt[nod*2],mt[nod*2+1]);}inline void xg(int x,int nod,int w){//线段树修改 if(lt[nod]==rt[nod]){t[nod]=mt[nod]=w;return;} int mid=lt[nod]+rt[nod]>>1; if(x<=mid)xg(x,nod*2,w); else xg(x,nod*2+1,w); t[nod]=t[nod*2]+t[nod*2+1];mt[nod]=max(mt[nod*2],mt[nod*2+1]); }inline int smax(int i,int j,int nod){//线段树找最大 if(lt[nod]>=i&&rt[nod]<=j)return mt[nod]; int mid=lt[nod]+rt[nod]>>1,ans=-1e9; if(i<=mid)ans=max(ans,smax(i,j,nod*2)); if(j>mid)ans=max(ans,smax(i,j,nod*2+1)); return ans;}inline int fmax(int x,int y){//树链上找最大 int fx=top[x],fy=top[y],ans=-1e9; while(fx!=fy){ if(deep[fx]<deep[fy])swap(fx,fy),swap(x,y); ans=max(ans,smax(sx[fx],sx[x],1)); x=fa[fx];fx=top[x]; } if(deep[x]>deep[y])swap(x,y); ans=max(ans,smax(sx[x],sx[y],1)); return ans;}inline int ssum(int i,int j,int nod){//线段树求和 if(lt[nod]>=i&&rt[nod]<=j)return t[nod]; int mid=lt[nod]+rt[nod]>>1,ans=0; if(i<=mid)ans+=ssum(i,j,nod*2); if(j>mid)ans+=ssum(i,j,nod*2+1); return ans;}inline int fsum(int x,int y){//树链上求和 int fx=top[x],fy=top[y],ans=0; while(fx!=fy){ if(deep[fx]<deep[fy])swap(fx,fy),swap(x,y); ans+=ssum(sx[fx],sx[x],1); x=fa[fx];fx=top[x]; } if(deep[x]>deep[y])swap(x,y); ans+=ssum(sx[x],sx[y],1); return ans;}int main(){ scanf("%d",&n); for(int i=1;i<n;i++){ int x,y;scanf("%d%d",&x,&y); addedge(x,y);addedge(y,x); } for(int i=1;i<=n;i++)scanf("%d",&v[i]); dfs(1,0,1);dfss(1,1);build(1,n,1); int m;scanf("%d",&m); for(int i=1;i<=m;i++){ int x,y;scanf("%s%d%d",c+1,&x,&y); if(c[1]=='C')xg(sx[x],1,y); if(c[2]=='M')printf("%d\n",fmax(x,y)); if(c[2]=='S')printf("%d\n",fsum(x,y)); } return 0;}
阅读全文
1 0
- 链剖——BZOJ1036/Luogu2590 [ZJOI2008]树的统计
- 【bzoj1036】 ZJOI2008—树的统计Count
- bzoj1036: [ZJOI2008]树的统计Count——树链剖分
- ZJOI2008树的统计bzoj1036
- BZOJ1036 ZJOI2008 树的统计
- 【BZOJ1036】【ZJOI2008】树的统计
- [Bzoj1036][ZJOI2008]树的统计Count
- BZOJ1036 [ZJOI2008]树的统计Count 树链剖分
- bzoj1036 codevs2460: [ZJOI2008]树的统计 树链剖分
- 【BZOJ1036】[ZJOI2008]树的统计Count 树链剖分
- 【BZOJ1036】【ZJOI2008】树的统计Count 树链剖分裸题
- 【ZJOI2008】【BZOJ1036】树的统计Count
- [BZOJ1036][ZJOI2008]树的统计Count && LCT
- bzoj1036: [ZJOI2008]树的统计Count
- [BZOJ1036][ZJOI2008]树的统计Count
- 【bzoj1036】【树链剖分】【ZJOI2008】树的统计Count
- 【BZOJ1036】【ZJOI2008】【树的统计count】【树链剖分】
- bzoj1036: [ZJOI2008]树的统计Count
- 收口思路
- apache错误:server's fully qualified domain name, using 127.0.0.1. Set the 'ServerName'
- Apache web服务
- 自定义进度条PictureProgressBar——从开发到开源发布全过程
- c++异常处理机制示例及讲解
- 链剖——BZOJ1036/Luogu2590 [ZJOI2008]树的统计
- 信息安全工程师 学习笔记(三十九)
- C语言调用C++函数
- 两条线段求交点编程
- Unity3D —— 手机设备震动接口
- iOS中tableview的执行顺序
- 各种聚类算法的比较
- 4 种不适合使用 NoSQL 数据库的场景
- sas infile和filename