【[SDOI2014]旅行】
听说这是动态开点主席树的板子题,但是发现我还不会,于是就来写一写
其实跟主席树一个样子的
这里就是存个板子吧
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define re register
#define max(a,b) ((a)>(b)?(a):(b))
#define maxn 100005
struct E
{
int v,nxt;
}e[maxn<<1];
inline int read()
{
char c=getchar();
int x=0;
while(c<'0'||c>'9') c=getchar();
while(c>='0'&&c<='9')
x=(x<<3)+(x<<1)+c-48,c=getchar();
return x;
}
int head[maxn],deep[maxn],top[maxn],son[maxn],sz[maxn],fa[maxn],to[maxn],b[maxn],c[maxn],a[maxn],_[maxn];
int rt[maxn];
int l[maxn*20],r[maxn*20],sum[maxn*20],d[maxn*20];
int cnt,num,tot,n,Q;
void dfs1(int x)
{
sz[x]=1;
int maxx=-1;
for(re int i=head[x];i;i=e[i].nxt)
if(!deep[e[i].v])
{
deep[e[i].v]=deep[x]+1;
fa[e[i].v]=x;
dfs1(e[i].v);
sz[x]+=sz[e[i].v];
if(sz[e[i].v]>maxx) maxx=sz[e[i].v],son[x]=e[i].v;
}
}
void dfs2(int x,int topf)
{
to[x]=++tot;
b[tot]=a[x];
c[tot]=_[x];
top[x]=topf;
if(!son[x]) return;
dfs2(son[x],topf);
for(re int i=head[x];i;i=e[i].nxt)
if(deep[e[i].v]>deep[x]&&son[x]!=e[i].v) dfs2(e[i].v,e[i].v);
}
inline void add_edge(int x,int y)
{
e[++num].v=y;
e[num].nxt=head[x];
head[x]=num;
}
int change(int now,int x,int y,int pos,int val)
{
int root=now;
if(!now) root=++cnt;
if(x==y)
{
d[root]=val,sum[root]=val;
return root;
}
int mid=x+y>>1;
if(pos<=mid) l[root]=change(l[root],x,mid,pos,val);
else r[root]=change(r[root],mid+1,y,pos,val);
d[root]=max(d[l[root]],d[r[root]]);
sum[root]=sum[l[root]]+sum[r[root]];
return root;
}
int query(int now,int x,int y,int ll,int rr)
{
if(!now) return 0;
if(ll<=x&&rr>=y) return sum[now];
int mid=x+y>>1;
if(rr<=mid) return query(l[now],x,mid,ll,rr);
if(ll>mid) return query(r[now],mid+1,y,ll,rr);
return query(l[now],x,mid,ll,rr)+query(r[now],mid+1,y,ll,rr);
}
int ask(int now,int x,int y,int ll,int rr)
{
if(!now) return 0;
if(ll<=x&&rr>=y) return d[now];
int mid=x+y>>1;
if(rr<=mid) return ask(l[now],x,mid,ll,rr);
if(ll>mid) return ask(r[now],mid+1,y,ll,rr);
return max(ask(l[now],x,mid,ll,rr),ask(r[now],mid+1,y,ll,rr));
}
inline int ask_sum(int x,int y,int col)
{
int s=0;
while(top[x]!=top[y])
{
if(deep[top[x]]<deep[top[y]]) std::swap(x,y);
s+=query(rt[col],1,n,to[top[x]],to[x]);
x=fa[top[x]];
}
if(deep[x]>deep[y]) std::swap(x,y);
s+=query(rt[col],1,n,to[x],to[y]);
return s;
}
inline int ask_max(int x,int y,int col)
{
int s=0;
while(top[x]!=top[y])
{
if(deep[top[x]]<deep[top[y]]) std::swap(x,y);
s=max(s,ask(rt[col],1,n,to[top[x]],to[x]));
x=fa[top[x]];
}
if(deep[x]>deep[y]) std::swap(x,y);
s=max(s,ask(rt[col],1,n,to[x],to[y]));
return s;
}
int main()
{
n=read(),Q=read();
for(re int i=1;i<=n;i++)
{
a[i]=read(),_[i]=read();
if(!rt[_[i]]) rt[_[i]]=++cnt;
}
int x,y;
for(re int i=1;i<n;i++)
x=read(),y=read(),add_edge(x,y),add_edge(y,x);
deep[1]=1,dfs1(1),dfs2(1,1);
for(re int i=1;i<=n;i++)
change(rt[c[i]],1,n,i,b[i]);
char opt[2];
int val;
while(Q--)
{
scanf("%s",opt);
if(opt[0]=='Q'&&opt[1]=='S')
{
x=read(),y=read();
printf("%d\n",ask_sum(x,y,_[x]));
}
if(opt[0]=='Q'&&opt[1]=='M')
{
x=read(),y=read();
printf("%d\n",ask_max(x,y,_[x]));
}
if(opt[0]=='C'&&opt[1]=='W')
{
x=read(),val=read();
b[to[x]]=val;
change(rt[_[x]],1,n,to[x],val);
}
if(opt[0]=='C'&&opt[1]=='C')
{
x=read(),val=read();
change(rt[val],1,n,to[x],b[to[x]]);
change(rt[_[x]],1,n,to[x],0);
_[x]=val;
}
}
return 0;
}
【[SDOI2014]旅行】的更多相关文章
- BZOJ 3531: [Sdoi2014]旅行 [树链剖分]
3531: [Sdoi2014]旅行 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 1685 Solved: 751[Submit][Status] ...
- 【BZOJ3531】[Sdoi2014]旅行 树链剖分+动态开点线段树
[BZOJ3531][Sdoi2014]旅行 Description S国有N个城市,编号从1到N.城市间用N-1条双向道路连接,满足从一个城市出发可以到达其它所有城市.每个城市信仰不同的宗教,如飞天 ...
- bzoj 3531 [Sdoi2014]旅行(树链剖分,线段树)
3531: [Sdoi2014]旅行 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 876 Solved: 446[Submit][Status][ ...
- [SDOI2014]旅行
洛谷 P3313 [SDOI2014]旅行 https://www.luogu.org/problem/show?pid=3313 题目描述 S国有N个城市,编号从1到N.城市间用N-1条双向道路连接 ...
- B20J_3231_[SDOI2014]旅行_树链剖分+线段树
B20J_3231_[SDOI2014]旅行_树链剖分+线段树 题意: S国有N个城市,编号从1到N.城市间用N-1条双向道路连接,城市信仰不同的宗教,为了方便,我们用不同的正整数代表各种宗教. S国 ...
- [luogu P3313] [SDOI2014]旅行
[luogu P3313] [SDOI2014]旅行 题目描述 S国有N个城市,编号从1到N.城市间用N-1条双向道路连接,满足从一个城市出发可以到达其它所有城市.每个城市信仰不同的宗教,如飞天面条神 ...
- 洛谷 P3313 [SDOI2014]旅行 解题报告
P3313 [SDOI2014]旅行 题目描述 S国有N个城市,编号从1到N.城市间用N-1条双向道路连接,满足从一个城市出发可以到达其它所有城市.每个城市信仰不同的宗教,如飞天面条神教.隐形独角兽教 ...
- bzoj 3531 [Sdoi2014]旅行 (树剖+线段树 动态开点)
3531: [Sdoi2014]旅行 Time Limit: 40 Sec Memory Limit: 512 MBSubmit: 2984 Solved: 1312[Submit][Status ...
- P3313 [SDOI2014]旅行
P3313 [SDOI2014]旅行 树链剖分+动态线段树(并不是lct) 显然的,我们对于每一个宗教都要维护一个线段树. (那么空间不是爆炸了吗) 在这里引入:动态开点线段树 就是需要的点开起来,不 ...
- 3531: [Sdoi2014]旅行
3531: [Sdoi2014]旅行 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 1731 Solved: 772 [Submit][Statu ...
随机推荐
- API 接口设计工具 --Swagger
swagger-editor,无法启动GUI软件,在线版的FQ也打不开 null
- shell 函数与内置变量
1,特殊shell变量 $# 传递到脚本的参数个数 $* 以一个单字符串显示所有向脚本传递的参数 $$ 脚本运行的当前进程ID号 $! 后台运行的最后一个进程的ID号 $@ 与$*相同,但是使用时加引 ...
- (转)netstat 命令详解
netstat 命令详解 原文:https://www.cnblogs.com/xieshengsen/p/6618993.html netstat命令是一个监控TCP/IP网络的非常有用的工具,它 ...
- 推荐文章unity框架与工具
https://www.indienova.com/u/kuaile/blogread/1343
- Excel&&word&&PPT
1. Excel 1.1 制作下拉框 选中单元格或列--> 菜单"数据" --> "数据验证"-->"设置" --> ...
- matlab安装过程的被要求的配置程序
顺序是这样的: 网址的顺序是这样的: 1. http://cn.mathworks.com/support/compilers/R2015b/index.html?sec=win64&s_ci ...
- JavaScript字符串去除空格
/*字符串去除空格*/ String.prototype.Trim = function() { return this.replace(/(^\s*)|(\s*$)/g, "") ...
- c# winform 欢迎界面时加载数据
using System; using System.Collections.Generic; using System.ComponentModel; using System.Data; usin ...
- java常用API之System类
System中代表程序所在系统,提供了对应的一些系统属性信息,和系统操作.System类不能手动创建对象,因为构造方法被private修饰,阻止外界创建对象.System类中的都是static方法,类 ...
- 【数据库】4.0 MySQL入门学习(四)——linux系统环境下MySQL安装
1.0 我的操作系统是CentOS Linux release 7.6.1810 (Core) 系统详细信息如下: Linux version 3.10.0-957.1.3.el7.x86_64 ( ...