Helvetic Coding Contest 2019 online mirror C2 - Heidi and the Turing Test(线段树/扫描线+坐标系旋转)

题目

给你n(n<=3e5)个点，让你自选(x,y)，可以不是这n个点的点，也可以不是整点

画一个以其为中心，距其曼哈顿距离不超过r(1<=r<=1e6)的图形

要求圈住最多的点，输出圈住的点的个数，-1e6<=xi,yi<=1e6

思路来源

http://dj3500.webfactional.com/helvetic-coding-contest-2019-editorial.pdf

http://blog.sina.com.cn/s/blog_3fd642cf0101cc8w.html

题解

官方题解，先坐标系旋转逆时针-45度+伸缩变换√2倍+y轴镜像，把正方形边长从√2r搞成2r，然后扫描线裸题

我只把坐标轴逆时针转了-45度+伸缩变换√2倍，

所以，(x,y)新坐标rotate一下是(x-y,x+y)，离散化竖线y的坐标建线段树跑扫描线即可

 

坐标系旋转，记向逆时针旋转度数为正，，

记矩阵形式就好，系数矩阵的模为1，还是比较好记的

代码

#include<bits/stdc++.h>
using namespace std;
const int N=6e5+10;
//插入点(x,y) 插入竖线 [y,y+h 区间+1
//删除点(x+w,y) 删除竖线 [y,y+h] 区间-1
//对纵坐标离散化后建树 
//按x坐标增序插入点 询问某一时刻区间最大值  
struct node
{
	int x,y,yy,v;//一条(x,y)到(x,yy)的竖直线 +1 -1 
}b[N];
bool operator<(node a,node b)
{
	if(a.x==b.x)return a.v>b.v;//先插入 后删除
	return a.x<b.x;
}
void rotate(int &x,int &y)//坐标系旋转-45度 长度伸长为原来根号2倍 (x-y,x+y) 
{
	int xx=x,yy=y;
	x=xx-yy;
	y=xx+yy;
}
int n,k,x,y,xx,yy;
int Y[N],cnt;
int ans;
int mx[4*N],cov[4*N];
void pushdown(int p,int l,int r)
{
	if(cov[p])
	{
		mx[p<<1]+=cov[p];
		cov[p<<1]+=cov[p];
		mx[p<<1|1]+=cov[p];
		cov[p<<1|1]+=cov[p];
		cov[p]=0;
	} 
}
void update(int p,int l,int r,int ql,int qr,int v)
{
	if(ql<=l&&r<=qr)
	{
		mx[p]+=v;
		cov[p]+=v;
		return;
	}
	pushdown(p,l,r);
	int mid=(l+r)/2;
	if(ql<=mid)update(p<<1,l,mid,ql,qr,v);
	if(qr>mid)update(p<<1|1,mid+1,r,ql,qr,v);
	mx[p]=max(mx[p<<1],mx[p<<1|1]);
}
int ask(int p,int l,int r,int ql,int qr)
{
	if(ql<=l&&r<=qr)return mx[p];
	pushdown(p,l,r);
	int mid=(l+r)/2,ans=0;
	if(ql<=mid)ans=max(ans,ask(p<<1,l,mid,ql,qr));
	if(qr>mid)ans=max(ans,ask(p<<1|1,mid+1,r,ql,qr));
	return ans;
}
int main()
{
	scanf("%d%d",&n,&k);
	for(int i=1;i<=n;++i)
	{
		scanf("%d%d",&x,&y);
		xx=x-k;yy=y;//正方形下顶点 
		rotate(xx,yy);//坐标系变换 逆时针转轴45度 伸缩根号2倍 
		b[2*i-1]=node{xx,yy,yy+2*k,1};
		xx=x+k;yy=y;//正方形上顶点 
		rotate(xx,yy);
		b[2*i]=node{xx,yy-2*k,yy,-1};
		Y[++cnt]=yy-2*k;
		Y[++cnt]=yy;
	}
	sort(Y+1,Y+cnt+1);
	cnt=unique(Y+1,Y+cnt+1)-(Y+1);
	for(int i=1;i<=2*n;++i)
	{
		b[i].y=lower_bound(Y+1,Y+cnt+1,b[i].y)-Y;
		b[i].yy=lower_bound(Y+1,Y+cnt+1,b[i].yy)-Y;
	}
	sort(b+1,b+2*n+1);
	for(int i=1;i<=2*n;++i)
	{
		update(1,1,cnt,b[i].y,b[i].yy,b[i].v);
		ans=max(ans,ask(1,1,cnt,1,cnt));
	}
	printf("%d\n",ans); 
	return 0;
}