bzoj 4538: [Hnoi2016]网络

题意：

给一棵树，可以将某条路径打上/删除一个标记，询问一个点，所有打了标记的不经过该点的路径中的最大值。

题解：

一开始想到拆成两条链，dfs序一下，但好像只会维护和。
orz，发现是树链剖分+线段树套堆。
于是就可做了，树剖是将不包含的线段标记出来，在线段树上的每个点上开两个优先队列，分别表示加的值和删的值，然后两个一起pop。标记永久化。
还有些玄学地方，具体见bzoj discuss。
code：

#include<queue>#include<cstdio>#include<cstdlib>#include<cstring>#include<iostream>#include<algorithm>using namespace std;struct node{    int y,next;}a[200010];int len=0,last[100010];struct TRnode{    int fa,top,son,tot,dep;}TR[100010];struct trnode{    int lc,rc;    priority_queue<int> q1,q2;}tr[200010];int tot=0,root=0;int ys[100010],z=0;int n,m;void ins(int x,int y){    a[++len].y=y;    a[len].next=last[x];last[x]=len;}void pre_node(int x,int fa){    TR[x].dep=TR[fa].dep+1;TR[x].fa=fa;TR[x].tot=1;    for(int i=last[x];i;i=a[i].next)    {        int y=a[i].y;        if(y==fa) continue;        pre_node(y,x);        if(TR[TR[x].son].tot<=TR[y].tot) TR[x].son=y;        TR[x].tot+=TR[y].tot;    }}void pre_edge(int x,int top){    TR[x].top=top;ys[x]=++z;    if(TR[x].son) pre_edge(TR[x].son,top);    for(int i=last[x];i;i=a[i].next)    {        int y=a[i].y;        if(y==TR[x].fa||y==TR[x].son) continue;        pre_edge(y,y);    }}struct query{    int a,b,v;}q[200010];query t[1000];int num;bool cmp(query a,query b){return a.a<b.a;}void change(int &x,int l,int r,int fl,int fr,int k,int op){    if(!x) x=++tot;    if(l==fl&&r==fr)    {        if(op==1) tr[x].q1.push(k);        else tr[x].q2.push(k);        return;    }    int mid=(l+r)/2;    if(fr<=mid) change(tr[x].lc,l,mid,fl,fr,k,op);    else if(fl>mid) change(tr[x].rc,mid+1,r,fl,fr,k,op);    else change(tr[x].lc,l,mid,fl,mid,k,op),change(tr[x].rc,mid+1,r,mid+1,fr,k,op);}int findans(int x,int l,int r,int k){    if(!x) return -1;    while(!tr[x].q1.empty()&&!tr[x].q2.empty()&&tr[x].q1.top()==tr[x].q2.top())        tr[x].q1.pop(),tr[x].q2.pop();    int mid=(l+r)/2,K;    if(tr[x].q1.empty()) K=-1;    else K=tr[x].q1.top();    if(l==r) return K;    if(k<=mid) return max(K,findans(tr[x].lc,l,mid,k));    return max(K,findans(tr[x].rc,mid+1,r,k));}void solve(int x,int y,int k,int op){    num=0;    int tx=TR[x].top,ty=TR[y].top;    while(tx!=ty)    {        if(TR[tx].dep<TR[ty].dep) swap(x,y),swap(tx,ty);        t[++num].a=ys[tx];t[num].b=ys[x];        x=TR[tx].fa;tx=TR[x].top;    }    if(TR[x].dep<TR[y].dep) swap(x,y);    t[++num].a=ys[y];t[num].b=ys[x];    sort(t+1,t+num+1,cmp);    int p=1;    for(int i=1;i<=num;i++)    {        if(p<t[i].a) change(root,1,n,p,t[i].a-1,k,op);        p=t[i].b+1;    }    if(p<=n) change(root,1,n,p,n,k,op);}int main(){    scanf("%d %d",&n,&m);    for(int i=1;i<n;i++)    {        int x,y;scanf("%d %d",&x,&y);        ins(x,y);ins(y,x);    }    pre_node(1,0);pre_edge(1,1);    for(int i=1;i<=m;i++)    {        int op;scanf("%d",&op);        if(op==0)        {            scanf("%d %d %d",&q[i].a,&q[i].b,&q[i].v);            solve(q[i].a,q[i].b,q[i].v,1);        }        if(op==1)        {            int t;scanf("%d",&t);            solve(q[t].a,q[t].b,q[t].v,2);        }        if(op==2)        {            int x;scanf("%d",&x);            printf("%d\n",findans(root,1,n,ys[x]));        }    }}