#include <iostream>
#include <cmath>
#include <algorithm>
#define pi 3.141592653589793238
using namespace std;
int main()
{
	freopen("in.txt","r",stdin);
	double x,y,r,ax,ay,temp;
	int dot,t,i,j,max,sum,sumofcentre;
	while(cin>>x>>y>>r&&r>=0)
	{
		double node[1000]={0};
		cin>>dot;
		for(t=sumofcentre=0;dot;--dot)
		{
			cin>>ax>>ay;
			if((ax-x)*(ax-x)+(ay-y)*(ay-y)<=r*r)
			{
				temp=atan((ay-y)/(ax-x));
				if(ax-x>0)
				{
					if(ay-y>=0)node[t++]=temp;
					else node[t++]=2*pi+temp;
				}
				else if(ax-x==0)
				{
					if(ay-y>0)node[t++]=pi/2;
					else if(ay-y==0)++sumofcentre;
					else node[t++]=3*pi/2;
				}
				else node[t++]=pi+temp;
			}
		}
		sort(node,node+t);
		for(i=max=0;i<t-1;++i)
		{
			for(j=i+1,sum=1+sumofcentre;j<t;++j)
			{
				if(node[j]-node[i]<=pi)
					++sum;
				else break;
			}
			if(sum>max)max=sum;
		}
		cout<<max<<endl;
	}
	return 0;
}

• 用向量叉积来做,不会有浮点误差. (1033) 00448242 2005-07-24 12:01:34
  • 他自己写了一个atan2 (4) frkstyc 2005-07-24 12:02:54
    • 嗯,能用向量积的,打死我也不用三角函数. (0) 00448242 2005-07-24 12:04:06
      • 来啊，拖出去打！ (0) frkstyc 2005-07-24 12:05:06
        • 啊~~~~什么时候成县太爷了~~~~ (0) 00448242 2005-07-24 12:05:58
          • 是Stephen Chow (0) frkstyc 2005-07-24 12:08:54
  • Re:用向量叉积来做,不会有浮点误差. (47) 00448242 2005-07-24 12:03:08
    • 我也这么写了，不过显然我没这么交 (51) frkstyc 2005-07-24 12:03:46
• 最近几何题比较火 (1031) sunmoonstar_love 2005-07-24 12:53:08
  • 某bt布置了1X个几何题，不然才不跟他折腾呢 (0) frkstyc 2005-07-24 12:54:57
    • 几何题往往是弱项，多练练也是有好处，一见几何就犯晕，现在能分析一下了，虽然很多分析好了，也不能AC (0) sunmoonstar_love 2005-07-24 13:01:05
      • re,一见几何就犯晕 (0) 00448242 2005-07-24 13:08:54
        • re (0) frkstyc 2005-07-24 13:09:39
        • 我看这道题目还是很规矩的啊!用ACOS(X)还是挺舒服的! 比图论的题目强多了!! (0) ook 2005-07-24 14:02:52