1495 -- 【BZOJ2157】旅游
Description
Ray 乐忠于旅游,这次他来到了T 城。T 城是一个水上城市,一共有 N 个景点,有些景点之间会用一座桥连接。为了方便游客到达每个景点但又为了节约成本,T 城的任意两个景点之间有且只有一条路径。换句话说, T 城中只有N ? 1 座桥。Ray 发现,有些桥上可以看到美丽的景色,让人心情愉悦,但有些桥狭窄泥泞,令人烦躁。于是,他给每座桥定义一个愉悦度w,也就是说,Ray 经过这座桥会增加w 的愉悦度,这或许是正的也可能是负的。有时,Ray 看待同一座桥的心情也会发生改变。现在,Ray 想让你帮他计算从u 景点到v 景点能获得的总愉悦度。有时,他还想知道某段路上最美丽的桥所提供的最大愉悦度,或是某段路上最糟糕的一座桥提供的最低愉悦度。
Input
输入的第一行包含一个整数N,表示T 城中的景点个数。景点编号为 0...N ? 1。
接下来N ? 1 行,每行三个整数u、v 和w,表示有一条u 到v,使 Ray 愉悦度增加w 的桥。桥的编号为1...N ? 1。|w| <= 1000。
输入的第N + 1 行包含一个整数M,表示Ray 的操作数目。
接下来有M 行,每行描述了一个操作,操作有如下五种形式: C i w,表示Ray 对于经过第i 座桥的愉悦度变成了w。 N u v,表示Ray 对于经过景点u 到v 的路径上的每一座桥的愉悦度都变成原来的相反数。 SUM u v,表示询问从景点u 到v 所获得的总愉悦度。 MAX u v,表示询问从景点u 到v 的路径上的所有桥中某一座桥所提供的最大愉悦度。 MIN u v,表示询问从景点u 到v 的路径上的所有桥中某一座桥所提供的最小愉悦度。测试数据保证,任意时刻,Ray 对于经过每一座桥的愉悦度的绝对值小于等于1000。
Output
对于每一个询问(操作S、MAX 和MIN),输出答案。
Sample Input
30
1 1
1 2 2
8
SUM 0 2
MAX 0 2
N 0 1
SUM 0 2
MIN 0 2
C 1 3
SUM 0 2
MAX 0 2
Sample Output
3
2
1
-1
5
3
Hint
Description
Ray 乐忠于旅游,这次他来到了T 城。T 城是一个水上城市,一共有 N 个景点,有些景点之间会用一座桥连接。为了方便游客到达每个景点但又为了节约成本,T 城的任意两个景点之间有且只有一条路径。换句话说, T 城中只有N ? 1 座桥。Ray 发现,有些桥上可以看到美丽的景色,让人心情愉悦,但有些桥狭窄泥泞,令人烦躁。于是,他给每座桥定义一个愉悦度w,也就是说,Ray 经过这座桥会增加w 的愉悦度,这或许是正的也可能是负的。有时,Ray 看待同一座桥的心情也会发生改变。现在,Ray 想让你帮他计算从u 景点到v 景点能获得的总愉悦度。有时,他还想知道某段路上最美丽的桥所提供的最大愉悦度,或是某段路上最糟糕的一座桥提供的最低愉悦度。
Input
输入的第一行包含一个整数N,表示T 城中的景点个数。景点编号为 0...N ? 1。
接下来N ? 1 行,每行三个整数u、v 和w,表示有一条u 到v,使 Ray 愉悦度增加w 的桥。桥的编号为1...N ? 1。|w| <= 1000。
输入的第N + 1 行包含一个整数M,表示Ray 的操作数目。
接下来有M 行,每行描述了一个操作,操作有如下五种形式: C i w,表示Ray 对于经过第i 座桥的愉悦度变成了w。 N u v,表示Ray 对于经过景点u 到v 的路径上的每一座桥的愉悦度都变成原来的相反数。 SUM u v,表示询问从景点u 到v 所获得的总愉悦度。 MAX u v,表示询问从景点u 到v 的路径上的所有桥中某一座桥所提供的最大愉悦度。 MIN u v,表示询问从景点u 到v 的路径上的所有桥中某一座桥所提供的最小愉悦度。测试数据保证,任意时刻,Ray 对于经过每一座桥的愉悦度的绝对值小于等于1000。
Output
对于每一个询问(操作S、MAX 和MIN),输出答案。
Sample Input
30
1 1
1 2 2
8
SUM 0 2
MAX 0 2
N 0 1
SUM 0 2
MIN 0 2
C 1 3
SUM 0 2
MAX 0 2
Sample Output
3
2
1
-1
5
3
Hint
一共有10 个数据,对于第i (1 <= i <= 10) 个数据, N = M = i * 2000。
树链剖分模板题,都是很基础的操作..
好久没写过树剖了,调了半天,发现一开始size没设为1,还少写了一句DFS2,有负边权maxx初值竟然设为0...
一脸尴尬,本来准备快速AC的,惨遭打脸。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
//#include<cmath>
#define L(x) (x<<1)
#define R(x) (x<<1|1)
#define INF 0x3fffffff
using namespace std;
struct seg
{
int l,r,maxx,minn,sum,lazy;
}tree[100005];
struct TREE
{
int x,y,val;
}T[50005];
struct node
{
int to,next,val;
}w[1000005];
int h[1000005]={0},cnt=0,dep[200005]={0},fa[1000005]={0},size[1000005]={0},top[1000005]={0},tid[1000005]={0},tim=0,n,son[100005]={0};
void Addarc(int x,int y,int val)
{
cnt++;w[cnt].to=y;w[cnt].next=h[x];h[x]=cnt;
}
void DFS1(int x,int f,int d)
{
dep[x]=d;
fa[x]=f;
size[x]=1;
for(int i=h[x];i;i=w[i].next)
{
int to=w[i].to;
if(to==f)continue;
DFS1(to,x,d+1);
size[x]+=size[to];
if(son[x]==0||(size[son[x]]<size[to]))son[x]=to;
}
}
void DFS2(int x,int tp)
{
top[x]=tp;
tid[x]=++tim;
if(son[x]==0)return;
DFS2(son[x],tp);
for(int i=h[x];i;i=w[i].next)
{
int to=w[i].to;
if(to==fa[x]||to==son[x])continue;
DFS2(to,to);
}
}
void Pushup(int root)
{
tree[root].sum=tree[L(root)].sum+tree[R(root)].sum;
tree[root].maxx=max(tree[L(root)].maxx,tree[R(root)].maxx);
tree[root].minn=min(tree[L(root)].minn,tree[R(root)].minn);
}
void Pushdown(int root)
{
if(tree[root].lazy)
{
tree[root].lazy=0;
tree[L(root)].lazy^=1;
tree[R(root)].lazy^=1;
tree[L(root)].sum*=-1;
tree[R(root)].sum*=-1;
tree[L(root)].maxx*=-1;tree[L(root)].minn*=-1;swap(tree[L(root)].maxx,tree[L(root)].minn);
tree[R(root)].maxx*=-1;tree[R(root)].minn*=-1;swap(tree[R(root)].maxx,tree[R(root)].minn);
}
}
void Build(int l,int r,int root)
{
tree[root].l=l;tree[root].r=r;
tree[root].maxx=-INF;tree[root].minn=INF;tree[root].lazy=0;tree[root].sum=0;
if(l==r)return;
int mid=(l+r)>>1;
Pushdown(root);
Build(l,mid,L(root));
Build(mid+1,r,R(root));
Pushup(root);
}
void Change(int x,int data,int root)
{
if(tree[root].l>x||tree[root].r<x)return;
if(tree[root].l==tree[root].r&&tree[root].l==x)
{
tree[root].maxx=data;
tree[root].minn=data;
tree[root].sum=data;
return;
}
Pushdown(root);
Change(x,data,L(root));
Change(x,data,R(root));
Pushup(root);
}
void Neg(int l,int r,int root)
{
if(tree[root].l>r||tree[root].r<l)return;
if(l<=tree[root].l&&tree[root].r<=r)
{
tree[root].lazy^=1;
tree[root].maxx*=-1;
tree[root].minn*=-1;
tree[root].sum*=-1;
swap(tree[root].minn,tree[root].maxx);
return;
}
Pushdown(root);
Neg(l,r,L(root));
Neg(l,r,R(root));
Pushup(root);
}
int Ask_sum(int l,int r,int root)
{
if(tree[root].l>r||tree[root].r<l)return 0;
if(l<=tree[root].l&&tree[root].r<=r)return tree[root].sum;
int sum=0;
Pushdown(root);
sum+=Ask_sum(l,r,L(root));
sum+=Ask_sum(l,r,R(root));
Pushup(root);
return sum;
}
int Ask_max(int l,int r,int root)
{
if(tree[root].l>r||tree[root].r<l)return -INF;
if(l<=tree[root].l&&tree[root].r<=r)return tree[root].maxx;
int maxx=-INF;
Pushdown(root);
maxx=max(maxx,Ask_max(l,r,L(root)));
maxx=max(maxx,Ask_max(l,r,R(root)));
Pushup(root);
return maxx;
}
int Ask_min(int l,int r,int root)
{
if(tree[root].l>r||tree[root].r<l)return INF;
if(l<=tree[root].l&&tree[root].r<=r)return tree[root].minn;
int minn=INF;
Pushdown(root);
minn=min(minn,Ask_min(l,r,L(root)));
minn=min(minn,Ask_min(l,r,R(root)));
Pushup(root);
return minn;
}
void Doneg(int x,int y)
{
int t1=top[x],t2=top[y];
while(t1!=t2)
{
if(dep[t1]<dep[t2]){swap(x,y);swap(t1,t2);}
Neg(tid[t1],tid[x],1);
x=fa[t1];t1=top[x];
}
if(x==y)return;
if(dep[x]>dep[y])swap(x,y);
Neg(tid[son[x]],tid[y],1);
}
int Getsum(int x,int y)
{
int t1=top[x],t2=top[y],ans=0;
while(t1!=t2)
{
if(dep[t1]<dep[t2]){swap(x,y);swap(t1,t2);}
ans+=Ask_sum(tid[t1],tid[x],1);
x=fa[t1];t1=top[x];
}
if(x==y)return ans;
if(dep[x]>dep[y])swap(x,y);
ans+=Ask_sum(tid[son[x]],tid[y],1);
return ans;
}
int Getmax(int x,int y)
{
int t1=top[x],t2=top[y],maxx=-INF;
while(t1!=t2)
{
if(dep[t1]<dep[t2]){swap(x,y);swap(t1,t2);}
maxx=max(maxx,Ask_max(tid[t1],tid[x],1));
x=fa[t1];t1=top[x];
}
if(x==y)return maxx;
if(dep[x]>dep[y])swap(x,y);
maxx=max(maxx,Ask_max(tid[son[x]],tid[y],1));
return maxx;
}
int Getmin(int x,int y)
{
int t1=top[x],t2=top[y],minn=INF;
while(t1!=t2)
{
if(dep[t1]<dep[t2]){swap(x,y);swap(t1,t2);}
minn=min(minn,Ask_min(tid[t1],tid[x],1));
x=fa[t1];t1=top[x];
}
if(x==y)return minn;
if(dep[x]>dep[y])swap(x,y);
minn=min(minn,Ask_min(tid[son[x]],tid[y],1));
return minn;
}
void init()
{
scanf("%d",&n);
for(int i=1;i<n;i++)
{
int x,y,z;
scanf("%d%d%d",&x,&y,&z);
x++;y++;
T[i].x=x;T[i].y=y;T[i].val=z;
Addarc(T[i].x,T[i].y,T[i].val);
Addarc(T[i].y,T[i].x,T[i].val);
}
DFS1(1,0,0);
DFS2(1,1);
Build(1,n,1);
for(int i=1;i<n;i++)
{
if(dep[T[i].x]>dep[T[i].y])swap(T[i].x,T[i].y);
Change(tid[T[i].y],T[i].val,1);
}
}
void work()
{
int m,a,b;
string s;
scanf("%d",&m);
for(int i=1;i<=m;i++)
{
cin>>s;
scanf("%d%d",&a,&b);
if(s=="C")Change(tid[T[a].y],b,1);
if(s=="N")Doneg(++a,++b);
if(s=="SUM")printf("%d\n",Getsum(++a,++b));
if(s=="MAX")printf("%d\n",Getmax(++a,++b));
if(s=="MIN")printf("%d\n",Getmin(++a,++b));
}
}
int main()
{
init();
work();
return 0;
}
/*
3
0 1 1
1 2 2
8
SUM 0 2
MAX 0 2
N 0 1
SUM 0 2
MIN 0 2
C 1 3
SUM 0 2
MAX 0 2
*/
稍微注意下这个细节,因为是边权所以记录在儿子头上,所以最后一步修改是xxx(tid[son[x]],tid[y],1)
还有就是题目是边权的情况下不要忘了x==y返回这个条件。