bzoj1984 月下“毛景树”

http://www.elijahqi.win/archives/811
Description

毛毛虫经过及时的变形,最终逃过的一劫,离开了菜妈的菜园。 毛毛虫经过千山万水,历尽千辛万苦,最后来到了小小的绍兴一中的校园里。爬啊爬~爬啊爬毛毛虫爬到了一颗小小的“毛景树”下面,发现树上长着他最爱吃的毛毛果~ “毛景树”上有N个节点和N-1条树枝,但节点上是没有毛毛果的,毛毛果都是长在树枝上的。但是这棵“毛景树”有着神奇的魔力,他能改变树枝上毛毛果的个数:  Change k w:将第k条树枝上毛毛果的个数改变为w个。  Cover u v w:将节点u与节点v之间的树枝上毛毛果的个数都改变为w个。  Add u v w:将节点u与节点v之间的树枝上毛毛果的个数都增加w个。 由于毛毛虫很贪,于是他会有如下询问:  Max u v:询问节点u与节点v之间树枝上毛毛果个数最多有多少个。

Input

第一行一个正整数N。 接下来N-1行,每行三个正整数Ui,Vi和Wi,第i+1行描述第i条树枝。表示第i条树枝连接节点Ui和节点Vi,树枝上有Wi个毛毛果。 接下来是操作和询问,以“Stop”结束。

Output

对于毛毛虫的每个询问操作,输出一个答案。

Sample Input

4
1 2 8
1 3 7
3 4 9
Max 2 4
Cover 2 4 5
Add 1 4 10
Change 1 16
Max 2 4
Stop
Sample Output

9
16
【Data Range】
1<=N<=100,000,操作+询问数目不超过100,000。
保证在任意时刻,所有树枝上毛毛果的个数都不会超过10^9个。

不得不说写程序看出个人的十分不严谨,对于add和change操作我觉得比较像,直接抄过来,结果两个add操作,有一个忘改了,查错又是良久

大致说一下这个边权树链剖分和点权的区别吧
做边权的时候,采取的策略是将 我的每个边权下放到子节点去,这样的话我的root就没有值了 不过没关系,对于每两个节点的的lca我们也不应该处理,因为这个点的值不属于这条链上的边权
我们采取的策略就是当跳不动的时候,说明在一条链上了,那么我们不妨把id小的加一做,这样就直接跳过了lca
另外这题有一个修改操作,是修改边,而其他操作都是点与点之间,我多记录一个son1[]数组解决这个问题

