树剖

树剖

只是想来放一波代码= =

单点修改+区间求和+区间最大值

区间修改就在线段树上加个Lazy操作就可以了

注意:代码里面,权值都是在点上的

#include<iostream> 
#include<cstdio> 
#define N 300000 
#define ll long long 
using namespace std; 
ll maxxx=0; 
ll V[N]; 
ll Last[N],Next[N],End[N],tot; 
ll Size[N],Fa[N],Dep[N],Son[N],Top[N],Dfn[N],Old[N],cnt; 
//all of the arrays use old ids 
struct node{ll l,r,v,Max;}tree[N]; 
void Find(ll u,ll uu,ll dep){ 
    Size[u]=1,Fa[u]=uu,Dep[u]=dep;
    for(ll i=Last[u];i;i=Next[i]){  
        ll v=End[i]; 
        if(v==Fa[u])continue;//Go Wrong 1
        Find(v,u,dep+1);  
        Size[u]+=Size[v]; 
        if(Size[v]>Size[Son[u]])Son[u]=v;
    } 
} 
void Connect(ll u,ll anc){ 
    Top[u]=anc,Dfn[u]=++cnt,Old[cnt]=u; 
    if(Son[u])Connect(Son[u],anc); 
    for(ll i=Last[u];i;i=Next[i]){ 
        ll v=End[i]; 
        if(v!=Son[u] && v!=Fa[u])Connect(v,v); 
    } 
} 
void Make(ll p,ll x,ll y){ 
    tree[p].l=x,tree[p].r=y; 
    if(x>=y){ 
        tree[p].v=V[Old[x]];//Go Wrong 2:tree[p].v=V[x];
        tree[p].Max=V[Old[x]];//Go Wrong 2
        return; 
    } 
    ll mid=(x+y)>>1; 
    Make(p<<1,x,mid),Make((p<<1)+1,mid+1,y); 
    tree[p].v=tree[p<<1].v+tree[(p<<1)+1].v; 
    tree[p].Max=max(tree[p<<1].Max,tree[(p<<1)+1].Max); 
} 
void Modify(ll p,ll u,ll t){ 
    if(tree[p].l==tree[p].r){ 
        tree[p].v=t,tree[p].Max=t; 
        return; 
    }
    int mid=(tree[p].l+tree[p].r)>>1;
    if(u<=mid)Modify(p<<1,u,t);
    else Modify((p<<1)+1,u,t); 
    tree[p].v=tree[p<<1].v+tree[(p<<1)+1].v; 
    tree[p].Max=max(tree[p<<1].Max,tree[(p<<1)+1].Max); 
} 
ll TreeSum(ll p,ll x,ll y){ 
    if(tree[p].r<x || tree[p].l>y)return 0; 
    if(tree[p].l>=x && tree[p].r<=y)return tree[p].v; 
    return TreeSum(p<<1,x,y)+TreeSum((p<<1)+1,x,y); 
} 
ll GetSum(ll x,ll y){ 
    ll Ans=0; 
    while(Top[x]!=Top[y]){ 
        if(Dep[Top[x]]<Dep[Top[y]])swap(x,y); 
        Ans+=TreeSum(1,Dfn[Top[x]],Dfn[x]); 
        x=Fa[Top[x]];//Go Wrong 3
    } 
    if(Dep[x]>Dep[y])swap(x,y); 
    Ans+=TreeSum(1,Dfn[x],Dfn[y]); 
    return Ans; 
} 
ll TreeMax(ll p,ll x,ll y){ 
    if(tree[p].l>y || tree[p].r<x)return -1e9; 
    if(tree[p].l>=x && tree[p].r<=y)return tree[p].Max; 
    return max(TreeMax(p<<1,x,y),TreeMax((p<<1)+1,x,y)); 
} 
ll GetMax(ll x,ll y){ 
    ll Ans=-1e9; 
    while(Top[x]!=Top[y]){ 
        if(Dep[Top[x]]<Dep[Top[y]])swap(x,y); 
        Ans=max(Ans,TreeMax(1,Dfn[Top[x]],Dfn[x])); 
        x=Fa[Top[x]];//Go Wrong 3
    } 
    if(Dep[x]>Dep[y])swap(x,y); 
    Ans=max(Ans,TreeMax(1,Dfn[x],Dfn[y])); 
    return Ans; 
} 
void Ins(ll x,ll y){ 
    End[++tot]=y; 
    Next[tot]=Last[x]; 
    Last[x]=tot; 
} 
int main(){ 
    ll n,q;scanf("%lld",&n); 
    for(ll i=1;i<n;i++){ 
        ll a,b;scanf("%lld%lld",&a,&b); 
        Ins(a,b),Ins(b,a); 
    } 
    for(ll i=1;i<=n;i++)scanf("%lld",&V[i]); 
    Find(1,0,1); 
    Connect(1,1);
    Make(1,1,n);
    scanf("%lld",&q); 
    for(ll i=1;i<=q;i++){ 
        char s[10];ll a,b; 
        scanf("%s%lld%lld",s,&a,&b); 
        if(s[1]=='H')Modify(1,Dfn[a],b);
        if(s[1]=='M')printf("%lld\n",GetMax(a,b)); 
        if(s[1]=='S')printf("%lld\n",GetSum(a,b));
    }
    return 0; 
}  

批注了Go Wrong的地方都是WA了之后才发现的错误之处。希望以后不会再在这些地方出错。

Something Goes Wrong !

Go Wrong 1

没有跳过v==father[u]的时候的情况。由于连的是无向边,所以在建立搜索树的时候必须判断。

