HDU 3966 Aragorn's Story (树链剖分)
发布时间：2020-12-13 20:43:21 所属栏目：PHP教程 来源：网络整理
导读：题目地址：HDU 3966 树链剖分模板题。不多说、、 代码以下： #include iostream #include string.h #include math.h #include queue #include algorithm #include stdlib.h #include map #include set #include stdio.h using namespace std ; #define LL lo
|
题目地址：HDU 3966 #include <iostream>
#include <string.h>
#include <math.h>
#include <queue>
#include <algorithm>
#include <stdlib.h>
#include <map>
#include <set>
#include <stdio.h>
using namespace std;
#define LL long long
#define pi acos(⑴.0)
#pragma comment(linker,"/STACK:1024000000,1024000000")
const int mod=1e9+7;
const int INF=0x3f3f3f3f;
const double eqs=1e⑶;
const int MAXN=50000+10;
#define root 1,n,1
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
int head[MAXN],cnt,a[MAXN],n;
int dep[MAXN],siz[MAXN],son[MAXN],w[MAXN],top[MAXN],tot,fa[MAXN];
int sum[MAXN<<2],lazy[MAXN<<2];
struct node
{
int u,v,next;
}edge[MAXN<<1];
void add(int u,int v)
{
edge[cnt].v=v;
edge[cnt].u=u;
edge[cnt].next=head[u];
head[u]=cnt++;
}
void init()
{
memset(head,-1,sizeof(head));
cnt=0;
memset(dep,0,sizeof(dep));
memset(son,sizeof(son));
tot=0;
memset(sum,sizeof(sum));
memset(lazy,sizeof(lazy));
}
void dfs1(int u,int p)
{
siz[u]=1;
for(int i=head[u];i!=-1;i=edge[i].next){
int v=edge[i].v;
if(v==p) continue ;
fa[v]=u;
dep[v]=dep[u]+1;
dfs1(v,u);
siz[u]+=siz[v];
if(siz[son[u]]<siz[v]) son[u]=v;
}
}
void dfs2(int u,int tp)
{
top[u]=tp;w[u]=++tot;
if(son[u]) dfs2(son[u],tp);
for(int i=head[u];i!=-1;i=edge[i].next){
int v=edge[i].v;
if(v!=fa[u]&&v!=son[u])
dfs2(v,v);
}
}
struct Line_Tree
{
void PushUp(int rt)
{
sum[rt]=sum[rt<<1]+sum[rt<<1|1];
}
void PushDown(int l,int r,int rt)
{
int mid=l+r>>1;
if(lazy[rt]){
sum[rt<<1]+=(mid-l+1)*lazy[rt];
sum[rt<<1|1]+=(r-mid)*lazy[rt];
lazy[rt<<1]+=lazy[rt];
lazy[rt<<1|1]+=lazy[rt];
lazy[rt]=0;
}
}
void Update(int ll,int rr,int x,int l,int rt)
{
if(ll<=l&&rr>=r){
lazy[rt]+=x;
sum[rt]+=x*(r-l+1);
return ;
}
int mid=l+r>>1;
PushDown(l,rt);
if(ll<=mid) Update(ll,rr,x,lson);
if(rr>mid) Update(ll,rson);
PushUp(rt);
}
int Query(int p,int rt)
{
if(l==r){
return sum[rt];
}
int mid=l+r>>1;
PushDown(l,rt);
if(p<=mid) return Query(p,lson);
return Query(p,rson);
}
}lt;
void Update(int u,int v,int x)
{
int f1=top[u],f2=top[v];
while(f1!=f2){
if(dep[f1]<dep[f2]){
swap(u,v);
swap(f1,f2);
}
lt.Update(w[f1],w[u],root);
u=fa[f1];f1=top[u];
}
if(dep[u]<dep[v]){
swap(u,v);
}
lt.Update(w[v],root);
}
int main()
{
int m,q,i,j,u,x;
char s[10];
while(scanf("%d%d%d",&n,&m,&q)!=EOF){
for(i=1;i<=n;i++){
scanf("%d",&a[i]);
}
init();
while(m--){
scanf("%d%d",&u,&v);
add(u,v);
add(v,u);
}
dfs1(1,-1);
dfs2(1,1);
for(i=1;i<=n;i++){
lt.Update(w[i],w[i],a[i],root);
}
while(q--){
scanf("%s",s);
if(s[0]=='I'){
scanf("%d%d%d",&v,&x);
Update(u,x);
}
else if(s[0]=='D'){
scanf("%d%d%d",-x);
}
else{
scanf("%d",&x);
printf("%d
",lt.Query(w[x],root));
}
}
}
return 0;
}
（编辑：李大同） 【声明】本站内容均来自网络，其相关言论仅代表作者个人观点，不代表本站立场。若无意侵犯到您的权利，请及时与联系站长删除相关内容! |