Contest Hunter 4301 Can you answer on these queries III (单点修改+查询区间最大子段和)

根据lyd 的模板,我们只需要在每次更新信息的时候 同时维护

t[p].sum=t[p*2].sum+t[p*2+1].sum;
t[p].lx=max(t[p*2].lx, t[p*2].sum+t[p*2+1].lx);

t[p].rx=max(t[p*2+1].rx, t[p*2].rx+t[p*2+1].sum);
t[p].dat=Max(t[p*2].dat, t[p*2+1].dat, t[p*2].rx+t[p*2+1].lx);

p节点所在区间的和,p节点所在区间的紧靠左边的最大字段和, 紧靠右边的最大字段和,区间内的最大字段和

每次change单点更新节点的时候 从下往上更新这些值

buil 建树的时候也需要从下往上更新这些值

问题来了区间查询ask (lyd没讲清楚,wa了一下午,最后还是去查了题解,巨菜ε=(´ο`*)))唉)

对于区间l~r 我们如果要查询l~r 的最大连续字段和, 那么肯定是分三种情况,l~r在左子区间,l~r在右子区间

l~r跨越两个区间,前两种情况我们只要递归查询左右子区间就行了。 

第三种情况,这时候就要用到我们之前一直维护的那些信息了。

首先我们递归找l~mid 左子区间 让他返回上述那些 信息,然后找mid~r 右子区间 同时也让他返回上面那些信息,

要同时返回这么多信息,肯定是不能直接返回某个值,于是我们用返回结构体(a和b)来保存这些信息,然后

就根据上面的递推式得到

        c.sum=a.sum+b.sum;
        c.lx=max(a.lx, a.sum+b.lx);
        c.rx=max(b.rx, a.rx+b.sum);
        c.dat=Max(a.dat, b.dat, a.rx+b.lx);

返回这个c 就是我们想要得到的

刚开始思路 以为可以直接得到最大连续字段和,而第三种情况必须根据上面的信息推导出来,所以每次递归都需要

返回这些信息,最后根据这些信息得出结果。

#include<cstdio>
#include<iostream>
#include<cmath>
#include<queue>
#include<cstring>
#define ll long long
using namespace std;
const int SIZE=500005;
ll a[SIZE];
int n,m;
ll Max(ll a, ll b, ll c)
{
    return max(max(a,b) , c);
}
struct segment_tree
{
    int l,r;
    ll sum;
    ll lx;
    ll rx;
    ll dat;
}t[SIZE*4];
void build(int p, int l, int r)
{
    t[p].l=l;
    t[p].r=r;
    if(l==r)
    {
        t[p].sum=a[l];
        t[p].lx=a[l];
        t[p].rx=a[l];
        t[p].dat=a[l];
        return;
    }
    int mid=(l+r)>>1;
    build(p*2,l,mid);
    build(p*2+1,mid+1,r);
    t[p].sum=t[p*2].sum+t[p*2+1].sum;
    t[p].lx=max(t[p*2].lx, t[p*2].sum+t[p*2+1].lx);
    t[p].rx=max(t[p*2+1].rx, t[p*2].rx+t[p*2+1].sum);
    t[p].dat=Max(t[p*2].dat, t[p*2+1].dat, t[p*2].rx+t[p*2+1].lx);
}
void change(int p, int x, int v)
{
    if(t[p].l==t[p].r)
    {
        t[p].sum=v;
        t[p].lx=v;
        t[p].rx=v;
        t[p].dat=v;
        //cout<<p<<"ass "<<t[p].dat<<endl;
        return ;
    }
    int mid=(t[p].l+t[p].r)>>1;
    if(x<=mid)change(p*2, x, v);
    else change(p*2+1,x,v);
    t[p].sum=t[p*2].sum+t[p*2+1].sum;
    t[p].lx=max(t[p*2].lx, t[p*2].sum+t[p*2+1].lx);
    t[p].rx=max(t[p*2+1].rx, t[p*2].rx+t[p*2+1].sum);
    t[p].dat=Max(t[p*2].dat, t[p*2+1].dat, t[p*2].rx+t[p*2+1].lx);
    //cout<<p<<"dd "<<t[p].dat<<endl;
}

segment_tree ask(int p, int l, int r)
{
    if(l<=t[p].l&&r>=t[p].r)return t[p];
    int mid=(t[p].l+t[p].r)>>1;
    if(r<=mid)return ask(p*2,l,r);
    else if(l>mid)return ask(p*2+1,l,r);
    else
    {
        segment_tree a=ask(p*2,l,r);
        segment_tree b=ask(p*2+1,l,r);
        segment_tree c;
        c.sum=a.sum+b.sum;
        c.lx=max(a.lx, a.sum+b.lx);
        c.rx=max(b.rx, a.rx+b.sum);
        c.dat=Max(a.dat, b.dat, a.rx+b.lx);
        return c;
    }
}
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)scanf("%lld",&a[i]);
    build(1,1,n);
    int k;
    int x,y;
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d%d",&k,&x,&y);
        if(k==2)
        {
            change(1,x,y);
        }
        if(k==1)
        {
            int temp1=min(x,y);
            int temp2=max(x,y);
            x=temp1;
            y=temp2;
            segment_tree ans=ask(1,x,y);
            printf("%lld\n",ans.dat);
        }
    }
    return 0;
}

写线段树一定要心平气和!!!理解每次递归所返回的信息能不能对所求的东西有所帮助。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值