[SHOI2007] 园丁的烦恼

• 二维数点问题
• 我们通过扫描线可以将静态的二维问题转换为动态的一维问题
• 维护动态的一维问题就使用数据结构维护序列

点击查看代码

#include <bits/stdc++.h>
using namespace std;
struct t1
{
	int x,y;
}t[500005];
struct q1
{
	int x1,y1,x2,y2;
}q[500005];
struct l1
{
	int x,y1,y2,id,opt;
}l[1000005];
int c[10000005],ans[500005];
int read1()
{
	char cc=getchar();
	while(!(cc>=48&&cc<=57))
	{
		if(cc=='-')
		{
			break;
		}
		cc=getchar();
	}
	bool f=false;
	int s=0;
	if(cc=='-')
	{
		f=true;
	}
	else
	{
		s=cc-48;
	}
	while(1)
	{
		cc=getchar();
		if(cc>=48&&cc<=57)
		{
			s=s*10+cc-48;
		}
		else
		{
			break;
		}
	}
	if(f==true)
	{
		s=-s;
	}
	return s;
}
int lowbit(int n)
{
	return n&(-n);
}
void add(int n1)
{
	if(n1==0)
	{
		c[0]++;
		return;
	}
	while(n1<=10000000)
	{
		c[n1]++;
		n1+=lowbit(n1);
	}
}
int ask(int n1)
{
	int s=0;
	while(n1>0)
	{
		s+=c[n1];
		n1-=lowbit(n1);
	}
	return s+c[0]*(n1>=0);
}
bool cmpt(t1 a,t1 b)
{
	return a.x<b.x;
}
bool cmpl(l1 a,l1 b)
{
	return a.x<b.x;
}
int main()
{
	int n,m;
	cin>>n>>m;
	for(int i=1;i<=n;i++)
	{
		t[i].x=read1();
		t[i].y=read1();
	}
	for(int i=1;i<=m;i++)
	{
		q[i].x1=read1();
		q[i].y1=read1();
		q[i].x2=read1();
		q[i].y2=read1();
		l[i].x=q[i].x1-1;
		l[i+m].x=q[i].x2;
		l[i].y1=l[i+m].y1=q[i].y1;
		l[i].y2=l[i+m].y2=q[i].y2;
		l[i].id=l[i+m].id=i;
		l[i].opt=-1;
		l[i+m].opt=1;
	}
	sort(t+1,t+n+1,cmpt);
	sort(l+1,l+2*m+1,cmpl);
	int j=0;
	for(int i=1;i<=2*m;i++)
	{
		while(j+1<=n&&t[j+1].x<=l[i].x)
		{
			j++;
			add(t[j].y);
		}
		ans[l[i].id]+=(l[i].opt*(ask(l[i].y2)-ask(l[i].y1-1)));
	}
	for(int i=1;i<=m;i++)
	{
		printf("%d\n",ans[i]);
	}
	return 0;
}