记录编号	26919	评测结果	AAAAAAAAAA
题目名称	打蚊子	最终得分	100
用户昵称	苏轼	是否通过	通过
代码语言	C++	运行时间	1.914 s
提交时间	2011-07-29 16:22:57	内存使用	0.30 MiB

#include <cstdio>
#include <cmath>
#include <vector>
#include <algorithm>
using namespace std;

typedef long long Int;
const int MAXN=5005;
const double PI=acos(-1);
typedef pair<double,int> Pair;
#define pb push_back
#define mp make_pair

struct Point
{
	int x,y;
}P[MAXN];

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

Int dis2(const Point &a,const Point &b)
{
	return sqr(a.x-b.x)+sqr(a.y-b.y);
}

int N,R;
int re;
vector<Pair> v;

int main()
{
	freopen("fight.in","r",stdin);
	freopen("fight.out","w",stdout);
	scanf("%d%d",&N,&R);
	for(int i=0;i<N;i++)
		scanf("%d%d",&P[i].x,&P[i].y);
	random_shuffle(P,P+N);
	for(int i=0;i<N;i++)
	{
		v.clear();
		double t;
		int x=0;
		for(int j=0;j<N;j++)
			if (i!=j && (t=dis2(P[i],P[j]))<=(Int)4*R*R)
			{
				double alpha=atan2(P[j].y-P[i].y,P[j].x-P[i].x);
				if (alpha<0)
					alpha+=2*PI;
				double beta=acos(sqrt(t)/(2.0*R));
				v.pb(mp(alpha-beta+2*PI,0));
				v.pb(mp(alpha+beta+2*PI,1));
				x++;
			}
		if (x<=re)
			continue;
		sort(v.begin(),v.end());
		int sum=0;
		for(vector<Pair>::iterator it=v.begin();it!=v.end();it++)
			if (it->second)
				sum--;
			else
			{
				sum++;
				if (sum>re)
					re=sum;
			}
	}
	printf("%d\n",re+1);
	return 0;
}