【bzoj】4538: [Hnoi2016]网络

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=4538

---

维护一个数据结构支持对于一颗树的操作，需要支持：

1.对于树上的一条路径上的每个点上放一个值。

2.撤销某次操作的路劲放。

3.查询除了经过这个点的路径的最大值。

---

往一个路径上丢值相当于往不经过条路径的所有点上丢值。

用一个树链剖分即可维护，对于操作区间取反。

直接查询单点最大值即可。

为了维护单点最大值，线段树中的每一个点对应两个堆，用于维护插入誉删除。

防止爆空间，所以标记永久化即可。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<queue>
using namespace std;
#define maxn 400100
#define llg int
#define yyj(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
llg n,m,tail,deep[maxn],top[maxn],size[maxn],son[maxn],dad[maxn],id[maxn],dfsn,T,type;
vector<llg>a[maxn];

struct node
{
    priority_queue<llg>q1,q2;
    
    void push1(llg x){ q1.push(x);}
    void push2(llg x){ q2.push(x);}

    llg top()
    {
        llg anss=-1;
        while (!q2.empty() && (q2.top()==q1.top())) q1.pop(),q2.pop();
        if (!q1.empty()) anss=q1.top();
        return anss;
    }
}val[maxn*4];

struct data
{
    llg l,r;
}dl[maxn];

struct data_
{
    llg x,y,z;
}ask[maxn];

bool cmp(const data&a,const data&b) {return a.l<b.l;}

void init()
{
    llg x,y;
    cin>>n>>m;
    for (llg i=1;i<n;i++)
    {
        scanf("%d%d",&x,&y);
        a[x].push_back(y),a[y].push_back(x);
    }    
}

void dfs_1(llg x,llg fa)
{
    size[x]=1;
    llg w=a[x].size(),v;
    for (llg i=0;i<w;i++)
    {
        v=a[x][i];
        if (v==fa) continue;
        deep[v]=deep[x]+1;
        dfs_1(v,x);
        dad[v]=x; size[x]+=size[v];  
        if (size[v]>size[son[x]]) son[x]=v;
    }
}

void dfs_2(llg x,llg fa)
{
    id[x]=++dfsn;
    if (son[x]) {top[son[x]]=top[x]; dfs_2(son[x],x);}
    llg w=a[x].size(),v;
    for (llg i=0;i<w;i++)
    {
        v=a[x][i];
        if (v==fa || v==son[x]) continue;
        top[v]=v; dfs_2(v,x);
    }
}

void add(llg o,llg l,llg r,llg L,llg R,llg V)
{
    if (l>=L && r<=R)
    {
        if (T) val[o].push1(V);
        else val[o].push2(V);
        return ;
    }
    llg lc=o<<1,rc=o<<1|1,mid=(l+r)>>1;
    if (L<=mid) add(lc,l,mid,L,R,V);
    if (R>mid) add(rc,mid+1,r,L,R,V);
}

void solve(llg x,llg y,llg V)
{
    llg f1=top[x],f2=top[y];
    tail=0;
    while (f1!=f2)
    {
        if (deep[f1]<deep[f2]) swap(x,y),swap(f1,f2);
        dl[++tail].l=id[f1]; dl[tail].r=id[x];
        x=dad[f1]; f1=top[x];
    }
    if (deep[x]<deep[y]) swap(x,y);
    dl[++tail].l=id[y],dl[tail].r=id[x];
    sort(dl+1,dl+tail+1,cmp);
    llg l=0,r,last=0;
    for (llg i=1;i<=tail;i++)
    {
        l=last+1; r=dl[i].l-1;
        if (l<=r) add(1,1,n,l,r,V);
        last=dl[i].r;
    }
    l=last+1; r=n;
    if (l<=r) add(1,1,n,l,r,V);
}

llg maxl_(llg o,llg l,llg r,llg wz)
{
    llg ans=val[o].top();
    if (l==r) return ans;
    llg mid=(l+r)>>1,lc=o<<1,rc=o<<1|1;
    if (wz<=mid) ans=max(ans,maxl_(lc,l,mid,wz));
    if (wz>mid) ans=max(ans,maxl_(rc,mid+1,r,wz));
    return ans;
}

int main()
{
    yyj("a");
    init();
    dfs_1(1,-1);
    top[1]=1;
    dfs_2(1,-1);
    for (llg o=1;o<=m;o++)
    {
        scanf("%d",&type);
        if (type==0)
        {
            scanf("%d%d%d",&ask[o].x,&ask[o].y,&ask[o].z);
            T=1;
            solve(ask[o].x,ask[o].y,ask[o].z);
        }
        if (type==1)
        {
            llg x;
            scanf("%d",&x);
            T=0;
            solve(ask[x].x,ask[x].y,ask[x].z);
        }
        if (type==2)
        {
            llg x;
            scanf("%d",&x);
            printf("%d
",maxl_(1,1,n,id[x]));
        }
    }
    return 0;
}

---