poj 3273 Tree(树链剖分)

来源：互联网发布：ps网络培训班编辑：程序博客网时间：2024/06/15 21:16

题意：给出一棵n个结点的树，有三个操作：把第i条边的权值改为v，把a到b路径上的边的权值取反，求a到b的路径上边的权值最大值。

思路：这题除了比较麻烦一点外，就是个裸的树链剖分。用线段树维护链的时候维护最小值，最大值，还有个取反标记就行了，取反的话最大值等于负的最小值，最小值等于负的最大值……由于一些2B错误改了半天才过。。。唉

代码：

#include<iostream>#include<cstdio>#include<cstring>#include<string>#include<algorithm>#include<map>#include<queue>#include<set>#include<stack>#include<cmath>#include<vector>#define inf 0x3f3f3f3f#define Inf 0x3FFFFFFFFFFFFFFFLL#define eps 1e-9#define pi acos(-1.0)using namespace std;typedef long long ll;const int maxn=10000+10;struct Edge{    int u,v,w,next;    Edge(){};    Edge(int uu,int vv,int ww,int nx){u=uu;v=vv;w=ww;next=nx;}}edges[maxn<<1];int head[maxn],d[maxn],nEdge,N;int pa[maxn],ft[maxn],childs[maxn],id[maxn<<1],val[maxn];bool flag[maxn<<1];void AddEdge(int u,int v,int w){    edges[++nEdge]=Edge(u,v,w,head[u]);    head[u]=nEdge;    edges[++nEdge]=Edge(v,u,w,head[v]);    head[v]=nEdge;}void Init(){    memset(head,0xff,sizeof(head));    memset(ft,0xff,sizeof(ft));    memset(flag,0,sizeof(flag));    memset(d,0,sizeof(d));    nEdge=-1;N=0;}int dfs(int u,int fa){    childs[u]=0;    int cnt=1,mx=-1,p;    for(int k=head[u];k!=-1;k=edges[k].next)    {        int v=edges[k].v;        if(v==fa) continue;        childs[u]++;d[v]=d[u]+1;        pa[v]=k;        int ct=dfs(v,u);        if(ct>mx) {mx=ct;p=k;}        cnt+=ct;    }    if(mx!=-1) flag[p]=flag[p^1]=true;    return cnt;}void dfs2(int u,int fa){    if(fa==-1) ft[u]=-1;    else    {        if(flag[fa])        {            id[fa]=id[fa^1]=++N;val[N]=edges[fa].w;            if(childs[edges[fa].u]<=1&&ft[edges[fa].u]!=-1) ft[u]=ft[edges[fa].u];            else ft[u]=fa;        }        else ft[u]=-1;    }    for(int k=head[u];k!=-1;k=edges[k].next)    {        if(k!=(fa^1))            dfs2(edges[k].v,k);    }}int maxv[maxn<<2],minv[maxn<<2],neg[maxn<<2];void PushUp(int rt){    maxv[rt]=max(maxv[rt<<1],maxv[rt<<1|1]);    minv[rt]=min(minv[rt<<1],minv[rt<<1|1]);}void PushDown(int rt){    if(neg[rt])    {        neg[rt<<1]^=1;        neg[rt<<1|1]^=1;        int tmp=minv[rt<<1];        minv[rt<<1]=-maxv[rt<<1];        maxv[rt<<1]=-tmp;        tmp=minv[rt<<1|1];        minv[rt<<1|1]=-maxv[rt<<1|1];        maxv[rt<<1|1]=-tmp;        neg[rt]=0;    }}void build(int l,int r,int rt){    neg[rt]=0;    if(l==r)    {        minv[rt]=maxv[rt]=val[l];        return ;    }    int m=(l+r)>>1;    build(l,m,rt<<1);    build(m+1,r,rt<<1|1);    PushUp(rt);}int Query(int L,int R,int l,int r,int rt){    if(l>=L&&r<=R)    {        return maxv[rt];    }    PushDown(rt);    int m=(l+r)>>1;    if(m>=R) return Query(L,R,l,m,rt<<1);    else if(m<L) return Query(L,R,m+1,r,rt<<1|1);    return max(Query(L,R,l,m,rt<<1),Query(L,R,m+1,r,rt<<1|1));}void Update(int p,int l,int r,int rt,int v){    if(l==r)    {        neg[rt]=0;        minv[rt]=maxv[rt]=v;        return ;    }    PushDown(rt);    int m=(l+r)>>1;    if(m>=p) Update(p,l,m,rt<<1,v);    else Update(p,m+1,r,rt<<1|1,v);    PushUp(rt);}void Negate(int L,int R,int l,int r,int rt){    if(l>=L&&r<=R)    {        neg[rt]^=1;        int tmp=minv[rt];        minv[rt]=-maxv[rt];        maxv[rt]=-tmp;        return;    }    PushDown(rt);    int m=(l+r)>>1;    if(m>=L) Negate(L,R,l,m,rt<<1);    if(m<R) Negate(L,R,m+1,r,rt<<1|1);    PushUp(rt);}int findmax(int x,int y){    if(x==y) return 0;    int maxres=-inf,L,R;    while(x!=y)    {        if(d[x]<d[y]) swap(x,y);        if(ft[x]!=-1)        {            L=id[ft[x]];R=id[pa[x]];            if(ft[x]==ft[y])            {                L=id[pa[y]]+1;                x=y;            }            else x=edges[ft[x]].u;            maxres=max(maxres,Query(L,R,1,N,1));        }        else        {            maxres=max(maxres,edges[pa[x]].w);            x=edges[pa[x]].u;        }    }    return maxres;}void negseg(int x,int y){    if(x==y) return ;    int L,R;    while(x!=y)    {        if(d[x]<d[y]) swap(x,y);        if(ft[x]!=-1)        {            L=id[ft[x]];R=id[pa[x]];            if(ft[x]==ft[y])            {                L=id[pa[y]]+1;                x=y;            }            else x=edges[ft[x]].u;            Negate(L,R,1,N,1);        }        else        {            edges[pa[x]].w*=-1;            x=edges[pa[x]].u;        }    }}void changes(int x,int v){    if(flag[x*2]) Update(id[x*2],1,N,1,v);    else edges[x*2].w=edges[(x*2)^1].w=v;}int main(){    //freopen("in.txt","r",stdin);    //freopen("out.txt","w",stdout);    int t,n;    scanf("%d",&t);    char str[10];    while(t--)    {        Init();        scanf("%d",&n);        int u,v,w;        for(int i=1;i<n;++i)        {            scanf("%d%d%d",&u,&v,&w);            AddEdge(u,v,w);        }        dfs(1,-1);dfs2(1,-1);        if(N) build(1,N,1);        while(true)        {            scanf("%s",str);            if(str[0]=='D') break;            scanf("%d%d",&u,&v);            if(str[0]=='Q')            {                w=findmax(u,v);                printf("%d\n",w);            }            else if(str[0]=='N')                negseg(u,v);            else changes(u-1,v);        }    }    return 0;}

0 0