poj 1981 Circle and Points poj 2693 Chocolate Chip Cookies

题意：两个题都是让你求一个半径给定的圆最多能框住多少个点。

思路：枚举任意两个点，当做圆上的点，现在已经知道两个点，和圆的半径，可以求出圆心，然后再枚举每个点到圆心的距离是不是小于半径，统计个数。求圆心的时候用到了一些简单的解析几何的知识。都是很基础的。

两个题代码差不多我只贴1981的：

/* * File: main.cpp * Author: Mi * * Created on 2011年5月6日, 下午10:38 */ #include <stdio.h> #include <string.h> #include <algorithm> #include <math.h> #include <map> #include <vector> #include <stdlib.h> #define N 305 #define eps 1e-8 #define pi acos(double(-1)) using namespace std; /* * */ struct Point { double x,y; }p[N]; double distance(double x1,double y1,double x2,double y2) { return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)); } Point Get_center(Point p1,Point p2,double r) { Point p,center; double a,h,tx,ty; p.x=(p1.x+p2.x)/2.0,p.y=(p1.y+p2.y)/2.0; tx=p1.x-p2.x,ty=p1.y-p2.y; if(fabs(ty)<eps) a=pi/2; else a=atan(-tx/ty); h=sqrt(r*r-distance(p.x,p.y,p1.x,p1.y)*distance(p.x,p.y,p1.x,p1.y)); center.x=p.x+h*cos(a); center.y=p.y+h*sin(a); return center; } int main(int argc, char** argv) { int i=0,ans=1,cnt,n; while(scanf("%d",&n),n) { ans=1; for(i=0;i<n;i++) scanf("%lf%lf",&p[i].x,&p[i].y); for(i=0;i<n;i++) { for(int j=i+1;j<n;j++) { if(distance(p[i].x,p[i].y,p[j].x,p[j].y)-2.0>eps) continue; else { cnt=0; Point temp=Get_center(p[i],p[j],1.0); for(int k=0;k<n;k++) if(distance(temp.x,temp.y,p[k].x,p[k].y)-1.0<eps) cnt++; } if(cnt>ans) ans=cnt; } } printf("%d/n",ans); } return 0; }