【bzoj2157】旅游

2015.02.09 15:47 Mon | 3次阅读 | 旧日oi | 固定链接 | 源码

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

巨长……自己看去

题解

很裸的树链剖分,边权赋给深度大的点,操作就那么几种,注意修改的时候是对边修改,需要记录每个边对应的点。时间长不写了,调了好久发现有一个rank【x】写成x了。。
狂吐槽样例……

我的程序

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define maxn 20005
#define inf 0x7fffffff
using namespace std;
int n,m,x,y,w,cnt=-1;
struct edge
{
    int to;
    int next;
    int val;
    int id;
}e[maxn<<1];
int h[maxn],tp;
inline void ae(int u,int v,int w,int t)
{
    e[++tp].to=v;e[tp].next=h[u];
    e[tp].val=w;e[tp].id=t;h[u]=tp;
}
char s[10];
int fa[maxn],dep[maxn],siz[maxn],son[maxn];
int rank[maxn],val[maxn],top[maxn],id[maxn],tid[maxn];
void dfs1(int u,int f,int d)
{
    fa[u]=f;dep[u]=d;siz[u]=1;
    for(int v,i=h[u];i;i=e[i].next)
    {
        v=e[i].to;
        if(v==f) continue;
        val[v]=e[i].val;
        id[e[i].id]=v;
        dfs1(v,u,d+1);
        if(son[u]==-1||siz[son[u]]<siz[v]) son[u]=v;
        siz[u]+=siz[v];
    }
}
void dfs2(int u,int tp)
{
    top[u]=tp;
    rank[u]=++cnt;
    tid[rank[u]]=u;
    if(son[u]==-1) return;
    dfs2(son[u],tp);
    for(int v,i=h[u];i;i=e[i].next)
    {
        v=e[i].to;
        if(v==fa[u]||son[u]==v) continue;
        dfs2(v,v);
    }
}
void tree_cut()
{
    memset(son,-1,sizeof(son));
    dfs1(0,0,0);
    dfs2(0,0);
}
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
int mx[maxn<<2],mn[maxn<<2],sum[maxn<<2];
bool flag[maxn<<2];
inline void push_up(int rt)
{
    sum[rt]=sum[rt<<1]+sum[rt<<1|1];
    mx[rt]=max(mx[rt<<1],mx[rt<<1|1]);
    mn[rt]=min(mn[rt<<1],mn[rt<<1|1]);
}
void rev(int rt)
{
    sum[rt]*=-1;
    mx[rt]*=-1;
    mn[rt]*=-1;
    swap(mx[rt],mn[rt]);
    flag[rt]^=1;
}
void push_down(int rt)
{
    if(flag[rt])
    {
        rev(rt<<1);
        rev(rt<<1|1);
        flag[rt]=0;
    }
}
void build(int l,int r,int rt)
{
    if(l==r)
    {
        mx[rt]=mn[rt]=sum[rt]=val[tid[l]];
        flag[rt]=0;
        return;
    }
    int mid=(l+r)>>1;
    build(lson);
    build(rson);
    push_up(rt);
}
void update(int l,int r,int rt,int x,int y)
{
    if(l==r)
    {
        mx[rt]=mn[rt]=sum[rt]=y;
        flag[rt]=0;
        return;
    }
    push_down(rt);
    int mid=(l+r)>>1;
    if(x<=mid) update(lson,x,y);
    else update(rson,x,y);
    push_up(rt);
}
void reverse(int l,int r,int rt,int L,int R)
{
    if(L<=l&&R>=r)
    {
        rev(rt);
        return;
    }
    push_down(rt);
    int mid=(l+r)>>1;
    if(L<=mid) reverse(lson,L,R);
    if(R>mid) reverse(rson,L,R);
    push_up(rt);
}
inline void do_rev(int x,int y)
{
    if(x==y) return;
    while(top[x]!=top[y])
    {
        if(dep[top[x]]<dep[top[y]]) swap(x,y);
//      for(int i=rank[top[x]];i<=rank[x];i++) cout<<val[tid[i]]<<" "<<tid[i]<<endl;
        reverse(1,n-1,1,rank[top[x]],rank[x]);
        x=fa[top[x]];
    }
    if(x==y) return;
    if(dep[x]<dep[y]) swap(x,y);
    //for(int i=rank[y]+1;i<=rank[x];i++) cout<<val[tid[i]]<<" "<<tid[i]<<endl;
    reverse(1,n-1,1,rank[son[y]],rank[x]);
}
inline int change(int ret,int x,int sig)
{
    if(sig==1) return ret+x;
    else if(sig==2) return max(ret,x);
    else return min(ret,x);
}
int query(int l,int r,int rt,int L,int R,int sig)
{
    int ret=0;
    if(sig==2) ret=-inf;
    else if(sig==3) ret=inf;
    if(L<=l&&R>=r)
    {
        if(sig==1) return sum[rt];
        if(sig==2) return mx[rt];
        return mn[rt];
    }
    push_down(rt);
    int mid=(l+r)>>1;
    if(L<=mid) ret=change(ret,query(lson,L,R,sig),sig);
    if(R>mid) ret=change(ret,query(rson,L,R,sig),sig); 
    push_up(rt);
    return ret;
}
inline int get_ans(int x,int y,int sig)
{
    int ret=0;
    if(sig==2) ret=-inf;
    else if(sig==3) ret=inf;
    while(top[x]!=top[y])
    {
        if(dep[top[x]]<dep[top[y]]) swap(x,y);
        ret=change(ret,query(1,n-1,1,rank[top[x]],rank[x],sig),sig);
        x=fa[top[x]];
    }
    if(x==y) return ret;
    if(dep[x]<dep[y]) swap(x,y);
    //for(int i=rank[y]+1;i<=rank[x];i++) cout<<val[tid[i]]<<" "<<tid[i]<<endl;
    ret=change(ret,query(1,n-1,1,rank[son[y]],rank[x],sig),sig);
    return ret;
}
int main()
{
    //freopen("stdin.in","r",stdin);
    //freopen("me.out","w",stdout);
    cin>>n;
    for(int i=1;i<n;i++)
    {
        scanf("%d%d%d",&x,&y,&w);
        ae(x,y,w,i);
        ae(y,x,w,i);
    }
    tree_cut();
    memset(mx,0xef,sizeof(mx));
    memset(mn,0x7f,sizeof(mn));
    build(1,n-1,1);
    cin>>m;
    for(int i=1;i<=m;i++)
    {
        scanf("%s%d%d",s,&x,&y);
        switch(s[0])
        {
            case 'C':update(1,n-1,1,rank[id[x]],y);break;
            case 'N':do_rev(x,y);break;
            case 'S':printf("%d\n",get_ans(x,y,1));break;
            default: printf("%d\n",s[1]=='A'?get_ans(x,y,2):get_ans(x,y,3));break;
        }
    }
    return 0;
}```