HDU 4173（计算几何

题意：给出若干个点，求半径为2.5的圆最多包含多少个点。

思路：查到的一篇题解是枚举三个点的外接圆，这个对于我等计算几何小白来说还是不太好做。。。。后来听大神讲的，因为半径已知，所以只需要枚举经过两点的半径为2.5的所有圆就好（两个点对应两个圆），这样不仅好写，而且复杂度也降低了。

（感觉不用一大串模板都不好意思叫计算几何了。。。。。）

#include<iostream>
#include<map>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<vector>
#include<queue>
#include<stack>
#include<functional>
#include<set>
#include<cmath>
#define pb push_back
#define fs first
#define se second
#define sq(x) (x)*(x)
#define eps 0.0000000001
using namespace std;
typedef long long ll;
typedef pair<ll,ll> P;
const int maxv=250;
double x[maxv],y[maxv];
int n;
double cycx[maxv*maxv+30];
double cycy[maxv*maxv+30];
int cych=0;
double dis(double x1,double y1,double x2,double y2){
    return sqrt(sq(x1-x2)+sq(y1-y2));
}
void getcenter(double x1,double y1,double x2,double y2){
    double mx=(x1+x2)/2,my=(y1+y2)/2;
    double k=(y1-y2)/(x1-x2);
    k=-1/k;
    if(x1-x2==0) k=0;
    double dis1=dis(x1,y1,x2,y2);
    dis1/=2;
    double dis2=sqrt(sq(2.5)-sq(dis1));
    double dx1=1./sqrt(1+k*k)*dis2,dy1=k/sqrt(1+k*k)*dis2;
    double dx2=-dx1,dy2=-dy1;
    cycx[cych++]=mx+dx1,cycy[cych-1]=my+dy1;
    cycx[cych++]=mx+dx2,cycy[cych-1]=my+dy2;
}
int main(){
    freopen("/home/files/CppFiles/in","r",stdin);
/*    std::ios::sync_with_stdio(false);
    std::cin.tie(0);*/
    while(cin>>n){
        for(int i=0;i<n;i++) scanf("%lf%lf",x+i,y+i);
        cych=0;
        for(int i=0;i<n;i++){
            for(int j=i+1;j<n;j++){
                if(dis(x[i],y[i],x[j],y[j])<5+eps);
                getcenter(x[i],y[i],x[j],y[j]);
            }
        }
        int ans=1;
        for(int i=0;i<cych;i++){
            int cont=0;
            for(int j=0;j<n;j++){
                if(dis(x[j],y[j],cycx[i],cycy[i])<2.5+eps){
                    cont++;
                }
            }
            ans=max(ans,cont);
        }
        cout<<ans<<endl;
    }
    return 0;
}