听说这是动态开点主席树的板子题,但是发现我还不会,于是就来写一写

其实跟主席树一个样子的

这里就是存个板子吧

#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]旅行】的更多相关文章

  1. BZOJ 3531: [Sdoi2014]旅行 [树链剖分]

    3531: [Sdoi2014]旅行 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 1685  Solved: 751[Submit][Status] ...

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

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

  3. bzoj 3531 [Sdoi2014]旅行(树链剖分,线段树)

    3531: [Sdoi2014]旅行 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 876  Solved: 446[Submit][Status][ ...

  4. [SDOI2014]旅行

    洛谷 P3313 [SDOI2014]旅行 https://www.luogu.org/problem/show?pid=3313 题目描述 S国有N个城市,编号从1到N.城市间用N-1条双向道路连接 ...

  5. B20J_3231_[SDOI2014]旅行_树链剖分+线段树

    B20J_3231_[SDOI2014]旅行_树链剖分+线段树 题意: S国有N个城市,编号从1到N.城市间用N-1条双向道路连接,城市信仰不同的宗教,为了方便,我们用不同的正整数代表各种宗教. S国 ...

  6. [luogu P3313] [SDOI2014]旅行

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

  7. 洛谷 P3313 [SDOI2014]旅行 解题报告

    P3313 [SDOI2014]旅行 题目描述 S国有N个城市,编号从1到N.城市间用N-1条双向道路连接,满足从一个城市出发可以到达其它所有城市.每个城市信仰不同的宗教,如飞天面条神教.隐形独角兽教 ...

  8. bzoj 3531 [Sdoi2014]旅行 (树剖+线段树 动态开点)

    3531: [Sdoi2014]旅行 Time Limit: 40 Sec  Memory Limit: 512 MBSubmit: 2984  Solved: 1312[Submit][Status ...

  9. P3313 [SDOI2014]旅行

    P3313 [SDOI2014]旅行 树链剖分+动态线段树(并不是lct) 显然的,我们对于每一个宗教都要维护一个线段树. (那么空间不是爆炸了吗) 在这里引入:动态开点线段树 就是需要的点开起来,不 ...

  10. 3531: [Sdoi2014]旅行

    3531: [Sdoi2014]旅行 Time Limit: 20 Sec  Memory Limit: 512 MB Submit: 1731  Solved: 772 [Submit][Statu ...

随机推荐

  1. 2.3 if switch for等流程控制

    if条件中可以写多个语句,语句的作用域仅限于if,不可在if之外的地方使用 package main import ( "fmt" "io/ioutil" ) ...

  2. JQuery脚本-通过禁用按钮防止表单重复提交

    <script type="text/javascript"> /* jquer 脚本,避免重复提交 隐藏点击的submit,后在他之后插入同名button伪装成被隐藏 ...

  3. C# 批量生成随机密码,必须包含数字和字母,并用加密算法加密

    要求:密码必须包含数字和字母 思路:1.列出数字和字符. 组成字符串 :chars 2.利用randrom.Next(int i)返回一个小于所指定最大值的非负随机数. 3. 随机取不小于chars长 ...

  4. CF 303C——Minimum Modular——————【剪枝】

    Minimum Modular time limit per test 2 seconds memory limit per test 256 megabytes input standard inp ...

  5. MySql存储引擎MyISAM和InnoDB的区别

    1.MySQL默认采用的是MyISAM. 2.MyISAM不支持事务,而InnoDB支持.InnoDB的AUTOCOMMIT默认是打开的,即每条SQL语句会默认被封装成一个事务,自动提交,这样会影响速 ...

  6. String常用操作

    常量池: 字符串一旦被初始化就不会被改变 String s="123"; s="abc"; System.out.print(s); 这段代码看上去s的值是被改 ...

  7. ThreadPoolExecutor实现原理

    转载:https://blog.csdn.net/yanyan19880509/article/details/52718497 前言 做java开发的,一般都避免不了要面对java线程池技术,像to ...

  8. css中的单位和css中的颜色表示方法

    css中颜色的表示方式: 图片来源http://www.w3school.com.cn

  9. Django组件——forms组件

    一.校验字段功能 通过注册用户这个实例来学习校验字段功能. 1.模型:models.py from django.db import models # Create your models here. ...

  10. mysql 常用操作语句

    1 根据表中的其中一个字段的值来修改同行某字段的值 UPDATE  radar a INNER JOIN radar b ON a.id=b.id SET a.letter=LEFT(b.filena ...