HDU3966 树链剖分

题目：HDU3966

题意：一棵树上每个节点有一些值，三个操作（其实是俩）
操作I是增加一条链所有点的值
操作D是减少一条链所有点的值
操作Q是查询某个点的值
网上说得手工扩栈，我没扩栈也A掉了

//#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <cstdio>
#include <iostream>
#include <cstring>
#include <vector>
using namespace std;
#define REP(I,N) for (I=0;I<N;I++)
#define rREP(I,N) for (I=N-1;I>=0;I--)
#define rep(I,S,N) for (I=S;I<N;I++)
#define rrep(I,S,N) for (I=N-1;I>=S;I--)
#define FOR(I,S,N) for (I=S;I<=N;I++)
typedef unsigned long long ULL;
typedef long long LL;
const int INF=0x3f3f3f3f;
const LL INFF=0x3f3f3f3f3f3f3f3fll;
const LL M=1e9+7;
const LL maxn=50007;
const double eps=0.00000001;
LL gcd(LL a,LL b){return b?gcd(b,a%b):a;}
template<typename T>inline T abs(T a,T b) {return a>0?a:-a;}int lazy[maxn*4],a[maxn];//用lazy就好吧。。这个不会pushup
void pushdown(int x){if (lazy[x]){lazy[x<<1]+=lazy[x];lazy[x<<1|1]+=lazy[x];lazy[x]=0;}
}
void build(int x,int l,int r){lazy[x]=0;if (l==r) {lazy[x]=a[l];return;}int mid=(l+r)/2;build(x<<1,l,mid);build(x<<1|1,mid+1,r);
}
void update(int x,int l,int r,int val,int L,int R){if (l<=L&&R<=r){lazy[x]+=val;return;}pushdown(x);int mid=(L+R)/2;if (mid>=l) update(x<<1,l,r,val,L,mid);if (r>mid) update(x<<1|1,l,r,val,mid+1,R);
}
int query(int x,int pos,int L,int R){if (L==R) return lazy[x];pushdown(x);int mid=(L+R)/2;if (mid>=pos) return query(x<<1,pos,L,mid);return query(x<<1|1,pos,mid+1,R);
}int n,q,s,T;
int i,j,k;
int u,v,len;
vector<int> edge[maxn];
int fa[maxn],son[maxn],sz[maxn],top[maxn],id[maxn],dep[maxn];
int b[maxn];
int tot;
void dfs1(int u,int from,int depth){int v,i,mx=-1;sz[u]=1;fa[u]=from;dep[u]=depth;son[u]=0;REP(i,edge[u].size()){v=edge[u][i];if (v==from) continue;dfs1(v,u,depth+1);sz[u]+=sz[v];if (sz[v]>mx) mx=sz[v],son[u]=v;}
}
void dfs2(int u,int x){int v,i;top[u]=x;id[u]=++tot;if (son[u]) dfs2(son[u],x);REP(i,edge[u].size()){v=edge[u][i];if (v==fa[u]||v==son[u]) continue;dfs2(v,v);}
}
void Change(int x,int y,int val){while (top[x]!=top[y]){if (dep[top[x]]<dep[top[y]]) swap(x,y);update(1,id[top[x]],id[x],val,1,tot);x=fa[top[x]];}if (dep[x]>dep[y]) swap(x,y);update(1,id[x],id[y],val,1,tot);
}int main(){while (~scanf("%d%d%d",&n,&i,&q)){FOR(i,1,n) scanf("%d",&b[i]);FOR(i,1,n) edge[i].clear();REP(i,n-1){scanf("%d%d",&u,&v);edge[u].push_back(v);edge[v].push_back(u);}tot=0;dfs1(1,0,1);dfs2(1,1);FOR(i,1,n) a[id[i]]=b[i];//映射一下build(1,1,tot);REP(i,q){char c;int c1,c2,k;c=' ';while (c!='I'&&c!='D'&&c!='Q') c=getchar();if (c=='I'){scanf("%d%d%d",&c1,&c2,&k);Change(c1,c2,k);}if (c=='D'){scanf("%d%d%d",&c1,&c2,&k);Change(c1,c2,-k);}if (c=='Q'){scanf("%d",&c1);printf("%d\n",query(1,id[c1],1,tot));}}}
}
/*
*/

HDU3966 树链剖分相关推荐

HDU3966(树链剖分)
题目:Aragorn's Story 题意:给一棵树,并给定各个点权的值,然后有3种操作: I C1 C2 K: 把C1与C2的路径上的所有点权值加上K D C1 C2 K:把C1与C2的路径上的所有 ...
hdu3966 树链剖分点权模板+线段树区间更新/树状数组区间更新单点查询
点权树的模板题,另外发现树状数组也是可以区间更新的.. 注意在对链进行操作时方向不要搞错线段树版本 #include<bits/stdc++.h> using namespace std ...
[hdu3966 Aragorn's Story]树链剖分
题意:要求在一棵N(<=50000)个带权节点的树上支持3种操作 (1)I u v w,u到v的路径上每个节点权值增加w (2)D u v w,u到v的路径上每个节点权值减少w (3)Q u,求 ...
SPOJ375(树链剖分)
题目:Query on a tree 题意:给定一棵树,告诉了每条边的权值,然后给出两种操作: (1)把第i条边的权值改为val (2)询问a,b路径上权值最大的边分析:本题与HDU3966差不多, ...
SPOJ- QTREE+HDU 3966（树链剖分裸题
题目:维护一个树,支持修改边长,查询任意两点间距离. 思路:学习了一下树链剖分,还是看了一阵子才看懂的(大概两个多小时orz...)其实总的来说并不是很难,主要是数组比较多,然后看的那篇博客大概是为了 ...
【ZJOI2008】树的统计（树链剖分）
传送门 Solution: 就是树链剖分入门题啦~ // luogu-judger-enable-o2 #include<bits/stdc++.h> #define N 30005 #d ...
树链剖分+线段树 HDOJ 4897 Little Devil I（小恶魔）
题目链接题意: 给定一棵树,每条边有黑白两种颜色,初始都是白色,现在有三种操作: 1 u v:u到v路径(最短)上的边都取成相反的颜色 2 u v:u到v路径上相邻的边都取成相反的颜色(相邻即仅有一 ...
P1505 [国家集训队]旅游树链剖分
题目链接题意:树上更新某一点权值,更新两点简单路径权值,查询最大,最小,和思路:思路应该比较简单,就是树链剖分后用线段树维护区间最大最小以及区间和. 但是本题比较特殊的是给的边权,转化为点权即可. ...
【模板】树链剖分 P3384
题目链接 //部分转自:https://www.luogu.org/problemnew/solution/P3384 初学树链剖分,感觉这个模板题还是容易理解的,但是实在是码量很大的. 知识点: 重 ...
树链剖分——线段树区间合并bzoj染色
线段树区间合并就挺麻烦了,再套个树链就更加鬼畜,不过除了代码量大就没什么其他的了.. 一些细节:线段树每个结点用结构体保存,pushup等合并函数改成返回一个结构体,这样好写一些 struct Seg ...

HDU3966 树链剖分

题目：HDU3966

HDU3966 树链剖分相关推荐

最新文章

热门文章