Go Wrong 2

线段树里面每个点的编号都是搜索树里的新编号,不是题目里直接给的编号。在赋值的时候应该对应新编号。

Go Wrong 3

当操作完该条重链之后,x应该跳到这条链顶端的点的父亲节点,而不是跳到自己的父亲节点(这里指的都是搜索树里面的父子关系)。x=father[x]是标准的错误写法,保证又WA又TLE,一定是x=father[top[x]]

LCA

void GetLCA(int x,int y){
    while(Top[x]!=Top[y]){
        if(Dep[Top[x]]<Dep[Top[y]])swap(x,y);
        x=Fa[Top[x]];
    }
    if(Dep[x]>Dep[y])swap(x,y);
    return x;
}
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
#include<bits/stdc++.h> using namespace std; const int mx=1e5+1; int n,Q,x,y,d[mx],fa[mx],siz[mx],ev[mx],a[mx],son[mx],dfn[mx],cnt,id[mx],top[mx],ans[mx]; struct edge{int c,w,id,u,v;}e[mx*2]; struct que{int u,v,x,y;}q[mx*2]; struct tree{int l,r,lzy1,lzy2;}t[mx*4]; vector<edge> v[mx]; vector<int> es[mx]; vector<int> qs[mx]; //以下树剖 void dfs1(int f,int u) { d[u]=d[f]+1,fa[u]=f,siz[u]=1; int len=v[u].size(); for(int i=0;i<len;i++) { edge next=v[u][i]; int nv=next.v; if(nv==f) continue; ev[next.id]=nv,a[nv]=next.w; dfs1(u,nv); siz[u]+=siz[nv]; if(siz[nv]>siz[son[u]]) son[u]=nv; } } void dfs2(int f,int u) { dfn[u]=++cnt,id[cnt]=u,top[u]=f; if(son[u]) dfs2(f,son[u]); int len=v[u].size(); for(int i=0;i<len;i++) { int nv=v[u][i].v; if(nv==fa[u] || nv==son[u]) continue; dfs2(nv,nv); } } //以上树剖 //以下线段树 void pushup1(int x){t[x].lzy1=t[x<<1].lzy1+t[x<<1|1].lzy1;} void pushup2(int x){t[x].lzy2=t[x<<1].lzy2+t[x<<1|1].lzy2;} void build(int x,int l,int r) { t[x].l=l,t[x].r=r; if(l==r) { t[x].lzy1=a[id[l]],t[x].lzy2=0; return; } int mid=(l+r)/2; build(x<<1,l,mid);build(x<<1|1,mid+1,r); pushup1(x); } void chang1(int x,int obx,int w) { if(t[x].l==t[x].r){t[x].lzy1=w;return;} int mid=(t[x].l+t[x].r)>>1; if(obx<=mid) chang1(x<<1,obx,w); else chang1(x<<1|1,obx,w); pushup1(x); } void chang2(int x,int obx,int w) { if(t[x].l==t[x].r){t[x].lzy2=w;return;} int mid=(t[x].l+t[x].r)>>1; if(obx<=mid) chang2(x<<1,obx,w); else chang2(x<<1|1,obx,w); pushup2(x); } int find1(int x,int l,int r) { if(l<=t[x].l && r>=t[x].r) return t[x].lzy1; int mid=(l+r)>>1,s=0; if(l<=mid) s+=find1(x<<1,l,r); if(r>mid) s+=find1(x<<1|1,l,r); return s; } int find2(int x,int l,int r) { if(l<=t[x].l && r>=t[x].r) return t[x].lzy2; int mid=(l+r)>>1,s=0; if(l<=mid) s+=find2(x<<1,l,r); if(r>mid) s+=find2(x<<1|1,l,r); return s; } //以上线段树 int fans(int x,int y,int k) { int ans=0; while(top[x]!=top[y]) { if(d[top[x]]<d[top[y]]) swap(x,y); ans+=find1(1,dfn[top[x]],dfn[x]); ans+=find2(1,dfn[top[x]],dfn[x]); x=fa[top[x]]; } if(d[x]>d[y]) swap(x,y); if(x!=y) { ans+=find1(1,dfn[x]+1,dfn[y]); ans+=k*find2(1,dfn[x]+1,dfn[y]); } return ans; } int main() { cin >> n >> Q; for(int i=1;i<n;i++) { cin >> e[i].u >> e[i].v >> e[i].c >> e[i].w; e[i].id=i; v[e[i].u].push_back({e[i].u,e[i].v,e[i].c,e[i].w,e[i].id}); v[e[i].v].push_back({e[i].v,e[i].u,e[i].c,e[i].w,e[i].id}); es[e[i].c].push_back(i); } for(int i=1;i<=Q;i++) { cin >> q[i].x >> q[i].y >> q[i].u >> q[i].v; qs[q[i].x].push_back(i); } dfs1(1,1);dfs2(1,1);build(1,1,n); for(int i=1;i<n;i++) { int len=es[i].size(); for(int j=0;j<len;j++) { int k=ev[es[i][j]]; find1(1,dfn[k],0); find2(1,dfn[k],1); } for(int j=0;j<len;j++) { int k=qs[i][j]; ans[k]=fans(q[k].u,q[k].v,q[k].y); } for(int j=0;j<len;j++) { int k=ev[es[i][j]]; find1(1,dfn[k],e[es[i][j]].w); find2(1,dfn[k],0); } } for(int i=1;i<=Q;i++) cout<<ans[i]<<"\n"; return 0; }
最新发布
07-14

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值