hdu 6638 Snowy Smile

原文链接：https://blog.csdn.net/wu_tongtong/article/details/73385029

线段树的区间最大子段和：

线段树需要维护的是： 
左端点 x 
右端点 y  （本人喜欢直接维护端点） 
[x,y]内的最大子段和 ms 
[x,y]的区间和 s 
[x,y]内的紧靠左端点的最大子段和 ls 
[x,y]内的紧靠右端点的最大子段和 rs

困难就是，update和ask（l,r）询问[l,r]区间内的最大子段和 
那我们一步一步来

s的维护很常规， 
ls：有两种情况： 
1.该区间内的ls是ta左儿子的ls 
2.该区间内的ls是左儿子的s+右儿子的ls 
同理，rs：有两种情况： 
1.该区间内的rs是ta右儿子的rs 
2.该区间内的rs是右儿子的s+左儿子的rs 
而ms有三种情况： 
1.该区间内的ms是左儿子的ms 
2.该区间内的ms是右儿子的ms 
3.该区间内的ms是左儿子的rs+右儿子的ls 

这个题离散化坐标，以y为轴，然后再插入x就好了。

看了这个大佬的博客https://blog.csdn.net/A_Thinking_Reed_/article/details/98778260

 

#include <bits/stdc++.h>
using namespace std;
#define  LL long long
const int maxn=2005;
struct node{
    LL x,y,val;
}p[maxn];
struct tree{
    LL sum,sum_max,ls_max,rs_max;
}tree[maxn*4];

int t,n; vector<int>vx,vy;

int get_x(int x){ return (lower_bound(vx.begin(),vx.end(),x)-vx.begin())+1;}
int get_y(int x){ return (lower_bound(vy.begin(),vy.end(),x)-vy.begin())+1;}

bool cmp(node a,node b){
    if(a.y==b.y)
        return a.x<b.x;
    else
        return a.y<b.y;
}

void build_tree(int rt,int sx,int ex){
    tree[rt].ls_max=tree[rt].rs_max=tree[rt].sum=tree[rt].sum_max=0;
    if(sx==ex)
        return ;
    int mid=(sx+ex)>>1;
    build_tree(rt<<1,sx,mid);
    build_tree(rt<<1|1,mid+1,ex);
}
void pushup(int st){
    int rs=st<<1|1,ls=st<<1;
    tree[st].sum=tree[rs].sum+tree[ls].sum;
    tree[st].sum_max=max(max(tree[ls].sum_max,tree[rs].sum_max) ,tree[ls].rs_max+tree[rs].ls_max );
    tree[st].ls_max=max(tree[ls].ls_max , tree[ls].sum+tree[rs].ls_max);
    tree[st].rs_max=max(tree[rs].rs_max , tree[rs].sum+tree[ls].rs_max);

}

void update(int st,int l,int r,int pos,LL val){

    if(l==r){
        tree[st].sum+=val;
        tree[st].ls_max=tree[st].rs_max=tree[st].sum_max=tree[st].sum;
        return ;
    }

    int mid=(l+r)>>1;
    if(pos<=mid)
        update(st<<1,l,mid,pos,val);
    else
        update(st<<1|1,mid+1,r,pos,val);
    pushup(st);

}


int main( ) {
    freopen("in.txt","r",stdin);
    ios::sync_with_stdio(false);
    cin>>t;
    while(t--){
        cin>>n;
        vx.clear(); vy.clear();
        for(int i=1;i<=n;i++){
            cin>>p[i].x>>p[i].y>>p[i].val;
            vx.push_back(p[i].x);vy.push_back(p[i].y);
        }
        sort(vx.begin(),vx.end()),sort(vy.begin(),vy.end());
        vx.erase(unique(vx.begin(),vx.end()),vx.end());
        vy.erase(unique(vy.begin(),vy.end()),vy.end());

        for(int i=1;i<=n;i++){
            p[i].x=get_x(p[i].x);
            p[i].y=get_y(p[i].y);
        }
        sort(p+1,p+n+1,cmp);

        LL sum_max=0,last=-1,nx=vx.size();
        for(int i=1;i<=n;i++){ //遍历Y的底边
            if(p[i].y==last) continue;//相同的底边跳过
            build_tree(1,1,nx);
            for(int j=i,k=j;j<=n;j=k){//Y的上端
                while(k<=n&&p[k].y==p[j].y){
                    update(1,1,nx,p[k].x,p[k].val);
                    k++;
                }
                sum_max=max(sum_max,tree[1].sum_max);
            }
            last=p[i].y;
        }
        cout<<sum_max<<endl;

    }Snowy Smile

    return 0;
}