求一个覆盖所有点的最小圆 最小圆覆盖

 1 float disy(float y,float x)
 2 {
 3     int n=pvec.size();
 4     float Rmin=0;
 5     
 6     for(int i=0;i<n;i++)
 7     {
 8         Rmin=max(Rmin,sqrt((pvec[i].y-y)*(pvec[i].y-y)+ (pvec[i].x-x)*(pvec[i].x-x)));
 9     }
10     return Rmin;
11 }
12 float dis(float x)
13 {
14     ymax=10000,ymin=0;
15     while(ymax-ymin>accu)
16     {
17         if(disy(ymin+(ymax-ymin)*2.0/3,x)<disy(ymin+(ymax-ymin)*1.0/3,x))
18         {
19             ymin=ymin+(ymax-ymin)*1.0/3;
20         }
21         else
22             ymax=ymin+(ymax-ymin)*2.0/3;
23     }
24     return disy((ymax+ymin)/2,x);
25 }
26 int findcenter(POINT &pstart)
27 {
28     xmax=0,xmin=10000;
29     int n=pvec.size();    
30     for(int i=0;i<n;i++)
31     {
32         xmax=max(xmax,pvec[i].x);
33         xmin=min(xmin,pvec[i].x);
34     }
35     while((xmax-xmin)/2>accu)
36     {
37         if(dis(xmin+(xmax-xmin)*1.0/3)<dis(xmin+(xmax-xmin)*2.0/3))
38         {
39             xmax=xmin+(xmax-xmin)*2.0/3;
40         }
41         else
42             xmin=xmin+(xmax-xmin)*1.0/3;
43     }
44     pstart.x=(xmax+xmin)/2;
45     pstart.y=(ymax+ymin)/2;
46     int r=disy((ymax+ymin)/2,(xmax+xmin)/2);
47     return r;
48 }