bzoj 1984: 月下“毛景树” 线段树+树链剖分

题意

给出一棵n个节点的树，每条边都有权值，要求资瓷以下操作：
Max x y表示查询x到y之间的最大权值
Cover x y z表示把x到y的权值赋为z
Change x y表示把第x条边的权值变成y
Add x y z表示把x到y的权值加上z
n+m<=100000

分析

一眼树剖+线段树。
把边权下放到两端深度较大的那个点即可。
线段树维护两个标记和最大值。
注意两个点的lca不用处理。

第一次把线段树写到struct里面，发现比写外面要快一丢丢。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define N 100005
using namespace std;

int n,last[N],cnt,pos[N],top[N],size[N],dep[N],fa[N],val[N],sz;
struct edge{int to,next,w;}e[N*2];

struct Sigma_Tree
{

    struct Sigma_tree{int l,r,tag1,tag2,mx;}t[N*5];

    void updata(int d)
    {
        t[d].mx=max(t[d*2].mx,t[d*2+1].mx);
    }

    void pushdown(int d)
    {
        if (t[d].l==t[d].r) return;
        if (t[d].tag1)
        {
            t[d*2].tag2=t[d*2+1].tag2=0;
            t[d*2].tag1=t[d*2+1].tag1=t[d].tag1;
            t[d*2].mx=t[d*2+1].mx=t[d].tag1;
            t[d].tag1=0;
        }
        if (t[d].tag2)
        {
            t[d*2].tag2+=t[d].tag2;t[d*2+1].tag2+=t[d].tag2;
            t[d*2].mx+=t[d].tag2;t[d*2+1].mx+=t[d].tag2;
            t[d].tag2=0;
        }
    }

    void build(int d,int l,int r)
    {
        t[d].l=l;
        t[d].r=r;
        if (l==r) return;
        int mid=(l+r)/2;
        build(d*2,l,mid);
        build(d*2+1,mid+1,r);
    }

    void cover(int d,int x,int y,int z)
    {
        pushdown(d);
        if (t[d].l==x&&t[d].r==y)
        {
            t[d].tag2=0;
            t[d].tag1=t[d].mx=z;
            return;
        }
        int mid=(t[d].l+t[d].r)/2;
        if (y<=mid) cover(d*2,x,y,z);
        else if (x>mid) cover(d*2+1,x,y,z);
        else
        {
            cover(d*2,x,mid,z);
            cover(d*2+1,mid+1,y,z);
        }
        updata(d);
    }

    void add(int d,int x,int y,int z)
    {
        pushdown(d);
        if (t[d].l==x&&t[d].r==y)
        {
            t[d].tag2=z;
            t[d].mx+=z;
            return;
        }
        int mid=(t[d].l+t[d].r)/2;
        if (y<=mid) add(d*2,x,y,z);
        else if (x>mid) add(d*2+1,x,y,z);
        else
        {
            add(d*2,x,mid,z);
            add(d*2+1,mid+1,y,z);
        }
        updata(d);
    }

    int query(int d,int x,int y)
    {
        pushdown(d);
        if (t[d].l==x&&t[d].r==y) return t[d].mx;
        int mid=(t[d].l+t[d].r)/2;
        if (y<=mid) return query(d*2,x,y);
        else if (x>mid) return query(d*2+1,x,y);
        else return max(query(d*2,x,mid),query(d*2+1,mid+1,y));
    }

}tree;

void addedge(int u,int v,int w)
{
    e[++cnt].to=v;e[cnt].w=w;e[cnt].next=last[u];last[u]=cnt;
    e[++cnt].to=u;e[cnt].w=w;e[cnt].next=last[v];last[v]=cnt;
}

void dfs1(int x)
{ 
    dep[x]=dep[fa[x]]+1;size[x]=1;
    for (int i=last[x];i;i=e[i].next)
    {
        if (e[i].to==fa[x]) continue;
        fa[e[i].to]=x;
        val[e[i].to]=e[i].w;
        dfs1(e[i].to);
        size[x]+=size[e[i].to];
    }
}

void dfs2(int x,int chain)
{
    pos[x]=++sz;top[x]=chain;
    int k=0;
    for (int i=last[x];i;i=e[i].next)
        if (e[i].to!=fa[x]&&size[e[i].to]>size[k]) k=e[i].to;
    if (!k) return;
    dfs2(k,chain);
    for (int i=last[x];i;i=e[i].next)
        if (e[i].to!=fa[x]&&e[i].to!=k) dfs2(e[i].to,e[i].to);
}

int query(int x,int y)
{
    int ans=0;
    while (top[x]!=top[y])
    {
        if (dep[top[x]]<dep[top[y]]) swap(x,y);
        ans=max(ans,tree.query(1,pos[top[x]],pos[x]));
        x=fa[top[x]];
    }
    if (x==y) return ans;
    if (dep[x]<dep[y]) swap(x,y);
    ans=max(ans,tree.query(1,pos[y]+1,pos[x]));
    return ans;
}

void cover(int x,int y,int z)
{
    while (top[x]!=top[y])
    {
        if (dep[top[x]]<dep[top[y]]) swap(x,y);
        tree.cover(1,pos[top[x]],pos[x],z);
        x=fa[top[x]];
    }
    if (x==y) return;
    if (dep[x]<dep[y]) swap(x,y);
    tree.cover(1,pos[y]+1,pos[x],z);
}

void add(int x,int y,int z)
{
    while (top[x]!=top[y])
    {
        if (dep[top[x]]<dep[top[y]]) swap(x,y);
        tree.add(1,pos[top[x]],pos[x],z);
        x=fa[top[x]];
    }
    if (x==y) return;
    if (dep[x]<dep[y]) swap(x,y);
    tree.add(1,pos[y]+1,pos[x],z);
}

int main()
{
    scanf("%d",&n);
    for (int i=1;i<n;i++)
    {
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        addedge(x,y,z);
    }
    dfs1(1);
    dfs2(1,1);
    tree.build(1,1,n);
    for (int i=2;i<=n;i++)
        tree.cover(1,pos[i],pos[i],val[i]);
    while (1)
    {
        char ch[10];
        scanf("%s",ch);
        if (ch[0]=='S') break;
        if (ch[0]=='M')
        {
            int x,y;
            scanf("%d%d",&x,&y);
            printf("%d\n",query(x,y));
        }else if (ch[0]=='C'&&ch[1]=='o')
        {
            int x,y,z;
            scanf("%d%d%d",&x,&y,&z);
            cover(x,y,z);
        }else if (ch[0]=='A')
        {
            int x,y,z;
            scanf("%d%d%d",&x,&y,&z);
            add(x,y,z);
        }else
        {
            int x,y;
            scanf("%d%d",&x,&y);
            if (dep[e[x*2].to]>dep[e[x*2-1].to]) tree.cover(1,pos[e[x*2].to],pos[e[x*2].to],y);
            else tree.cover(1,pos[e[x*2-1].to],pos[e[x*2-1].to],y);
        }
    }
    return 0;
}