题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3966

题意:

给出一棵树,并给定各个点权的值,然后有3种操作:

I C1 C2 K: 把C1与C2的路径上的所有点权值加上K

D C1 C2 K:把C1与C2的路径上的所有点权值减去K

Q C:查询节点编号为C的权值

分析:

典型的树链剖分,对节点进行操作,可以用树状数组或者线段树。

树链剖分+树状数组:

 #include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring> using namespace std;
const int maxn = ; struct Edge
{
int to,next;
}edge[maxn*];
int head[maxn],cnt,num;
int a[maxn],n,m,p,c[maxn];
int size[maxn],top[maxn],id[maxn],fa[maxn],son[maxn],dep[maxn]; void init()
{
memset(head,-,sizeof(head));
memset(son,-,sizeof(son));
memset(c,,sizeof(c));
cnt=;
num=;
}
void addedge(int u,int v)
{
edge[cnt].to = v;
edge[cnt].next = head[u];
head[u] = cnt++;
} void dfs_1(int u,int f,int d)
{
dep[u]=d;
size[u]=;
fa[u]=f;
for(int i=head[u];i!=-;i=edge[i].next)
{
int v=edge[i].to;
if(v==f)
continue;
dfs_1(v,u,d+);
size[u]+=size[v];
if(son[u]==-||size[son[u]]<size[v])
son[u]=v;
}
} void dfs_2(int u,int tp)
{
top[u] = tp;
id[u] = ++num;
if(son[u]!=-)
dfs_2(son[u],tp);
for(int i=head[u];i!=-;i=edge[i].next)
{
int v=edge[i].to;
if(v==fa[u]||v==son[u])
continue;
dfs_2(v,v);
}
} int lowbit(int x)
{
return x&-x;
} int sum(int x)
{
int res=;
while(x>)
{
res+=c[x];
x-=lowbit(x);
}
return res;
} void add(int x,int d)
{
while(x<=n)
{
c[x]+=d;
x+=lowbit(x);
}
} void change(int u,int v,int val)
{
int tp1=top[u],tp2=top[v];
while(tp1!=tp2)
{
if(dep[tp1]<dep[tp2])
{
swap(tp1,tp2);
swap(u,v);
}
add(id[tp1],val);
add(id[u]+,-val);
u=fa[tp1];
tp1=top[u];
}
if(dep[u]>dep[v])
swap(u,v);
add(id[u],val);
add(id[v]+,-val);
} int main()
{
while(~scanf("%d%d%d",&n,&m,&p))
{
init();
for(int i=;i<=n;i++)
scanf("%d",&a[i]);
for(int i=;i<=m;i++)
{
int u,v;
scanf("%d%d",&u,&v);
addedge(u,v);
addedge(v,u);
}
dfs_1(,,);
dfs_2(,);
for(int i=;i<=n;i++)
{
add(id[i],a[i]);
add(id[i]+,-a[i]);
}
char s[];
int c1,c2,k,c;
for(int i=;i<p;i++)
{
scanf("%s",s);
if(s[]=='I')
{
scanf("%d%d%d",&c1,&c2,&k);
change(c1,c2,k);
}
if(s[]=='D')
{
scanf("%d%d%d",&c1,&c2,&k);
change(c1,c2,-k);
}
if(s[]=='Q')
{
scanf("%d",&c);
cout<<sum(id[c])<<endl;
}
}
}
return ;
}

树链剖分+线段树:

 #include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
