luoguP1502 窗口的星星（扫描线）20180305 T1

题目

传送门
假定天空是一个平面，每个星星都有一个坐标(x,y)，每颗星星都有一个亮度C，代表它的亮度。窗户是长方形的，有固定的长和宽，边平行于x轴和y轴。你的任务是告诉我如何摆放窗户，才能获得在窗口内所有星星的亮度总和最大值。注意，边框的星星不算。窗口可以被平移，但不允许旋转。
输入
在输入中有几个测试用例。每一行的第一行包含3个整数：n，w，h，表示星形的数目，矩形窗口的水平长度和垂直高度。然后N行如下，有3个整数分别为：X，Y，C，表示星星位置（x，y）和每个恒星的亮度。没有两颗星星在同一点上。
天空中至少有1颗，最多有10000颗星星。1 =

题解

扫描线好题

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
const int maxn=100001;

int cnt,maxnow;
int n,w,h,le[maxn],ri[maxn],add[maxn<<2],v[maxn<<2];
struct Star{
    int x,y,dis,id;
    bool f;
}a[maxn],b[maxn];

bool cmp1(Star x,Star y) {return x.x<y.x;}
bool cmp2(Star x,Star y) {return x.y<y.y;}

void pushdown(int rt)
{
    if (add[rt])
    {
        add[rt<<1]+=add[rt];
        add[rt<<1|1]+=add[rt];
        v[rt<<1]+=add[rt];
        v[rt<<1|1]+=add[rt];
        add[rt]=0;
    }
}

void change(int L,int R,int C,int l,int r,int rt)
{
    if (L<=l && r<=R)
    {
        add[rt]+=C; v[rt]+=C; return;
    }
    int mid=(l+r)>>1;
    pushdown(rt);
    if (L<=mid) change(L,R,C,l,mid,rt<<1);
    if (R>mid) change(L,R,C,mid+1,r,rt<<1|1);
    v[rt]=max(v[rt<<1],v[rt<<1|1]);
}

int ques(int L,int R,int l,int r,int rt) 
{
    if (L<=l && r<=R) return v[rt];
    int mid=(l+r)>>1;
    pushdown(rt);
    int ans=0;
    if (L<=mid) ans=max(ans,ques(L,R,l,mid,rt<<1));
    if (R>mid) ans=max(ans,ques(L,R,mid+1,r,rt<<1|1));
    return ans;
}

int main()
{
    freopen("stars.in","r",stdin);
    freopen("stars.out","w",stdout);

    scanf("%d%d%d",&n,&w,&h);
    for (int i=1; i<=n; i++)
    {
        scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].dis);
        a[i].id=i; a[i].f=0;
        b[i].y=a[i].y; b[i].dis=a[i].dis; b[i].id=i;
    }
    int tot=n;
    for (int i=1; i<=n; i++)
    {
        a[++tot].x=a[i].x+w-1;
        a[tot].id=i;
        a[tot].f=1;
    }
    sort(a+1,a+1+tot,cmp1);
    for (int i=1; i<=tot; i++)
    {
        int tt=0;
        if (i!=1 && a[i].x==a[i-1].x) tt=cnt;
        else tt=++cnt;
        if (!a[i].f) le[a[i].id]=tt;
        else ri[a[i].id]=tt;
    }
//      for (int i=1; i<=n; i++) printf("%d %d\n",le[a[i].id],ri[a[i].id]);
    tot=n;
//      for (int i=1; i<=n; i++) printf("%d %d \n",b[i].y,b[i].id); 
    for (int i=1; i<=n; i++)
    {
        b[++tot].y=b[i].y+h-1;
        b[tot].id=i;
        b[tot].dis=-b[i].dis;//
    }
    sort(b+1,b+1+tot,cmp2);
//      for (int i=1; i<=tot; i++) printf("%d %d %d\n",b[i].y,b[i].id,b[i].dis);
//      printf("%d ",cnt); return 0;
    for (int i=1; i<=tot; i++)
    {
        if (b[i].dis>0)
        {
            int num=ques(le[b[i].id],ri[b[i].id],1,cnt,1);
            if (num+b[i].dis>maxnow) maxnow=num+b[i].dis;
        }
        change(le[b[i].id],ri[b[i].id],b[i].dis,1,cnt,1);
    }
    printf("%d\n",maxnow);
    return 0;
}

写过的题不能写不对