HDU 3308 LCIS 线段树区间更新

最近开始线段树一段时间了，也发现了不少大牛的博客比如HH大牛  ，小媛姐。这个题目是我在看HH大牛的线段树专题是给出的习题，（可以去他博客找找，真心推荐）原本例题是POJ3667 Hotel  这个题目，是一个求连续空区间的情况，而hdoj这个题目是求给定区间单调连续的最大区间长度，两个题目思路很相似，将节点rt用sum[rt],lsum[rt],rsum[rt]来描述，分别表示rt对应区间即[l,r]内满足条件的区间的最大长度，从左边端点l开始满足条件的最大区间长度，从右边r开始向左的满足条件的最大区间长度。

void PushUp(int rt,int m,int mid)//mid表示更新区间的中点，m表示长度
{
    lsum[rt]=lsum[rt<<1];
    rsum[rt]=rsum[rt<<1|1];
    int t=1;
    if(A[mid]<A[mid+1])
    {
        if(lsum[rt]==(m-(m>>1)))
            lsum[rt]+=lsum[rt<<1|1];
        if(rsum[rt]==(m>>1))
            rsum[rt]+=rsum[rt<<1];
        t= rsum[rt<<1]+lsum[rt<<1|1];
    }
       sum[rt]=max(t,max(sum[rt<<1],sum[rt<<1|1]));
}

PushUp函数和query函数是关键  ，PushUp函数向上更新时，先将lsum[rt]用

lsum[rt<<1]也就是左子区间左边开始的最大长度，如果这个区间长度刚好是l到mid之间的长度，说明区间已经穿过中点了，应该在加上右子区间lsum[rt<<1|1]这部分,类似的可以更新rsum[rt]。对于sum[rt]应该对应rsum[rt<<1]+lsum[rt<<1|1]，sum[rt<<1],sum[rt<<1|1]三种情况中的最大值，因为满足条件的最长单调区间可以是左半部分，右半部分，或者贯穿中点。

理解了这部分那么query函数可以按照类似的思路来写，但要注意的是对于贯穿中点的判断条件是A[mid]<A[mid+1]，这样才能将两部分加起来。

期间又一次犯二TLE了，原本是用C写的，max函数是用的宏定义，结果里面含有query函数，这样最终可能会两次调用query函数，不超时才怪，好在以前犯过一次，被基友发现了，改成algorithm里面的max就过了。哎，有个能debug基友就是好啊。hhha~

#include<cstdio>

#include<cstring> #include<algorithm> #define N 100010 #define lson l,m,rt<<1 #define rson m+1,r,rt<<1|1 using namespace std; int sum[N<<2],rsum[N<<2],lsum[N<<2],A[N]; void PushUp(int rt,int m,int mid) { lsum[rt]=lsum[rt<<1]; rsum[rt]=rsum[rt<<1|1]; int t=1; if(A[mid]<A[mid+1]) { if(lsum[rt]==(m-(m>>1))) lsum[rt]+=lsum[rt<<1|1]; if(rsum[rt]==(m>>1)) rsum[rt]+=rsum[rt<<1]; t= rsum[rt<<1]+lsum[rt<<1|1]; } sum[rt]=max(t,max(sum[rt<<1],sum[rt<<1|1])); } void build(int l,int r,int rt) { if(l==r) { sum[rt]=lsum[rt]=rsum[rt]=1; return; } int m=(l+r)>>1; build(lson); build(rson); PushUp(rt,r-l+1,m); } void update(int p,int l,int r,int rt) { if(l==r) return; int m=(l+r)>>1; if(p<=m) update(p,lson); else update(p,rson); PushUp(rt,r-l+1,m); } int query(int L,int R,int l,int r,int rt) { if(L<=l&&R>=r) { return sum[rt]; } int m=(l+r)>>1; if(l==r) return 1; int ans=0; if(R<=m) ans=max(ans,query(L,R,lson)); else if(L>m) ans=max(ans,query(L,R,rson)); else { ans=max(ans,query(L,R,lson)); ans=max(ans,query(L,R,rson)); int ll,rr; if(m-L+1>=rsum[rt<<1]) ll=rsum[rt<<1]; else ll=m-L+1; if(R-m>=lsum[rt<<1|1]) rr=lsum[rt<<1|1]; else rr=R-m; ans=max(ans,ll); ans=max(ans,rr); if(A[m+1]>A[m]) ans=max(ans,ll+rr); } return ans; } int main(void) { int T; int n,m,a,b; char op[3]; scanf("%d",&T); while(T--) { scanf("%d%d",&n,&m); for(int i=1; i<=n; i++) scanf("%d",A+i); build(1,n,1); while(m--) { scanf("%s%d%d",op,&a,&b); if(op[0]=='U') { A[a+1]=b; update(a+1,1,n,1); } else { int ans=query(a+1,b+1,1,n,1); printf("%d\n",ans); } } } return 0; }