uestc Another LCIS(区间更新,区间合并)

题意:

求最长上升子序列的最大长度。

两种操作:

q: 某一段区间的最长上升子序列的最大长度为多少

a: 某一段区间都加上val

解题思路:

如果单单求最长上升子序列,我们可以通过dp,或者分治的思想来求。因为此处操作多,所以上述方法不可行。

我们在线段树的结点里增加几个域:

lsum :  表示从当前结点的最左端开始的上升子序列的长度

rsum:  表示从当前结点的最右端开始的上升子序列的长度

a :       表示整体的增幅

lazy:   lazy标记,lazy>0 表示还未向下更新

la:    表示当前结点的最左端开始的上升子序列的值的增幅

ra:    表示从当前结点的最右端开始的上升子序列值的增幅

sum:表示当前结点上升子序列长度的最大值

所以我们

向下更新的时候:

只要讲lazy 向下传递 增幅向下传递即可

向上更新的时候:

我们比价值的大小时要加上增幅(详见代码)


注意:

左边的上升子序列在比较时必须加上左边的增幅,右边的也一样

题目地址:    http://acm.uestc.edu.cn/#/problem/show/360

#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#define MAXN0 100010
#define l1(x) (x)<<1
#define r1(x) (x)>>1
#define max(a,b) ((a)>(b)?(a):(b))
struct TRnode{
    int L,R;
    int la,ra,a,lazy;
    int lsum,rsum,sum;
};
TRnode TR[MAXN0<<2];
int S[MAXN0];
void pu(int k,int k1,int k11){
    int tmp0;
    TR[k].la = TR[k1].la;
    TR[k].ra = TR[k11].ra;
    TR[k].lsum = TR[k1].lsum;
    tmp0 = TR[k1].R - TR[k1].L + 1 ;
    if(TR[k].lsum==tmp0){
        if(S[TR[k1].R]+TR[k1].ra<S[TR[k11].L]+TR[k11].la){
            TR[k].lsum += TR[k11].lsum;
        }
    }
    TR[k].rsum = TR[k11].rsum;
    tmp0 = TR[k11].R - TR[k11].L + 1;
    if(tmp0==TR[k].rsum){
        if(S[TR[k1].R]+TR[k1].ra<S[TR[k11].L]+TR[k11].la){
            TR[k].rsum += TR[k1].rsum;
        }
    }
    TR[k].sum = max(TR[k].lsum,TR[k].rsum);
    TR[k].sum = max(TR[k].sum,max(TR[k1].sum,TR[k11].sum));
    if(S[TR[k1].R]+TR[k1].ra<S[TR[k11].L]+TR[k11].la){
        //TR[k].lsum += TR[k11].lsum;
        TR[k].sum = max(TR[k].sum,TR[k1].rsum+TR[k11].lsum);
    }
}
void pd(int k,int k1,int k11){
    //TR[k].lazy = 0;
    TR[k1].la += TR[k].a;
    TR[k1].ra += TR[k].a;
    TR[k1].a +=TR[k].a;
    
    TR[k11].la += TR[k].a;
    TR[k11].ra += TR[k].a;
    TR[k11].a +=TR[k].a;
    
    TR[k1].lazy = TR[k11].lazy = TR[k].lazy;
    TR[k].a = 0;
    TR[k].lazy = 0;
}
void buildTR(int L,int R,int k){
    //int L,R;
    TR[k].L = L;
    TR[k].R = R;
    TR[k].lazy = TR[k].a = 0;
    TR[k].la = TR[k].ra = 0;
    if(L==R){
        TR[k].lsum = TR[k].rsum = TR[k].sum = 1;
        return;
    }
    int k1 = l1(k);
    int k11 = k1 + 1;
    int mid = r1(L+R);
    buildTR(L,mid,k1);
    buildTR(mid+1,R,k11);
    pu(k,k1,k11);
}
void update(int L,int R,int k,int a){
    if(TR[k].L==L&&TR[k].R==R){
        TR[k].la += a;
        TR[k].ra += a;
        TR[k].a  += a;
        TR[k].lazy = 1;
        return;
    }
    int k1,k11,mid;
    k1 = l1(k);
    k11 = k1 + 1;
    if(TR[k].lazy){
        pd(k,k1,k11);
    }
    mid = r1(TR[k].L+TR[k].R);
    if(mid>=R){
        update(L, R, k1, a);
    }
    else if(mid<L){
        update(L, R, k11, a);
    }
    else {
        update(L,mid,k1,a);
        update(mid+1,R,k11,a);
    }
    pu(k,k1,k11);
}
int query(int L,int R,int k){
    if(TR[k].L==L&&TR[k].R==R){
        return TR[k].sum;
    }
    int k1,k11;
    k1 = l1(k);
    k11 = k1+1;
    if(TR[k].lazy){
        pd(k,k1,k11);
    }
    int mid = r1(TR[k].L+TR[k].R);
    if(mid>=R){
        return query(L, R, k1);
    }
    else if(mid<L){
        return  query(L, R, k11);
    }
    else {
        int tmp0,tmp1,tmp2,tmp3;
        tmp0 = query(L, mid, k1);
        tmp1 = query(mid+1, R, k11);
        tmp3 = 0,tmp2 = 0;
        if(S[TR[k1].R]+TR[k1].ra<S[TR[k11].L]+TR[k11].la){
            tmp2 = TR[k1].rsum;
            tmp3 = mid - L + 1;
            if(TR[k1].rsum>tmp3){
                tmp2 = tmp3;
            }
            tmp3 = R - mid;
            if(TR[k11].lsum<tmp3){
                tmp3 = TR[k11].lsum;
            }
            tmp2+=tmp3;
        }
        return max(tmp0,max(tmp1,tmp2));
    }
}
int main(){
    int T,L,R;
    int N,Q,a,cas;
    
    char w[2];
    scanf("%d",&T);
    cas = 0;
    while(T--){
        scanf("%d%d",&N,&Q);
        for(int i=1;i<=N;++i){
            scanf("%d",&S[i]);
        }
        buildTR(1,N,1);
        printf("Case #%d:\n",++cas);
        while(Q--){
            scanf("%s%d%d",w,&L,&R);
            if(w[0]=='q'){
                int ans = query(L, R, 1);
                printf("%d\n",ans);
            }
            else {
                scanf("%d",&a);
                update(L,R,1,a);
            }
        }
    }
    return 0;
}






评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值