(poj step8.1.1.1)POJ 1106 Transmitters(判断一条边在另一条边的哪一个方向——计算叉积)

/*
 * POJ_1106.cpp
 *
 *  Created on: 2013年11月16日
 *      Author: Administrator
 */

#include <iostream>
#include <cstdio>
#include <cmath>

using namespace std;

const double epsi = 1e-10;
const double pi = acos(-1.0);
const int maxn = 50005;

struct Point{
	double x,y;
	Point(double _x = 0,double _y = 0):x(_x),y(_y){

	}

	Point operator-(const Point& op2)const{//**要注意这种写法
		return Point(x - op2.x,y - op2.y);
	}

	double operator^(const Point& op2)const{
		return x*op2.y - y*op2.x;
	}
};

int sign(const double& x){//判断x是正数还是负数
	if(x > epsi){
		return 1;//传入的数是一个整数...
	}else if(x < -epsi){
		return -1;//传入的数是一个负数...
	}

	return 0;
}

double sqr(double x){//计算一个数的平方数
	return x*x;
}

double mul(const Point& p0,const Point& p1,const Point& p2){//计算p0p1与p1p2的叉积叉积
	return (p1-p0)^(p2-p0);
}

double dis2(const Point& p0,const Point& p1){//返回两个点的距离的平方
	return sqr(p0.x - p1.x) + sqr(p0.y - p1.y);
}

double dis(const Point& p0,const Point& p1){//返回两个点的距离
	return sqrt(dis2(p0,p1));
}

int n;
Point p[maxn],cp;
double r;

int main(){
	while(scanf("%lf%lf%lf",&cp.x,&cp.y,&r),r >= 0){
		scanf("%d",&n);

			int i;
			for(i = 0 ; i < n ; ++i){
				scanf("%lf%lf",&p[i].x,&p[i].y);
			}

			int j;
			int ans = 0;
			for(i = 0 ; i < n ; ++i){
				int temp = 0 ;
				for(j = 0 ; j < n ; ++j){
					if(sign(dis(p[j],cp) - r) != 1){//如果在辐射范围内
						if(sign(mul(cp,p[i],p[j])) != 1){//如果在辐射的一边,这里等于-1也是行的
							temp++;
						}
					}
				}

				ans = max(ans,temp);
			}

			printf("%d\n",ans);
	}

	return 0;
}