const int maxn = ;
int n,m,p;
int val[maxn],a[maxn];
struct Edge
{
int to,next;
}edge[maxn*];
int head[maxn],add[maxn*];
int dep[maxn],fa[maxn],size[maxn],son[maxn],top[maxn],id[maxn];
int cnt,num;
void init()
{
memset(head,-,sizeof(head));
memset(son,-,sizeof(son));
memset(add,,sizeof(add));
cnt=;
num=;
} void addedge(int u,int v)
{
edge[cnt].to=v;
edge[cnt].next=head[u];
head[u]=cnt++;
} void dfs_1(int u,int f,int d)
{
dep[u]=d;
size[u]=;
fa[u]=f;
for(int i=head[u];i!=-;i=edge[i].next)
{
int v=edge[i].to;
if(v==f)
continue;
dfs_1(v,u,d+);
size[u]+=size[v];
if(son[u]==-||size[son[u]]<size[v])
son[u]=v;
}
} void dfs_2(int u,int tp)
{
top[u] = tp;
id[u] = ++num;
if(son[u]!=-)
dfs_2(son[u],tp);
for(int i=head[u];i!=-;i=edge[i].next)
{
int v = edge[i].to;
if(v == fa[u] ||v == son[u])
continue;
dfs_2(v,v);
}
} struct Tree
{
int left,right;
int sum;
}tree[maxn*]; void pushup(int i)
{
tree[i].sum = tree[i*].sum + tree[i*+].sum;
} void build(int i,int begin,int end)
{
tree[i].left=begin;
tree[i].right=end;
if(begin==end)
{
tree[i].sum=val[begin];
return;
}
int mid=(begin+end)/;
build(i*,begin,mid);
build(i*+,mid+,end);
pushup(i);
} void pushdown(int i)
{
if(add[i])
{
add[i*] += add[i];
add[i*+] += add[i];
int mid=(tree[i].left+tree[i].right)/;
tree[i*].sum += add[i]*(mid-tree[i].left+);
tree[i*+].sum += add[i]*(tree[i].right-mid);
add[i]=;
}
} void update(int i,int begin,int end,int value)
{
if(tree[i].left>=begin&&tree[i].right<=end)
{
add[i]+=value;
tree[i].sum+=value*(tree[i].right-tree[i].left+);
return;
}
pushdown(i);
int mid=(tree[i].left+tree[i].right)/;
if(mid>=begin)
update(i*,begin,end,value);
if(mid<end)
update(i*+,begin,end,value);
pushup(i);
} void change(int u,int v,int value)
{
int tp1=top[u],tp2=top[v];
while(tp1!=tp2)
{
if(dep[tp1]<dep[tp2])
{
swap(tp1,tp2);
swap(u,v);
}
update(,id[tp1],id[u],value);
u = fa[tp1];
tp1 = top[u];
}
if(dep[u]>dep[v])
swap(u,v);
update(,id[u],id[v],value);
} long long query(int i,int begin,int end)
{
if(tree[i].left>=begin&&tree[i].right<=end)
return tree[i].sum;
pushdown(i);
int mid=(tree[i].left+tree[i].right)/;
long long ans=;
if(mid>=begin)
ans+=query(i*,begin,end);
if(mid<end)
ans+=query(i*+,begin,end);
return ans;
} int main()
{
while(~scanf("%d%d%d",&n,&m,&p))
{
init();
for(int i=;i<=n;i++)
scanf("%d",&a[i]);
for(int i=;i<=m;i++)
{
int u,v;
scanf("%d%d",&u,&v);
addedge(u,v);
addedge(v,u);
}
dfs_1(,,);
dfs_2(,); for(int i=;i<=n;i++)
val[id[i]]=a[i];
build(,,n);
char s[];
int c1,c2,k,c;
for(int i=;i<p;i++)
{
scanf("%s",s);
if(s[]=='I')
{
scanf("%d%d%d",&c1,&c2,&k);
change(c1,c2,k);
}
if(s[]=='D')
{
scanf("%d%d%d",&c1,&c2,&k);
change(c1,c2,-k);
}
if(s[]=='Q')
{
scanf("%d",&c);
cout<<query(,id[c],id[c])<<endl;
}
}
}
return ;
}

