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 }

 

posted @ 2017-03-26 20:58  wfj_2048  阅读(132)  评论(0)  编辑  收藏  举报