• 【BZOJ-1984】月下“毛景树” 树链剖分


    1984: 月下“毛景树”

    Time Limit: 20 Sec  Memory Limit: 64 MB
    Submit: 1314  Solved: 416
    [Submit][Status][Discuss]

    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个。

    HINT

    Source

    Solution

    很容易想到树链剖分,但此题权值在边上,所以思路应有所变化,即 边权下放至点 

    具体实现起来并不是很难。 把边权下放到较深的点,然后进行正常的链剖即可。

    维护两个标记,变换和增加,变换操作的时候,同时应该把当前的增加标记置成0,而增加操作时,则无需处理变换标记。

    在区间修改的时候,需要多进行一次判断;这是不同于普通点权链剖的地方。

    Code

    #include<iostream>
    #include<cstdio>
    #include<cstring>
    #include<algorithm>
    #include<cmath>
    using namespace std;
    inline int read()
    {
        int x=0,f=1;char ch=getchar();
        while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
        while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}
        return x*f;
    }
    #define maxn 400010
    int id[maxn<<2],va[maxn<<2],uu[maxn],vv[maxn],ww[maxn];
    int n;struct data{int next,to,w;}edge[maxn<<1];int head[maxn],cnt;
    void add(int u,int v,int w){cnt++;edge[cnt].w=w;edge[cnt].to=v;edge[cnt].next=head[u];head[u]=cnt;}
    void insert(int u,int v,int w){add(u,v,w); add(v,u,w);}
    //--------------------------------------------------------------------------------------
    int deep[maxn],fa[maxn],son[maxn],size[maxn],pl[maxn],sz,pr[maxn],top[maxn],pre[maxn];
    void dfs_1(int now)
    {
        size[now]=1;
        for (int i=head[now]; i; i=edge[i].next)
            if (edge[i].to!=fa[now])
                {
                    fa[edge[i].to]=now;
                    deep[edge[i].to]=deep[now]+1;
                    dfs_1(edge[i].to);
                    if (size[son[now]]<size[edge[i].to]) son[now]=edge[i].to;
                    size[now]+=size[edge[i].to];
                }
    }
    void 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 tree[maxn<<2],delc[maxn<<2],dela[maxn<<2];
    inline void update(int now){tree[now]=max(tree[now<<1],tree[now<<1|1]);}
    inline void pushdown(int now)
    {
        if (delc[now]!=-1)
            {
                int cc=delc[now]; delc[now]=-1;
                delc[now<<1]=cc; dela[now<<1]=0; tree[now<<1]=cc;
                delc[now<<1|1]=cc; dela[now<<1|1]=0; tree[now<<1|1]=cc;
            }
        if (dela[now])
            {
                int aa=dela[now]; dela[now]=0;
                dela[now<<1]+=aa; tree[now<<1]+=aa;
                dela[now<<1|1]+=aa; tree[now<<1|1]+=aa;
            }
    }
    void build(int now,int l,int r)
    {
        delc[now]=-1; dela[now]=0;
        if (l==r) {tree[now]=va[l];return;}
        int mid=(l+r)>>1;
        build(now<<1,l,mid); build(now<<1|1,mid+1,r);
        update(now);
    }
    void point_change(int now,int l,int r,int loc,int val)
    {
        if (l==r) {tree[now]=val;return;}
        int mid=(l+r)>>1;pushdown(now);
        if (loc<=mid) point_change(now<<1,l,mid,loc,val);
        else point_change(now<<1|1,mid+1,r,loc,val);  
        update(now);
    }
    void segment_change(int now,int l,int r,int L,int R,int val)
    {
        if (L<=l && R>=r) {tree[now]=val;delc[now]=val;dela[now]=0;return;}
        int mid=(l+r)>>1;pushdown(now);
        if (L<=mid) segment_change(now<<1,l,mid,L,R,val);
        if (R>mid) segment_change(now<<1|1,mid+1,r,L,R,val);
        update(now);
    }
    void segment_add(int now,int l,int r,int L,int R,int val)
    {
        if (L<=l && R>=r) {tree[now]+=val;dela[now]+=val;return;}
        int mid=(l+r)>>1;pushdown(now);
        if (L<=mid) segment_add(now<<1,l,mid,L,R,val);
        if (R>mid) segment_add(now<<1|1,mid+1,r,L,R,val);
        update(now);
    }
    int segment_ask(int now,int l,int r,int L,int R)
    {
        if (L<=l && R>=r) return tree[now];
        int mid=(l+r)>>1; int ans=0;pushdown(now);
        if (L<=mid) ans=max(segment_ask(now<<1,l,mid,L,R),ans);
        if (R>mid) ans=max(ans,segment_ask(now<<1|1,mid+1,r,L,R));
        return ans;
    }
    //--------------------------------------------------------------------------------------
    void Change(int x,int val)
    {
        point_change(1,1,n,pl[id[x]],val);
    }
    void Cover(int x,int y,int val)
    {
        while (top[x]!=top[y])
            {
                if (deep[top[x]]<deep[top[y]]) swap(x,y);
                segment_change(1,1,n,pl[top[x]],pl[x],val);
                x=fa[top[x]];
            }
        if (deep[x]>deep[y]) swap(x,y);
        if (x!=y) segment_change(1,1,n,pl[x]+1,pl[y],val);//这是有别于平常的判断
    }
    void Add(int x,int y,int val)
    {
        while (top[x]!=top[y])
            {
                if (deep[top[x]]<deep[top[y]]) swap(x,y);
                segment_add(1,1,n,pl[top[x]],pl[x],val);
                x=fa[top[x]];
            }
        if (deep[x]>deep[y]) swap(x,y);
        if (x!=y) segment_add(1,1,n,pl[x]+1,pl[y],val);
    }
    void Max(int x,int y)
    {
        int ans=0;
        while (top[x]!=top[y])
            {
                if (deep[top[x]]<deep[top[y]]) swap(x,y);
                ans=max(ans,segment_ask(1,1,n,pl[top[x]],pl[x]));
                x=fa[top[x]];
            }
        if (deep[x]>deep[y]) swap(x,y);
        if (x!=y) ans=max(ans,segment_ask(1,1,n,pl[x]+1,pl[y]));
        printf("%d
    ",ans);
    }
    //--------------------------------------------------------------------------------------
    int main()
    {
    //  freopen("msn.in","r",stdin);
    //  freopen("msn.out","w",stdout);
        n=read();
        for (int i=1; i<=n-1; i++)
            uu[i]=read(),vv[i]=read(),ww[i]=read(),insert(uu[i],vv[i],ww[i]);
        dfs_1(1); dfs_2(1,1);
        for (int i=1; i<=n-1; i++)
            if (deep[uu[i]]<deep[vv[i]]) id[i]=vv[i]; else id[i]=uu[i];
        for (int i=1; i<=n-1; i++) va[pl[id[i]]]=ww[i];//边权下放至点,记录所下放的点的id
        build(1,1,n);
        while (true)
            {
                char opt[10]; scanf("%s",opt);if(opt[1]=='t')break;
                int u=read(),v=read(),w;
                switch (opt[1])
                    {
                        case 'a': Max(u,v);break;
                        case 'o': w=read();Cover(u,v,w);break;
                        case 'd': w=read();Add(u,v,w);break;
                        case 'h': Change(u,v);break;
                    }
            }
        return 0;
    }

    找DCrusher和Char哥帮我调过...无果,结果还得靠自己...然后发现又是二逼错误...

  • 相关阅读:
    解决android运行速度慢的问题
    nginx配合tomcat、resin等java应用服务器提供java支持(转)
    部分面向对象习题,留给学生做不错!!!
    该内存不能为“read”或“written”
    ORACLE数据库管理系统体系结构详解
    Linux下设置ORACLE自启动
    Oracle静默安装模板
    python安装步骤
    Oracle删除重复索引
    windows操作系统下python IDLE安装
  • 原文地址:https://www.cnblogs.com/DaD3zZ-Beyonder/p/5334131.html
Copyright © 2020-2023  润新知