hdu 3966 Aragorn's Story(树链剖分+树状数组/线段树)的更多相关文章

  1. Qtree3题解(树链剖分(伪)+线段树+set)

    外话:最近洛谷加了好多好题啊...原题入口 这题好像是SPOJ的题,挺不错的.看没有题解还是来一篇... 题意: 很明显吧.. 题解: 我的做法十分的暴力:树链剖分(伪)+线段树+\(set\)... ...

  2. Luogu 2590 [ZJOI2008]树的统计 / HYSBZ 1036 [ZJOI2008]树的统计Count (树链剖分,LCA,线段树)

    Luogu 2590 [ZJOI2008]树的统计 / HYSBZ 1036 [ZJOI2008]树的统计Count (树链剖分,LCA,线段树) Description 一棵树上有n个节点,编号分别 ...

  3. 【bzoj4999】This Problem Is Too Simple! 树链剖分+动态开点线段树

    题目描述 给您一颗树,每个节点有个初始值. 现在支持以下两种操作: 1. C i x(0<=x<2^31) 表示将i节点的值改为x. 2. Q i j x(0<=x<2^31) ...

  4. [bzoj 3531][SDOI2014]旅行(树链剖分+动态开点线段树)

    题目:http://www.lydsy.com:808/JudgeOnline/problem.php?id=3531 分析: 对于每个颜色(颜色<=10^5)都建立一颗线段树 什么!那么不是M ...

  5. 洛谷P3313 [SDOI2014]旅行(树链剖分 动态开节点线段树)

    题意 题目链接 Sol 树链剖分板子 + 动态开节点线段树板子 #include<bits/stdc++.h> #define Pair pair<int, int> #def ...

  6. 刷题总结——骑士的旅行(bzoj4336 树链剖分套权值线段树)

    题目: Description 在一片古老的土地上,有一个繁荣的文明. 这片大地几乎被森林覆盖,有N座城坐落其中.巧合的是,这N座城由恰好N-1条双 向道路连接起来,使得任意两座城都是连通的.也就是说 ...

  7. BZOJ 3531 [Sdoi2014]旅行 树链剖分+动态开点线段树

    题意 S国有N个城市,编号从1到N.城市间用N-1条双向道路连接,满足从一个城市出发可以到达其它所有城市.每个城市信仰不同的宗教,如飞天面条神教.隐形独角兽教.绝地教都是常见的信仰. 为了方便,我们用 ...

  8. hdu3966 树链剖分点权模板+线段树区间更新/树状数组区间更新单点查询

    点权树的模板题,另外发现树状数组也是可以区间更新的.. 注意在对链进行操作时方向不要搞错 线段树版本 #include<bits/stdc++.h> using namespace std ...

  9. 【BZOJ3531】[Sdoi2014]旅行 树链剖分+动态开点线段树

    [BZOJ3531][Sdoi2014]旅行 Description S国有N个城市,编号从1到N.城市间用N-1条双向道路连接,满足从一个城市出发可以到达其它所有城市.每个城市信仰不同的宗教,如飞天 ...

  10. bzoj3531——树链剖分+动态开点线段树

    3531: [Sdoi2014]旅行 Time Limit: 20 Sec  Memory Limit: 512 MB Description S国有N个城市,编号从1到N.城市间用N-1条双向道路连 ...

随机推荐

  1. Pattern Recognition And Machine Learning (模式识别与机器学习) 笔记 (1)

    By Yunduan Cui 这是我自己的PRML学习笔记,目前持续更新中. 第二章 Probability Distributions 概率分布 本章介绍了书中要用到的概率分布模型,是之后章节的基础 ...

  2. 电赛总结(二)&mdash;&mdash;AD芯片总结之音频处理芯片ADC8009

    一.特性参数 1.专门用来用音频处理的AD芯片 2.内部固定好8K的采样速率 3.8位AD芯片 二.内部结构图 三.芯片管脚图 四.管脚功能说明 管脚名称 功能 IN0~IN7 数据输入端 ABC 数 ...

  3. Android drawBitmapMesh扭曲图像

    今天介绍一下在Android当中怎么扭曲图像,在Android系统中的Canvas提供了一个drawBitmapMesh(Bitmap bitmap, int meshWidth, int meshH ...

  4. 说说C#和.NET的关系

    .NET和C#的关系:C#语言可以通过.NET平台编写.部署.运行.NET应用程序..NET可以支持包括C#在内的多种语言编写的应用程序..NET平台主要包括FCL(框架类库)它是微软事先定义好的类的 ...

  5. android自动填充短信验证码

    自动拦截短信实际上就是在系统注册一个BroadcastReceiver,然后通过设置拦截短信的: filter.addAction("android.provider.Telephony.S ...

  6. swfupload操作手册

    SWFUpload SWFUpload 最初是Vinterwebb.se 开发的客户端文件上传工具.它联合javascript和flash,在浏览器中提供一个优于传统上传标签 <input ty ...

  7. 解决安装Visual Studio 2012后SQL Server 2008 远程过程调用失败的问题

    安装了Visual Studio 2012后,打开SQL Server 2008配置管理器,发现了一个问题.如下图 解决办法:

  8. jbpm4.3表结构和表字段说明

    整体说明: JBPM4_DEPLOYMENT   流程定义表JBPM4_DEPLOYPROP 流程定义属性表 JBPM4_EXECUTION  流程实例表JBPM4_HIST_ACTINST 流程活动 ...

  9. 搭建真正的zookeeper集群

    搭建zookeeper伪分布式集群 zookeeper是Hadop Ecosystem中非常重要的组件,它的主要功能是为分布式系统提供一致性协调服务, 提供的功能包括配置维护,域名服务,分布式同步和组 ...

  10. MFC单文档视图程序简介

    在视图应用程序中,应用程序的数据由文档对象代表,数据的视图由视图对象代表.MFC的Cdocument类是文档对象的基类,Cview类是视图对象的基类.应用程序的主窗口,其操作功能在MFC的Cframe ...