hdu--1077--Catching Fish

思路：
  1.枚举两点确定圆心，大于2不用考虑
  2.逐个判断判断距圆心的距离小于1.00001符合题意
这个题，主要在求圆心上废了不少功夫，但是仍存在问题

#include<iostream> 
#include<cmath>
#include<vector>
#include<cstdio>
const static double eps = 1e-6;
using namespace std;
struct point{
    double x,y;
    point():x(0),y(0){};
};
double getDistance(point,point);
point getCentral(point,point);
int main()
{    
    int t;
    cin>>t;
    while(t--){
        vector<point> vec;
        int n;cin>>n;
        for(int i=0;i<n;++i){
            point temp;
            cin>>temp.x>>temp.y;
            vec.push_back(temp);
        }
        int cut = 1;//捕鱼数 
        auto it=vec.begin();
        for(int i=0;i<vec.size()-1;++i,++it){
            auto iter=it;
            for(++iter;iter!=vec.end();++iter){
                point cen;//圆心
                if(getDistance((*it),(*iter)) >2.0)continue;
                cen = getCentral((*it),(*iter));
                int k = 0;
                for(auto p=vec.begin();p!=vec.end();++p){
                    if(getDistance(cen,(*p)) <= 1.0001)k++;
                }
                if(cut < k)cut=k;
            }
        }
        cout<<cut<<endl;
    }
    return 0;
}
double getDistance(point a,point b){
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
point getCentral(point a,point b){
    point o;
    point c;//a,b中点 
    c.x=(a.x+b.x)/2;
    c.y=(a.y+b.y)/2;
    double l=getDistance(c,o);//oc长度 
    double temp=(getDistance(a,b)/2);
    double len=sqrt(1.0-temp*temp);//圆心到直线AB的距离 
    point central;
//    if(fabs(a.y-b.y)<eps){
//        central.x=c.x;
//        central.y=c.y+len;
//    }
//    else {//大牛考虑了两个点几乎重合的情况，按照题意不必考虑 
        double ang=atan(-(a.x-b.x)/(a.y-b.y));
        central.x=c.x+cos(ang)*len;//cos*斜边长 
        central.y=c.y+sin(ang)*len;//sin*斜边长    
        /*
        最初自己用临边/斜边的方式模拟正余弦值，但是误差不小，
        甚至两点确定的圆，距离两点的距离大于1
        */
/*注意：这个圆心最后是用中点的坐标加上cos*斜边长，但是减去呢？明显不对了虽然A了，但枚举不是正确的方法*/
//     }
    return central;
}