BZOJ-3531 旅行 树链剖分+动态开点线段树
3531: [Sdoi2014]旅行
Time Limit: 20 Sec Memory Limit: 512 MB
Submit: 1097 Solved: 517
[Submit][Status][Discuss]
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。
Source
Round 1 Day 1
题解：
树链剖分无疑，但是需要一些技巧，同样需要注意一些细节
首先需要对每个宗教建一棵线段树，所以做法需要写技巧，动态的开点，动态的建树，对每个宗教记录三个值，root，ll，rr，root记录的是这个宗教的根，ll是左儿子，rr右儿子；这样每次向下询问的时候，以及向上更新的时候就不再是now<<1 和now<<1|1了
注意了这些问题，剩下的就是些写法上的技巧了。
CC操作：把x城市在原来的宗教中评级置成0，在新的宗教中加入评级，并记录宗教的变化
CW操作：直接修改评级，记录下
QS操作：区间查询和，不断跳跳就好
QM操作：区间查询最大，同上。
code：（自己YY的开点线段树，一个地方写的丑，找男神看了看….YM）
#includevoid dfs_2(int now,int chain)
{
pl[now]=++sz; pre[sz]=now; top[now]=chain;
if (son[now]) dfs_2(son[now],chain);
for (int i=head[now]; i; i=edge[i].next)
if (edge[i].to!=son[now] && edge[i].to!=fa[now])
dfs_2(edge[i].to,edge[i].to);
pr[now]=sz;
}
//--------------------------------------------------------------------------------------
int sum[maxn<<2],maxx[maxn<<2],tot,root[maxn<<2],ll[maxn<<2],rr[maxn<<2];
inline void update(int now)
{
sum[now]=sum[ll[now]]+sum[rr[now]];
if (maxx[ll[now]]>maxx[rr[now]])
maxx[now]=maxx[ll[now]]; else maxx[now]=maxx[rr[now]];
}
void point_change(int &now,int l,int r,int loc,int val)
{
if (!now) now=++tot;
if (l==r) {sum[now]=maxx[now]=val;return;}
int mid=(l+r)>>1;
if (loc<=mid) point_change(ll[now],l,mid,loc,val);
else point_change(rr[now],mid+1,r,loc,val);
update(now);
}
int segment_ask_sum(int now,int l,int r,int L,int R)
{
if (!now) return 0;
if(L<=l && R>=r) return sum[now];
int mid=(l+r)>>1;int ans=0;
if (L<=mid) ans+=segment_ask_sum(ll[now],l,mid,L,R);
if (R>mid) ans+=segment_ask_sum(rr[now],mid+1,r,L,R);
return ans;
}
int segment_ask_max(int now,int l,int r,int L,int R)
{
if (!now) return 0;
if (L<=l && R>=r) return maxx[now];
int mid=(l+r)>>1; int ans=-0x7fffffff;
if (L<=mid) ans=max(ans,segment_ask_max(ll[now],l,mid,L,R));
if (R>mid) ans=max(ans,segment_ask_max(rr[now],mid+1,r,L,R));
return ans;
}
//--------------------------------------------------------------------------------------
void CC(int x,int y)
{
point_change(root[city[x].c],1,n,pl[x],0);
point_change(root[y],1,n,pl[x],city[x].w);
city[x].c=y;
}
void CW(int x,int y)
{
point_change(root[city[x].c],1,n,pl[x],y);
city[x].w=y;
}
void QS(int x,int y)
{
int ans=0; int cc=city[x].c,tpx,tpy;
while (top[x]!=top[y])
{
tpx=top[x],tpy=top[y];
if (deep[tpx]1,n,pl[tpx],pl[x]);
x=fa[tpx];
}
//if (x==y) {printf("%d\n",ans);return;}
if (deep[x]>deep[y]) swap(x,y);
ans+=segment_ask_sum(root[cc],1,n,pl[x],pl[y]);
printf("%d\n",ans);
}
void QM(int x,int y)
{
int ans=-0x7fffffff; int cc=city[x].c,tpx,tpy;
while (top[x]!=top[y])
{
tpx=top[x],tpy=top[y];
if (deep[tpx]1,n,pl[tpx],pl[x]));
x=fa[tpx];
}
if (deep[x]>deep[y]) swap(x,y);
ans=max(ans,segment_ask_max(root[cc],1,n,pl[x],pl[y]));
printf("%d\n",ans);
}
//--------------------------------------------------------------------------------------
int main()
{
n=read(),q=read();
for (int i=1; i<=n; i++) city[i].w=read(),city[i].c=read();
for (int u,v,i=1; i<=n-1; i++)
u=read(),v=read(),insert(u,v);
dfs_1(1);dfs_2(1,1);
for (int i=1; i<=n; i++) point_change(root[city[i].c],1,n,pl[i],city[i].w);
for (int i=1; i<=q; i++)
{
char opt[10]; scanf("%s",opt); int x=read(),y=read();
switch (opt[1])
{
case 'C': CC(x,y); break;
case 'S': QS(x,y); break;
case 'M': QM(x,y); break;
case 'W': CW(x,y); break;
}
}
return 0;
}