BZOJ1756 小白逛公园 (线段树,区间最大子段和)

题意:

操作1:每次查询[l, r]内的最大子段和,操作2:修改a[x] 为 y。

分析:

复习了下最大子段和和最大子矩阵和,

先说最大子段和,可以进行DP,取dp[i] 表示到 i  为止的最大子段和,那么很容易得到状态转移方程dp[i]=max(dp[i-1]+a[i],a[i])o(n)的时间复杂度就可以解决。

再说最大子矩阵和,也可以进行DP,仿照上面最大子段和的思路,将矩阵转化为一维数组就好了,具体就是选择 i 行到 j 行的矩阵(列就是全部列),i-j的每一列求和然后压缩成一维数组,计算最大子段和就好了,所以时间复杂度o(m^2*n)

还有一种方法用线段树来做,时间复杂度\tiny o(n*m*logn) 根下面所讲的方法还有上面所讲的方法类似,从最后一行往上进行,每次加入一行,计算最大的子段和,类似的题目可以看HDU6638

再说区间最大子段和,因为还有修改操作所以很容易想到线段树,定义几个状态

sum 表示该区间的和。

lmx 表示该区间从最左边开始的最大子段和

rmx 表示该区间以最右边为结束的最大子段和

mx 表示该区间的最大子段和

那么很容易向上转移:

\tiny tree[i].sum=tree[i<<1].sum+tree[i<<1|1].sum;

\tiny tree[i].lmx=max(tree[i<<1].sum+tree[i<<1|1].lmx,tree[i<<1].lmx);

\tiny tree[i].rmx=max(tree[i<<1|1].sum+tree[i<<1].rmx,tree[i<<1|1].rmx);

\tiny tree[i].mx=max(max(tree[i<<1].mx,tree[i<<1|1].mx),tree[i<<1].rmx+tree[i<<1|1].lmx);

复杂度\tiny o(nlogn)

从而进行区间查询即可。

#include<bits/stdc++.h>

#define mm(a,b) memset(a,b,sizeof(a))
#define ACCELERATE (ios::sync_with_stdio(false),cin.tie(0))
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
#define PI acos(-1.0)
#define E exp(1.0)
//#define io
using namespace std;

const int inf=0x3f3f3f3f;

const int maxn=5e5+5;

struct node{
    int sum;
    int lmx,rmx;
    int mx;
}tree[maxn<<2];

int p[maxn];

void push_up(int i){
    tree[i].sum=tree[i<<1].sum+tree[i<<1|1].sum;
    tree[i].lmx=max(tree[i<<1].sum+tree[i<<1|1].lmx,tree[i<<1].lmx);
    tree[i].rmx=max(tree[i<<1|1].sum+tree[i<<1].rmx,tree[i<<1|1].rmx);
    tree[i].mx=max(max(tree[i<<1].mx,tree[i<<1|1].mx),tree[i<<1].rmx+tree[i<<1|1].lmx);
}

void build(int l,int r,int i){
    if(l==r){
        tree[i].lmx=tree[i].mx=tree[i].rmx=tree[i].sum=p[l];
        return ;
    }
    int mid=l+r>>1;
    build(l,mid,i<<1);
    build(mid+1,r,i<<1|1);
    push_up(i);
}

node query(int l,int r,int x,int y,int i){
    if(x<=l&&r<=y) return tree[i];
    int mid=l+r>>1;
    if(y<=mid){
        return query(l,mid,x,y,i<<1);
    }else if(x>mid){
        return query(mid+1,r,x,y,i<<1|1);
    }else{
        node t1=query(l,mid,x,mid,i<<1);
        node t2=query(mid+1,r,mid+1,y,i<<1|1);
        node ans;
        ans.lmx=max(t1.lmx,t1.sum+t2.lmx);
        ans.rmx=max(t2.rmx,t2.sum+t1.rmx);
        ans.sum=t1.sum+t2.sum;
        ans.mx=max(max(t1.mx,t2.mx),t1.rmx+t2.lmx);
        return ans;
    }
}

void change(int l,int r,int x,int y,int i){
    if(l==r){
        tree[i].lmx=tree[i].mx=tree[i].rmx=tree[i].sum=y;
        return ;
    }
    int mid=l+r>>1;
    if(mid>=x){
        change(l,mid,x,y,i<<1);
    }else{
        change(mid+1,r,x,y,i<<1|1);
    }
    push_up(i);
}

int main(){
    #ifdef io
    freopen("in.txt","r",stdin);
    freopen("out.txt","w",stdout);
    #endif

    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        scanf("%d",&p[i]);
    }
    build(1,n,1);
    while(m--){
        int op,x,y;
        scanf("%d%d%d",&op,&x,&y);
        if(op==1){
            if(x>y) swap(x,y);
            printf("%d\n",query(1,n,x,y,1).mx);
        }else{
            change(1,n,x,y,1);

        }
    }
    return 0;
}

 

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值