SPOJGSS7 Can you answer these queries VII

最新推荐文章于 2024-06-21 17:54:11 发布

原创最新推荐文章于 2024-06-21 17:54:11 发布 · 305 阅读

0 ·

CC 4.0 BY-SA版权

动态树同时被 2 个专栏收录

12 篇文章

订阅专栏

spoj

4 篇文章

订阅专栏

本文介绍了一种解决树状结构中寻找最大连续子段和问题的方法。通过使用线段树和动态规划技巧，文章详细阐述了如何高效地合并节点信息以找到最优解，并提供了完整的代码实现。

题意

树上最大连续子段和

题解

维护前后缀和总的最大连续子段和
合并即可

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#define L(i) (T[i].s[0])
#define R(i) (T[i].s[1])
#define F(i) (T[i].fa)
#define S(i) (T[i].sum)
#define sz(i) (T[i].siz)
#define Loc(i) (R(F(i))==i)
#define Ls(i) (T[i].ls)
#define Rs(i) (T[i].rs)
#define ms(i) (T[i].ms)
#define For(i,j,k) for(int i=(j);i<=(int)k;i++)
#define Forr(i,j,k) for(int i=(j);i>=(int)k;i--)
#define Set(a,b) memset(a,b,sizeof(a))
#define ll long long 
using namespace std;
const int N=100010,INF=0x3f3f3f3f;
int to[N*2],Next[N*2],Begin[N],e;
inline void add(int x,int y){
    to[++e]=y,Next[e]=Begin[x],Begin[x]=e;
}
inline void read(int &x){
    x=0;char c=getchar();int f(0);
    while(c<'0'||c>'9')f|=(c=='-'),c=getchar();
    while(c>='0'&&c<='9')x=x*10+c-'0',c=getchar();
    if(f)x=-x;
}
struct node{
    int s[2],fa,rev,val,mk,siz;
    ll sum,ls,rs,ms;
    node(){s[0]=s[1]=fa=rev=sum=val=ls=rs=ms=0,mk=-INF,siz=1;}
};
struct LCT{
    node T[N];
    LCT(){T[0].siz=0;}
    inline void pushup(int x){
        sz(x)=sz(L(x))+sz(R(x))+1;
        S(x)=S(L(x))+S(R(x))+T[x].val;
        Ls(x)=max(Ls(L(x)),S(L(x))+T[x].val+Ls(R(x)));
        Rs(x)=max(Rs(R(x)),S(R(x))+T[x].val+Rs(L(x)));
        ms(x)=max(max(ms(L(x)),ms(R(x))),Rs(L(x))+T[x].val+Ls(R(x)));
    }
    inline void Swap(int x){
        swap(Ls(x),Rs(x)),swap(L(x),R(x)),T[x].rev^=1;
    }
    inline void mark(int x,int v){
        T[x].mk=v,T[x].val=v,S(x)=v*sz(x);
        Ls(x)=Rs(x)=ms(x)=max(0ll,S(x));
    }
    inline void pushdown(int x){
        if(T[x].rev){
            T[x].rev=0;
            Swap(L(x)),Swap(R(x));
        }
        if(T[x].mk>-INF){
            if(L(x))mark(L(x),T[x].mk);
            if(R(x))mark(R(x),T[x].mk);
            T[x].mk=-INF;
        }
    }
    inline bool isrt(int x){
        return R(F(x))!=x&&L(F(x))!=x;
    }
    inline void Pushdown(int x){
        if(!isrt(x))Pushdown(F(x));
        pushdown(x);
    }
    inline void Rotate(int x){
        int A=F(x),B=F(A),l=Loc(x),r=l^1,d=Loc(A);
        if(!isrt(A))T[B].s[d]=x;F(x)=B;
        F(A)=x,F(T[x].s[r])=A,T[A].s[l]=T[x].s[r],T[x].s[r]=A;
        pushup(A);
    }   
    inline void splay(int x){
        Pushdown(x);
        while(!isrt(x)){
            if(!isrt(F(x)))Rotate(x);
            Rotate(x);
        }
                pushup(x);
    }
    inline void access(int x){
        for(int i=0;x;i=x,x=F(x))
            splay(x),R(x)=i,pushup(x);
    }
    inline void reverse(int x){
        access(x),splay(x),Swap(x);
    }
    inline void split(int x,int y){
        reverse(x),access(y),splay(y);
    }
    inline void modify(int x,int y){
        int v;read(v);
        split(x,y);mark(y,v);
    }
    inline void query(int x,int y){
        split(x,y);
        printf("%lld\n",ms(y));
    }
}t;
void dfs(int u,int fa){
    for(int i=Begin[u];i;i=Next[i]){
        int v=to[i];
        if(v==fa)continue;
        t.T[v].fa=u;
        dfs(v,u);
    }
}
int main(){
    int n,m,u,v,p;
    read(n);
    For(i,1,n){
        read(v);
        t.mark(i,v),t.T[i].mk=0;
    }
    For(i,1,n-1){
        read(u),read(v);
        add(u,v),add(v,u);
    }
    dfs(1,0);
    read(m);
    while(m--){
        read(p),read(u),read(v);
        if(p==1)t.query(u,v);
        else t.modify(u,v);
    }
    return 0;
}