ACdream 1127(Base Station-树状数组-2个约束条件)

Base Station

Time Limit: 20000/10000MS (Java/Others)Memory Limit: 512000/256000KB (Java/Others)

Submit Statistic Next
Problem

Problem Description

移动通信系统中。通信网的建立主要通过基站来完毕。

基站能够分为主基站和子基站。子基站和各个移动用户进行连接，子基站必须通过主基站来和外界实现通信。主基站能够覆盖到的范围是一个圆形区域，子基站和主基站的距离小于半径r才干被该主基站覆盖到。半径r由主基站的发射功率确定。

某个区域的移动通信网，包括2个主基站和N个子基站。它们的位置都能够相应到一个整数坐标上。假设子基站至少被一个主基站覆盖。则该子基站是激活的。

如今通信公司在调试设备，它们不停地改变主基站的发射功率。当两个主基站的覆盖半径分别为r1和r2时，须要知道有多少个子基站处于非激活状态。

Input

有若干组输入数据。

第一行是四个整数：x1、y1、x2、y2（1<=x1、y1、x2、y2<=10^9），表示两个主基站的坐标是（x1,y1）和（x2,y2）。

第二行是一个整数N（1<=N<=100000），表示有N个子基站。

接下来的N行。每行两个整数x、y(1<=x, y<=10^9)。表示每一个子基站的坐标。

接下来一行包括一个整数M（1<=M<=100000），表示有M个询问。

接下来的M行。每行两个整数r1、r2(1<=r1, r2<=10^9)。表示询问当两个主基站的覆盖半径为r1和r2时。处于非激活状态的子基站数。

Output

对每一个查询，输出答案。

Sample Input

Sample Output

Hint

Source

kuangbin

Manager

Dshawn

将基站到2个主站的距离^2表示成二维坐标

则本题的询问表示成半径^2

本题等价于“平面上有n个点，问n-横坐标<r1或纵坐标<r2的点数".

用树状数组维护。先依照x从小到大插入，询问y值，得到左下角点数->左上角点数。

接着用容斥。求右上角的点数

#include<cstdio>

#include<cstring>

#include<cstdlib>

#include<algorithm>

#include<functional>

#include<iostream>

#include<cmath>

#include<cctype>

#include<ctime>

using namespace std;

#define For(i,n) for(int i=1;i<=n;i++)

#define Fork(i,k,n) for(int i=k;i<=n;i++)

#define Rep(i,n) for(int i=0;i<n;i++)

#define ForD(i,n) for(int i=n;i;i--)

#define RepD(i,n) for(int i=n;i>=0;i--)

#define Forp(x) for(int p=pre[x];p;p=next[p])

#define Lson (x<<1)

#define Rson ((x<<1)+1)

#define MEM(a) memset(a,0,sizeof(a));

#define MEMI(a) memset(a,127,sizeof(a));

#define MEMi(a) memset(a,128,sizeof(a));

#define INF (2139062143)

#define F (100000007)

#define MAXN (100000+10)

#define MAXXi (1000000000)

#define MAXM (100000+10)

long long mul(long long a,long long b){return (a*b)%F;}

long long add(long long a,long long b){return (a+b)%F;}

long long sub(long long a,long long b){return (a-b+(a-b)/F*F+F)%F;}

typedef long long ll;

int n,X1,Y1,X2,Y2,m,tot;

int lowbit(int x){return x&(-x);}

ll sqr(ll x){return x*x;}

ll dis2(ll x1,ll y1,ll x2,ll y2){return sqr(x1-x2)+sqr(y1-y2);}

struct arr_tree

{

	int a[MAXN+MAXM*6];

	void reset(){ MEM(a) }

	void add1(int x)

	{

		for(;x<=tot;x+=lowbit(x)) a[x]++;

	}

	int sum(int x)

	{

		int ans=0;

		for(;x>0;x-=lowbit(x)) ans+=a[x];

		return ans;

	}

}T;

struct node

{

	ll x,y; //距离^2

	friend bool operator<(node a,node b){return a.x<b.x;	}

}a[MAXN*2];

struct comm

{

	ll r1,r2; //半径^2

	int i;

	friend bool operator<(comm a,comm b){return a.r1<b.r1;	}

}ask[MAXM*2];

int ans[MAXM*2];

ll y_[MAXN+MAXM*2];

int loc(ll y)

{

	return lower_bound(y_+1,y_+1+tot,y)-y_;

}

int main()

{

	while (scanf("%d%d%d%d",&X1,&Y1,&X2,&Y2)==4)

	{

		T.reset();

		scanf("%d",&n);

		For(i,n)

		{

			int x,y;

			scanf("%d%d",&x,&y);

			a[i].x=dis2(x,y,X1,Y1);a[i].y=dis2(x,y,X2,Y2);

			y_[i]=a[i].y;

		}

		sort(a+1,a+1+n);

		scanf("%d",&m);

		For(i,m)

		{

			int r1,r2;

			scanf("%d%d",&r1,&r2);

			ask[i].i=i;ask[i].r1=sqr(r1);ask[i].r2=sqr(r2);

			y_[n+i]=ask[i].r2;

		}

		sort(ask+1,ask+1+m);	

		sort(y_+1,y_+n+m+1);

		tot=unique(y_+1,y_+n+m+1)-(y_+1);

		int j=1;

		For(i,m)

		{

			while(j<=n&&a[j].x<ask[i].r1) T.add1(loc(a[j].y)),j++;

			ans[ask[i].i]=j-1-T.sum(loc(ask[i].r2)-1);

		}

		while(j<=n) T.add1(loc(a[j].y)),j++;

		For(i,m) ans[ask[i].i]+=T.sum(loc(ask[i].r2)-1);

		For(i,m) printf("%d\n",n-ans[i]);

	}

	return 0;

}

巴特西