Snowy Smile 【多校6 HDU6638】【线段树维护最长字段和】

题目链接

题目大意

有n个点，每个点有一个坐标和一个权值（权值可正可负），让你找出一个矩形，矩形要平行于坐标轴，矩形里所有的点的权值和最大。

解题思路

首先把坐标都离散化，因为之后2000个点，点的数目比较少。
因为不知道矩形的长和宽，所以我们开始枚举矩形的左右边界，把点按照x从小到大排好序，然后可以开始枚举它的上边界和下边界，之后我们把每次添加的边界里的点都加进去，用线段树维护边界里面的y值对应的权值，每次找到最大的权值和。
注意这个边应该是一排一排的判断，当x值相同的时候，这个边界里x相同的值加进去完了之后才可以判断，y也是一样。

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;

const int maxn=2005;
int yy[maxn],b[maxn];
struct node
{
    int x,y,w;
};
node a[maxn];
struct Node
{
    long long ll,lr,lm,sum;
    int l,r;
};
Node e[maxn*4];

int cmp(node n,node m)
{
    return n.x<m.x;
}
void build(int root,int l,int r)
{
    e[root].l=l;
    e[root].r=r;
    e[root].ll=e[root].lr=e[root].lm=e[root].sum=0;
    if(e[root].l==e[root].r)
        return ;
    build(root*2,l,(l+r)/2);
    build(root*2+1,(l+r)/2+1,r);
}
void updown(int root)
{
    e[root].sum=e[root*2].sum+e[root*2+1].sum;
    e[root].ll=max(e[root*2].ll,e[root*2].sum+e[root*2+1].ll);
    e[root].lr=max(e[root*2+1].lr,e[root*2+1].sum+e[root*2].lr);
    e[root].lm=max(max(e[root*2].lm,e[root*2+1].lm),e[root*2].lr+e[root*2+1].ll);
}
void insertt(int root,int x,int w)
{
    if(e[root].l==e[root].r)
    {
        e[root].ll=e[root].lr=e[root].lm=e[root].sum=e[root].lm+w;
        return ;
    }
    int mid=(e[root].l+e[root].r)/2;
    if(x<=mid)
        insertt(root*2,x,w);
    else
        insertt(root*2+1,x,w);
    updown(root);
}
int main()
{
    int T,n,x,y,z;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        for(int i=1; i<=n; i++)
        {
            scanf("%d %d %d",&a[i].x,&a[i].y,&a[i].w);
            yy[i]=a[i].y;
        }
        sort(a+1,a+1+n,cmp);
        sort(yy+1,yy+1+n);
        int mm=unique(yy+1,yy+1+n)-yy-1;
        long long ans=0;
        for(int i=1; i<=n; i++)
        {
            if(i!=1&&a[i].x==a[i-1].x)
                continue;
            build(1,1,mm);
            for(int j=i; j<=n; j++)
            {
                if(j!=i&&a[j].x!=a[j-1].x)//上一次x的值相同的一排都加进去之后再判断
                    ans=max(ans,e[1].lm);
                int aa=lower_bound(yy+1,yy+1+mm,a[j].y)-yy;
                insertt(1,aa,a[j].w);
            }
            ans=max(ans,e[1].lm);
        }
        printf("%lld\n",ans);
    }
    return 0;
}