[kuangbin带你飞]专题七 线段树Y - Stars in Your Window POJ - 2482 扫描线

Y - Stars in Your Window

 POJ - 2482

题目很美。。。

扫描线类型的题。

有一个比较好的线段树思想。

给你一个数列  让你求连续数字和（满足连续数字个数小于等于w）最大。是多少

对每个数a，（原始数列）坐标为id。。：插入a到线段树对应位置x中，再把-a插入到x+w上。

每次操作完记录st【1】的值求最大就行。

注意对每个位置  先插入负数，再插入正数（先消去 x-w位置数 的影响，因为这之前的数--不满足连续数字个数到当前位置小于w这个条件）

 

这一题也是这个思想，我们只考虑矩形的右上角坐标的位置。

对每个星星来说     如果矩形右上角坐标在（x,y）--》(x+w-1,x+h-1)这个范围内，那这个矩形就可以围住这个星星。

那么这一题就等价于 给定很多区域（每个星星有一个矩形区域），每个区域有一个权值，这些区域并中的坐标权值最大的是多少。

我们可以采用扫描线的思想（如不会自行百度），从左往右扫描，每次先加入正数，再加入负数（或者把矩形横坐标扩大1，先加入负数）总之这个顺序是为了消除边界的影响。。

然后对每条线，进行区间更新，整体区间求最大值即ans[1]。再求每次的最大即可。

#include<iostream>
#include<cmath>
#include<algorithm>
#include<queue>
#include<cstring>
#include<cstdio>
#include<vector> 
#include<map>
using namespace std;
typedef long long ll;
const int M = 100000+7;
struct node
{
	ll x,y1,y2;
	ll c;
}rem[M];
vector<ll>v;
int getid(ll x)
{
	return lower_bound(v.begin(),v.end(),x)-v.begin()+1;
}
bool cmp(node a, node b)
{
	if(a.x==b.x)
	return a.c>b.c;
	return a.x<b.x;
}
ll st[M<<2],lazy[M<<2];
void pushup(int l,int r,int rt)
{
	if(lazy[rt])
	{
		int mid=(l+r)>>1;
		lazy[rt<<1]+=lazy[rt];
		lazy[rt<<1|1]+=lazy[rt];
		st[rt<<1]+=lazy[rt];
		st[rt<<1|1]+=lazy[rt];
		lazy[rt]=0;
}
}
void update(int l,int r,int rt,int x,int y,ll c)
{
	//printf("%d  %d  %d  %d\n",l,r,x,y);
	if(x<=l&&r<=y)
	{
		lazy[rt]+=c;
		st[rt]+=c;
		return ;
	}
	pushup(l,r,rt);
	int mid=(l+r)>>1;
	if(x<=mid)update(l,mid,rt<<1,x,y,c);
	if(y>mid)update(mid+1,r,rt<<1|1,x,y,c);
	st[rt]=max(st[rt<<1],st[rt<<1|1]);
}
int main()
{
	ll n,w,h,x,y,c,cnt=0;
	while(~scanf("%lld%lld%lld",&n,&w,&h))
	{
		cnt=0;
		memset(lazy,0,sizeof(lazy));
		memset(st,0,sizeof(st));
		v.clear();
		for(int i=1;i<=n;i++)
			scanf("%lld%lld%lld",&x,&y,&c),v.push_back(y),v.push_back(y+h-1),
			rem[++cnt].x=x,rem[cnt].y1=y,rem[cnt].y2=y+h-1,rem[cnt].c=c,
			rem[++cnt].x=x+w-1,rem[cnt].y1=y,rem[cnt].y2=y+h-1,rem[cnt].c=-c;
		sort(v.begin(),v.end()),v.erase(unique(v.begin(),v.end()),v.end());//离散化去重，很方便 
		sort(rem+1,rem+1+cnt,cmp);
		ll ans=0;
		for(int i=1;i<=cnt;i++)
		{
			update(1,cnt,1,getid(rem[i].y1),getid(rem[i].y2),rem[i].c);
			//printf("%d----%d\n",st[1],rem[i].x);
			ans=max(ans,st[1]);
		}
		printf("%lld\n",ans);
	}
	return 0;
 }