poj-2482

最新推荐文章于 2021-02-20 10:33:09 发布

原创最新推荐文章于 2021-02-20 10:33:09 发布 · 247 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#C++ #POJ

【12】编程网站C++题解专栏收录该内容

36 篇文章

订阅专栏

本文介绍了一个基于扫描线算法结合线段树实现的二维区间更新查询问题解决方案。通过使用C++编程语言，该方案能够高效地处理一系列矩形区域的增加和查询操作，找到所有矩形覆盖下的最大值。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
using namespace std;
#define maxn 110000
#define LL long long
struct list
{
    LL l,r,x,val;
} node[maxn*6];
struct liss
{
    LL x;
    LL y1,y2;
    LL c,val;
} line[maxn*2];
LL lines;
LL yy[maxn];
LL a[maxn],b[maxn],c[maxn];
LL sk;
LL n,W,H;
LL maxx;
int cmp(struct liss a,struct liss b)
{
    if(a.x!=b.x)return a.x<b.x;
    return a.val>b.val;
}
LL search(LL x)
{
    LL l=0;
    LL r=sk;
    LL mid=(l+r)/2;
    while(l<r)
    {
        if(yy[mid]==x)return mid;
        if(yy[mid]>x)r=mid;
        if(yy[mid]<x)l=mid+1;
        mid=(l+r)/2;
    }


}
void build(LL l,LL r,int st)
{
    LL mid=(l+r)/2;
    node[st].l=l;
    node[st].r=r;
    node[st].x=0;
    node[st].val=0;
    if(l!=r)
    {
        build(l,mid,st*2);
        build(mid+1,r,st*2+1);
    }
}
void insert(LL l,LL r,LL val,int st)
{
    LL ll=node[st].l;
    LL rr=node[st].r;
    if(l>rr||r<ll)return;
    if(l<=ll&&rr<=r)
    {
        node[st].x+=val;
        node[st].val+=val;
        return ;
    }
    insert(l,r,val,st*2);
    insert(l,r,val,st*2+1);
    node[st].val=max(node[st*2].val,node[st*2+1].val)+node[st].x;
}
int main()
{
    LL i;
    while(~scanf("%lld%lld%lld",&n,&W,&H))
    {


        LL sy=0;
        LL yy1,yy2,xx1,xx2;
        for(i=0; i<n; i++)
        {
            scanf("%lld%lld%lld",&a[i],&b[i],&c[i]);
            yy1=b[i];
            yy[sy++]=yy1;
            yy2=b[i]+H-1;
            yy[sy++]=yy2;
        }
        if(W==0||H==0)
        {
            cout<<"0"<<endl;
            continue;
        }
        H--,W--;
        sort(yy,yy+sy);
        sk=1;
        for(i=1; i<sy; i++)if(yy[i]!=yy[i-1])yy[sk++]=yy[i];
        lines=0;
        for(i=0; i<n; i++)
        {
            yy2=b[i]+H;
            yy1=b[i];
            xx2=a[i]+W;
            xx1=a[i];
            yy1=search(yy1);
            yy2=search(yy2);
            line[lines].x=xx1;
            line[lines].y1=yy1;
            line[lines].y2=yy2;
            line[lines++].val=c[i];
            line[lines].x=xx2;
            line[lines].y1=yy1;
            line[lines].y2=yy2;
            line[lines++].val=-c[i];
        }
        sort(line,line+lines,cmp);
        build(0,sk,1);
        maxx=0;
        for(i=0; i<lines; i++)
        {
            insert(line[i].y1,line[i].y2,line[i].val,1);
            maxx=max(maxx,node[1].val);
        }
        cout<<maxx<<endl;
    }
    return 0;
}