#include<cstdio>
#define N 110000
inline int read(){
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9') {if (ch=='-') f=-1;ch=getchar();}
    while (ch<='9'&&ch>='0') {x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
struct node{
    int y,next,z;
}data[N<<1];
struct node1{
    int left,right,l,r,max,lazy,add;
}tree[N<<2];
inline void swap(int &x,int &y){
    x^=y;y^=x;x^=y;
}
int size[N],a[N],h[N],son[N],w[N],id[N],num,dep[N],fa[N],tp[N],n,son1[N],root;
void dfs1(int x){
    size[x]=1;
    for (int i=h[x];i;i=data[i].next){
        int y=data[i].y,z=data[i].z;
        if (fa[x]==y) continue;a[y]=z;son1[i>>1]=y;
        dep[y]=dep[x]+1;fa[y]=x;dfs1(y),size[x]+=size[y];
        if (size[y]>size[son[x]]) son[x]=y;
    }
}
void dfs2(int x,int top){
    id[x]=++num;tp[x]=top;w[num]=a[x];
    if (son[x]) dfs2(son[x],top);
    for (int i=h[x];i;i=data[i].next){
        int y=data[i].y;
        if (fa[x]==y||son[x]==y) continue;
        dfs2(y,y);
    }
}
inline int max(int x,int y){
    return x>y?x:y;
}
inline void update(int x){
    int l=tree[x].left,r=tree[x].right;
    tree[x].max=max(tree[l].max,tree[r].max);
}
inline void pushdown(int x){
    if (!tree[x].add&&!tree[x].lazy) return;
    int l=tree[x].left,r=tree[x].right;
    if (tree[x].lazy){
        tree[l].max=tree[r].max=tree[l].lazy=tree[r].lazy=tree[x].lazy;
        tree[l].add=tree[r].add=0;
        tree[x].lazy=0;
    }
    if (tree[x].add){
        tree[l].max+=tree[x].add;tree[r].max+=tree[x].add;
        tree[l].add+=tree[x].add;tree[r].add+=tree[x].add;
        tree[x].add=0;
    }
}
void build(int &x,int l,int r){
    x=++num;
    tree[x].l=l;tree[x].r=r;
    if (l==r) {tree[x].max=w[l];return;}
    int mid=l+r>>1;
    build(tree[x].left,l,mid);build(tree[x].right,mid+1,r);
    update(x);
}
void change(int x,int l,int r,int v){
    if (l<=tree[x].l&&r>=tree[x].r){
        tree[x].max=tree[x].lazy=v;tree[x].add=0;return;
    }
    int mid=(tree[x].l+tree[x].r)>>1;pushdown(x);
    if (l<=mid) change(tree[x].left,l,r,v);
    if (r>mid) change(tree[x].right,l,r,v);
    update(x);
}
void add(int x,int l,int r,int v){
    if (l<=tree[x].l&&r>=tree[x].r){
        tree[x].max+=v;tree[x].add+=v;return;
    }
    int mid=(tree[x].l+tree[x].r)>>1;pushdown(x);
    if (l<=mid) add(tree[x].left,l,r,v);
    if (r>mid) add(tree[x].right,l,r,v);
    update(x);
}
int query(int x,int l,int r){
    if (l<=tree[x].l&&r>=tree[x].r) return tree[x].max;
    int mid=(tree[x].l+tree[x].r)>>1;pushdown(x);int max1=-1;
    if (l<=mid) max1=max(max1,query(tree[x].left,l,r));
    if (r>mid) max1=max(max1,query(tree[x].right,l,r));
    return max1;
}
int main(){
    freopen("bzoj1984.in","r",stdin);
    n=read();num=1;
    for (int i=1;i<n;++i){
        int x=read(),y=read(),z=read();
        data[++num].y=y;data[num].z=z;data[num].next=h[x];h[x]=num;
        data[++num].y=x;data[num].z=z;data[num].next=h[y];h[y]=num;
    }dep[1]=1;num=0;
    dfs1(1);dfs2(1,1);num=0;
    build(root,1,n);char str1[10];
    while (1){
        scanf("%s",str1);
        if (str1[1]=='t') break;
        if (str1[1]=='o'){
            int x=read(),y=read(),z=read();
            while(tp[x]!=tp[y]){
                if(dep[tp[x]]<dep[tp[y]]) swap(x,y);
                change(root,id[tp[x]],id[x],z);
                x=fa[tp[x]];
            } 
            if (id[x]>id[y]) swap(x,y);
            if (id[x]==id[y]) continue;
            change(root,id[x]+1,id[y],z);
            continue;
        }
        if (str1[1]=='d'){
            int x=read(),y=read(),z=read();
            while(tp[x]!=tp[y]){
                if(dep[tp[x]]<dep[tp[y]]) swap(x,y);
                add(root,id[tp[x]],id[x],z);
                x=fa[tp[x]];
            } 
            if (id[x]>id[y]) swap(x,y);
            if (id[x]==id[y]) continue;
            add(root,id[x]+1,id[y],z);
            continue;
        }
        if (str1[1]=='h'){
            int x=read(),z=read();
            change(root,id[son1[x]],id[son1[x]],z);continue;
        }
        if (str1[1]=='a'){
            int x=read(),y=read();int max1=-1;
            while (tp[x]!=tp[y]){
                if (dep[tp[x]]<dep[tp[y]]) swap(x,y);
                max1=max(max1,query(root,id[tp[x]],id[x]));
                x=fa[tp[x]];
            }
            if (id[x]>id[y]) swap(x,y);
            if (id[x]==id[y]) {printf("%d\n",max1);continue;}
            max1=max(max1,query(root,id[x]+1,id[y]));
            printf("%d\n",max1);
        }
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值