首页 > 代码库 > bzoj3531 [Sdoi2014]旅行

bzoj3531 [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,并记下了途中留宿过的城市的评级最大值。
由于年代久远,旅行者记下的数字已经遗失了,但记录开始之前每座城市的信仰与评级,还有事件记录本身是完好的。请根据这些信息,还原旅行者记下的数字。

为了方便,我们认为事件之间的间隔足够长,以致在任意一次旅行中,所有城市的评级和信仰保持不变。

Input

    输入的第一行包含整数N,Q依次表示城市数和事件数。
    接下来N行,第i+l行两个整数Wi,Ci依次表示记录开始之前,城市i的评级和信仰。
    接下来N-1行每行两个整数x,y表示一条双向道路。
    接下来Q行,每行一个操作,格式如上所述。

Output

    对每个QS和QM事件,输出一行,表示旅行者记下的数字。

Sample Input

5 6
3 1
2 3
1 2
3 3
5 1
1 2
1 3
3 4
3 5
QS 1 5
CC 3 1
QS 1 5
CW 3 3
QS 1 5
QM 2 4

Sample Output

8
9
11
3

HINT

N,Q < =10^5    , C < =10^5

 

 数据保证对所有QS和QM事件,起点和终点城市的信仰相同;在任意时刻,城市的评级总是不大于10^4的正整数,且宗教值不大于C。

 

正解:树链剖分+线段树。

这题很水啊,直接树链剖分+线段树就过了。但是这是作业,所以还是来水一下。。不过这题要建n棵线段树,所以必须动态开结点才行。。修改和查询直接照着题目的意思去弄就行了。。

 

  1 //It is made by wfj_2048~  2 #include <algorithm>  3 #include <iostream>  4 #include <complex>  5 #include <cstring>  6 #include <cstdlib>  7 #include <cstdio>  8 #include <vector>  9 #include <cmath> 10 #include <queue> 11 #include <stack> 12 #include <map> 13 #include <set> 14 #define inf (1<<30) 15 #define N (100010) 16 #define il inline 17 #define RG register 18 #define ll long long 19 #define File(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout) 20  21 using namespace std; 22  23 struct edge{ int nt,to; }g[2*N]; 24  25 int head[N],top[N],fa[N],son[N],tid[N],dep[N],sz[N],w[N],c[N],n,Q,num,cnt,ccnt; 26 int sum[300*N],mx[300*N],ls[300*N],rs[300*N],rt[N]; 27 char s[5]; 28  29 il int gi(){ 30     RG int x=0,q=1; RG char ch=getchar(); while ((ch<0 || ch>9) && ch!=-) ch=getchar(); 31     if (ch==-) q=-1,ch=getchar(); while (ch>=0 && ch<=9) x=x*10+ch-48,ch=getchar(); return q*x; 32 } 33  34 il void insert(RG int from,RG int to){ g[++num]=(edge){head[from],to},head[from]=num; return; } 35  36 il void dfs1(RG int x,RG int p){ 37     fa[x]=p,dep[x]=dep[p]+1,sz[x]=1; RG int v; 38     for (RG int i=head[x];i;i=g[i].nt){ 39     v=g[i].to; if (v==p) continue; 40     dfs1(v,x); sz[x]+=sz[v]; 41     if (sz[son[x]]<=sz[v]) son[x]=v; 42     } 43     return; 44 } 45  46 il void dfs2(RG int x,RG int p,RG int anc){ 47     top[x]=anc,tid[x]=++cnt; RG int v; 48     if (son[x]) dfs2(son[x],x,anc); 49     for (RG int i=head[x];i;i=g[i].nt){ 50     v=g[i].to; if (v==p || v==son[x]) continue; 51     dfs2(v,x,v); 52     } 53     return; 54 } 55  56 il void update(RG int &x,RG int l,RG int r,RG int p,RG int v){ 57     if (!x) x=++ccnt; 58     if (l==r){ sum[x]=mx[x]=v; return; } 59     RG int mid=(l+r)>>1; 60     if (p<=mid) update(ls[x],l,mid,p,v); 61     else update(rs[x],mid+1,r,p,v); 62     sum[x]=sum[ls[x]]+sum[rs[x]]; 63     mx[x]=max(mx[ls[x]],mx[rs[x]]); 64     return; 65 } 66  67 il int querysum(RG int x,RG int l,RG int r,RG int xl,RG int xr){ 68     if (xl<=l && r<=xr) return sum[x]; 69     if (!x) return 0; RG int mid=(l+r)>>1; 70     if (xr<=mid) return querysum(ls[x],l,mid,xl,xr); 71     else if (xl>mid) return querysum(rs[x],mid+1,r,xl,xr); 72     else return querysum(ls[x],l,mid,xl,mid)+querysum(rs[x],mid+1,r,mid+1,xr); 73 } 74  75 il int querymax(RG int x,RG int l,RG int r,RG int xl,RG int xr){ 76     if (xl<=l && r<=xr) return mx[x]; 77     if (!x) return 0; RG int mid=(l+r)>>1; 78     if (xr<=mid) return querymax(ls[x],l,mid,xl,xr); 79     else if (xl>mid) return querymax(rs[x],mid+1,r,xl,xr); 80     else return max(querymax(ls[x],l,mid,xl,mid),querymax(rs[x],mid+1,r,mid+1,xr)); 81 } 82  83 il int Querysum(RG int c,RG int x,RG int y){ 84     RG int res=0; 85     while (top[x]!=top[y]){ 86     if (dep[top[x]]<dep[top[y]]) swap(x,y); 87     res+=querysum(rt[c],1,n,tid[top[x]],tid[x]); 88     x=fa[top[x]]; 89     } 90     if (dep[x]>dep[y]) swap(x,y); 91     res+=querysum(rt[c],1,n,tid[x],tid[y]); 92     return res; 93 } 94  95 il int Querymax(RG int c,RG int x,RG int y){ 96     RG int res=0; 97     while (top[x]!=top[y]){ 98     if (dep[top[x]]<dep[top[y]]) swap(x,y); 99     res=max(res,querymax(rt[c],1,n,tid[top[x]],tid[x]));100     x=fa[top[x]];101     }102     if (dep[x]>dep[y]) swap(x,y);103     res=max(res,querymax(rt[c],1,n,tid[x],tid[y]));104     return res;105 }106 107 il void work(){108     n=gi(),Q=gi(); for (RG int i=1;i<=n;++i) w[i]=gi(),c[i]=gi();109     for (RG int i=1,u,v;i<n;++i) u=gi(),v=gi(),insert(u,v),insert(v,u);110     dfs1(1,0),dfs2(1,0,1); for (RG int i=1;i<=n;++i) update(rt[c[i]],1,n,tid[i],w[i]);111     for (RG int i=1,x,y;i<=Q;++i){112     scanf("%s",s); x=gi(),y=gi();113     if (s[0]==C && s[1]==C){114         update(rt[c[x]],1,n,tid[x],0);115         update(rt[y],1,n,tid[x],w[x]);116         c[x]=y;117     }118     if (s[0]==C && s[1]==W)119         update(rt[c[x]],1,n,tid[x],y),w[x]=y;120     if (s[0]==Q && s[1]==S)121         printf("%d\n",Querysum(c[y],x,y));122     if (s[0]==Q && s[1]==M)123         printf("%d\n",Querymax(c[y],x,y));124     }125     return;126 }127 128 int main(){129     File("travel");130     work();131     return 0;132 }

 

bzoj3531 [Sdoi2014]旅行