首页 > 代码库 > [BZOJ 3531][Sdoi2014]旅行（树链剖分+线段树）
[BZOJ 3531][Sdoi2014]旅行（树链剖分+线段树）
Description
S国有N个城市，编号从1到N。城市间用N-1条双向道路连接，满足从一个城市出发可以到达其它所有城市。每个城市信仰不同的宗教，如飞天面条神教、隐形独角兽教、绝地教都是常见的信仰。为了方便，我们用不同的正整数代表各种宗教， S国的居民常常旅行。旅行时他们总会走最短路，并且为了避免麻烦，只在信仰和他们相同的城市留宿。当然旅程的终点也是信仰与他相同的城市。S国政府为每个城市标定了不同的旅行评级，旅行者们常会记下途中（包括起点和终点）留宿过的城市的评级总和或最大值。
在S国的历史上常会发生以下几种事件：
”CC x c”：城市x的居民全体改信了c教；
”CW x w”：城市x的评级调整为w;
”QS x y”：一位旅行者从城市x出发，到城市y，并记下了途中留宿过的城市的评级总和；
”QM x y”：一位旅行者从城市x出发，到城市y，并记下了途中留宿过
的城市的评级最大值。
由于年代久远，旅行者记下的数字已经遗失了，但记录开始之前每座城市的信仰与评级，还有事件记录本身是完好的。请根据这些信息，还原旅行者记下的数字。 为了方便，我们认为事件之间的间隔足够长，以致在任意一次旅行中，所有城市的评级和信仰保持不变。
Solution
只需要对每个宗教建一棵线段树就好啦
#include<iostream>#include<cstdlib>#include<cstdio>#include<cstring>#define MAXN 100005using namespace std;int n,q,w[MAXN],c[MAXN],head[MAXN],cnt=0;int father[MAXN],top[MAXN],maxv[MAXN],siz[MAXN],deep[MAXN],pos[MAXN],sz=0;int tot=0,rt[MAXN],ls[MAXN*100],rs[MAXN*100],maxn[MAXN*100],sum[MAXN*100];int read(){ int x=0,f=1;char c=getchar(); while(c<‘0‘||c>‘9‘){if(c==‘-‘)f=-1;c=getchar();} while(c>=‘0‘&&c<=‘9‘){x=x*10+c-‘0‘;c=getchar();} return x*f;}struct Node1{ int next,to;}Edges[MAXN*2];void addedge(int u,int v){ Edges[++cnt].next=head[u]; head[u]=cnt; Edges[cnt].to=v;}void dfs1(int u){ siz[u]=1; for(int i=head[u];~i;i=Edges[i].next) { int v=Edges[i].to; if(v==father[u])continue; father[v]=u,deep[v]=deep[u]+1; dfs1(v); siz[u]+=siz[v]; if(siz[v]>siz[maxv[u]])maxv[u]=v; }}void dfs2(int u,int t){ top[u]=t; ++sz,pos[u]=sz; if(maxv[u])dfs2(maxv[u],t); for(int i=head[u];~i;i=Edges[i].next) { int v=Edges[i].to; if(v==maxv[u]||v==father[u])continue; dfs2(v,v); }}void update(int idx){ sum[idx]=sum[ls[idx]]+sum[rs[idx]]; maxn[idx]=max(maxn[ls[idx]],maxn[rs[idx]]);}void change(int &idx,int l,int r,int p,int w){ if(!idx)++tot,idx=tot; if(l==r){maxn[idx]=sum[idx]=w;return;} int mid=(l+r)>>1; if(p<=mid)change(ls[idx],l,mid,p,w); else change(rs[idx],mid+1,r,p,w); update(idx);}int querysum(int idx,int l,int r,int L,int R){ if(!idx)return 0; if(L==l&&R==r)return sum[idx]; int mid=(l+r)>>1; if(R<=mid)return querysum(ls[idx],l,mid,L,R); else if(L>mid)return querysum(rs[idx],mid+1,r,L,R); else return querysum(ls[idx],l,mid,L,mid)+querysum(rs[idx],mid+1,r,mid+1,R);}int querymax(int idx,int l,int r,int L,int R){ if(!idx)return 0; if(L==l&&R==r)return maxn[idx]; int mid=(l+r)>>1; if(R<=mid)return querymax(ls[idx],l,mid,L,R); else if(L>mid)return querymax(rs[idx],mid+1,r,L,R); else return max(querymax(ls[idx],l,mid,L,mid),querymax(rs[idx],mid+1,r,mid+1,R));}int QS(int x,int y){ int t=c[x],res=0; while(top[x]!=top[y]) { if(deep[top[x]]<deep[top[y]])swap(x,y); res+=querysum(rt[t],1,sz,pos[top[x]],pos[x]); x=father[top[x]]; } if(deep[x]<deep[y])swap(x,y); res+=querysum(rt[t],1,sz,pos[y],pos[x]); return res;}int QM(int x,int y){ int t=c[x],res=0; while(top[x]!=top[y]) { if(deep[top[x]]<deep[top[y]])swap(x,y); res=max(res,querymax(rt[t],1,sz,pos[top[x]],pos[x])); x=father[top[x]]; } if(deep[x]<deep[y])swap(x,y); res=max(res,querymax(rt[t],1,sz,pos[y],pos[x])); return res;}int main(){ memset(head,-1,sizeof(head)); n=read(),q=read(); for(int i=1;i<=n;i++)w[i]=read(),c[i]=read(); for(int i=1;i<n;i++) { int x=read(),y=read(); addedge(x,y),addedge(y,x); } deep[1]=1;dfs1(1),dfs2(1,1); for(int i=1;i<=n;i++)change(rt[c[i]],1,sz,pos[i],w[i]); for(int i=1;i<=q;i++) { char opt[5]; scanf("%s",opt); int x=read(),y=read(); if(opt[1]==‘C‘) { change(rt[c[x]],1,sz,pos[x],0); change(rt[y],1,sz,pos[x],w[x]); c[x]=y; } else if(opt[1]==‘W‘) { change(rt[c[x]],1,sz,pos[x],y); w[x]=y; } else if(opt[1]==‘S‘) printf("%d\n",QS(x,y)); else printf("%d\n",QM(x,y)); } return 0;}
[BZOJ 3531][Sdoi2014]旅行（树链剖分+线段树）
声明：以上内容来自用户投稿及互联网公开渠道收集整理发布，本网站不拥有所有权，未作人工编辑处理，也不承担相关法律责任，若内容有误或涉及侵权可进行投诉： 投诉/举报 工作人员会在5个工作日内联系你，一经查实，本站将立刻删除涉嫌侵权内容。