旅游（国家集训队）

最新推荐文章于 2024-08-28 15:11:10 发布

Hunter_HunterL

最新推荐文章于 2024-08-28 15:11:10 发布

阅读量306

点赞数 1

CC 4.0 BY-SA版权

分类专栏：树剖文章标签：树剖

本文链接：https://blog.youkuaiyun.com/Hunter_HunterL/article/details/77206470

树剖专栏收录该内容

1 篇文章

订阅专栏

本文介绍了一种利用树链剖分和线段树处理路径操作的方法，包括节点值更新、区间求和、最大值和最小值查询等功能，并通过具体的代码实现了这些操作。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题面

树链剖分，边权存到点上，注意边界处理

#include<iostream>
#include<cstdio>
#include<cstring>
#define maxn 100005
using namespace std;
int n,m;
struct edge
{
    int to,ne,w;  
}b[maxn];
int k=0,head[maxn];
struct tree
{
    int l,r,mx,mi,sum; 
}t[maxn*4];
struct edge2
{
     int s,t,w; 
}a[maxn];
int c[maxn];
bool rev[maxn];
int max(int x,int y)
{
    return x >y ? x: y;
}
void swap(int &x,int &y)
{
    int z=x;x=y;y=z; 
}
inline int read()
{
     int x=0;char ch=getchar();
     while(ch<'0'||ch>'9') ch=getchar();
     while(ch>='0'&&ch<='9'){ x=x*10+ch-'0'; ch=getchar(); }
     return x;  
}
void add(int u,int v,int w)
{
    k++;
    b[k].to=v;b[k].ne=head[u];b[k].w=w;head[u]=k;
}
int fa[maxn],sz[maxn],son[maxn],d[maxn],top[maxn],id[maxn],pos[maxn];
void dfs1(int x)
{
    sz[x]=1; 
    for(int i=head[x];i!=-1;i=b[i].ne)
    if(b[i].to!=fa[x]){
       fa[b[i].to]=x; d[b[i].to]=d[x]+1;
       c[b[i].to]=b[i].w;
       dfs1(b[i].to);
       sz[x]+=sz[b[i].to];
       if(son[x]==-1) son[x]=b[i].to;
       else if(sz[b[i].to]>sz[son[x]]) son[x]=b[i].to;
    }
}
int cnt=0;
void dfs2(int x,int tp)
{
     id[x]=++cnt;  pos[cnt]=x;
     top[x]=tp;
     if(son[x]!=-1) dfs2(son[x],tp);
     for(int i=head[x];i!=-1;i=b[i].ne)
     if(b[i].to!=fa[x]&&b[i].to!=son[x])  dfs2(b[i].to,b[i].to);
}
struct  Tree
{
        void swp(int &x,int &y)
        {
             int z=x; x=-y; y=-z;
        }
        void update(int z)
        {
            t[z].mx=max(t[z<<1].mx,t[z<<1|1].mx);
            t[z].mi=min(t[z<<1].mi,t[z<<1|1].mi);
            t[z].sum=t[z<<1].sum+t[z<<1|1].sum; 
        }
        void build(int l,int r,int z)
        {
            t[z].l=l;t[z].r=r;
            if(l==r){
               t[z].mx=t[z].mi=t[z].sum=c[pos[l]];
               return ;
            }
            int mid=(l+r)>>1;
            build(l,mid,z<<1);
            build(mid+1,r,z<<1|1);
            update(z);
        }
        void xg(int x,int y,int z)
        {
            if(t[z].l==t[z].r){
               t[z].mx=t[z].mi=t[z].sum=y;
               return ;
            }
            int mid=(t[z].l+t[z].r)>>1;
            if(x<=mid) xg(x,y,z<<1);
            else xg(x,y,z<<1|1);
            update(z);
        }
        void Rev(int x,int y,int z)
        {
             if(t[z].l==t[z].r){
                t[z].mx*=-1;t[z].mi*=-1;t[z].sum*=-1;
                return ;
             }
             int mid=(t[z].l+t[z].r)>>1;
             if(x<=mid) Rev(x,y,z<<1);
             if(y>mid) Rev(x,y,z<<1|1);
             update(z);
        }
        int getsum(int x,int y,int z)
        {
             if(x<=t[z].l&&t[z].r<=y) return t[z].sum;
             int ans=0;
             int mid=(t[z].l+t[z].r)>>1;
             if(x<=mid) ans+=getsum(x,y,z<<1);
             if(y>mid) ans+=getsum(x,y,z<<1|1);
             return ans;
        }
        int getmax(int x,int y,int z)
        {
             if(x<=t[z].l&&t[z].r<=y) return t[z].mx;
             int ans=-0x7fffffff;
             int mid=(t[z].l+t[z].r)>>1;
             if(x<=mid) ans=max(ans,getmax(x,y,z<<1));
             if(y>mid) ans=max(ans,getmax(x,y,z<<1|1));
             return ans;
        }   
        int getmin(int x,int y,int z)
        {
             if(x<=t[z].l&&t[z].r<=y) return t[z].mi;
             int ans=0x7fffffff;
             int mid=(t[z].l+t[z].r)>>1;
             if(x<=mid) ans=min(ans,getmin(x,y,z<<1));
             if(y>mid) ans=min(ans,getmin(x,y,z<<1|1));
             return ans;
        }   
}tr;
void revc(int x,int y)
{
     int fx=top[x],fy=top[y];
     while(fx!=fy){
        if(d[fx]<d[fy])  { swap(fx,fy); swap(x,y); }
        tr.Rev(id[fx],id[x],1);
        x=fa[fx]; fx=top[x];
     }
     if(d[x]>d[y]) swap(x,y);
     tr.Rev(id[x]+1,id[y],1);
}
int getsum(int x,int y)
{
     int fx=top[x],fy=top[y];
     int ans=0;
     while(fx!=fy){
         if(d[fx]<d[fy]) { swap(fx,fy); swap(x,y); }
         ans+=tr.getsum(id[fx],id[x],1);
         x=fa[fx]; fx=top[x];
     }
     if(d[x]>d[y]) swap(x,y);
     ans+=tr.getsum(id[x]+1,id[y],1);
     return ans;
}
int getmax(int x,int y)
{
     int fx=top[x],fy=top[y];
     int ans=-0x7fffffff;
     while(fx!=fy){
         if(d[fx]<d[fy]) { swap(fx,fy); swap(x,y); }
         ans=max(ans,tr.getmax(id[fx],id[x],1));
         x=fa[fx]; fx=top[x];
     }
     if(d[x]>d[y]) swap(x,y);
     ans=max(ans,tr.getmax(id[x]+1,id[y],1));
     return ans;
}
int getmin(int x,int y)
{
     int fx=top[x],fy=top[y];
     int ans=0x7fffffff;
     while(fx!=fy){
         if(d[fx]<d[fy]) { swap(fx,fy); swap(x,y); }
         ans=min(ans,tr.getmin(id[fx],id[x],1));
         x=fa[fx]; fx=top[x];
     }
     if(d[x]>d[y]) swap(x,y);
     ans=min(ans,tr.getmin(id[x]+1,id[y],1));
     return ans;
}
int main()
{
    //freopen("in.txt","r",stdin);
    freopen("nt2011_travel.in","r",stdin);
    freopen("nt2011_travel.out","w",stdout);
    memset(head,-1,sizeof(head));
    memset(son,-1,sizeof(son));
    memset(fa,-1,sizeof(fa));
    scanf("%d",&n);
    int x,y,z;
    char type[5];
    for(int i=1;i<n;i++){
       scanf("%d%d%d",&x,&y,&z);
       add(x,y,z);  add(y,x,z);
       a[i].s=x;a[i].t=y;
    }
    dfs1(0);
    dfs2(0,0);
    tr.build(1,n+1,1);
    scanf("%d",&m);
    while(m--){
        scanf("%s",type); scanf("%d%d",&x,&y);
        if(type[0]=='C'){
           if(fa[a[x].s]==a[x].t) tr.xg(id[a[x].s],y,1);
           else  tr.xg(id[a[x].t],y,1);
        }
        if(type[0]=='N')  revc(x,y); 
        if(type[0]=='S')  printf("%d\n",getsum(x,y)); 
        if(type[0]=='M'&&type[1]=='A') printf("%d\n",getmax(x,y));
        if(type[0]=='M'&&type[1]=='I') printf("%d\n",getmin(x,y));
    }
    return 0;
}