BZOJ 1336&1337最小圆覆盖
来源:互联网 发布:知乎图片大小 编辑:程序博客网 时间:2024/05/17 01:51
思路:
http://blog.csdn.net/commonc/article/details/52291822
(照着算法步骤写……)
已知三点共圆 求圆心的时候 就设一下圆心坐标(x,y) 解个方程就好了
//By SiriusRen#include <cmath>#include <cstdio>#include <algorithm>using namespace std;int n;double R,tempx,tempy,tempz,tmpx,tmpy,tmpz;struct Point{double x,y;}point[100050],Ans;double Sqr(double x){return x*x;}double dis(Point a,Point b){return sqrt(Sqr(a.x-b.x)+Sqr(a.y-b.y));}bool in_circle(Point x){return dis(Ans,x)<=R;}int main(){ scanf("%d",&n); for(int i=1;i<=n;i++)scanf("%lf%lf",&point[i].x,&point[i].y); random_shuffle(point+1,point+n); for(int i=1;i<=n;i++)if(!in_circle(point[i])){ Ans.x=point[i].x,Ans.y=point[i].y,R=0; for(int j=1;j<i;j++)if(!in_circle(point[j])){ Ans.x=(point[i].x+point[j].x)/2; Ans.y=(point[i].y+point[j].y)/2; R=dis(Ans,point[j]); for(int k=1;k<j;k++)if(!in_circle(point[k])){ tempz=point[j].x-point[i].x; tempx=2*(point[i].y-point[j].y)/tempz; tempy=(Sqr(point[j].x)+Sqr(point[j].y)-Sqr(point[i].x)-Sqr(point[i].y))/tempz; tmpz=point[k].x-point[j].x; tmpx=2*(point[j].y-point[k].y)/tmpz; tmpy=(Sqr(point[k].x)+Sqr(point[k].y)-Sqr(point[j].x)-Sqr(point[j].y))/tmpz; Ans.y=(tmpy-tempy)/(tempx-tmpx); Ans.x=(tempx*Ans.y+tempy)/2; R=dis(Ans,point[j]); } } } printf("%f\n%f %f\n",R,Ans.x,Ans.y);}
0 0
- BZOJ 1336 & 1337 最小圆覆盖
- BZOJ 1336&1337最小圆覆盖
- BZOJ 1337 最小圆覆盖
- bzoj 1337: 最小圆覆盖 (最小圆覆盖)
- bzoj 1336 && bzoj 1337 最小圆覆盖 随机增量法
- 【BZOJ 1336&1337】 [Balkan2002]Alien最小圆覆盖
- BZOJ 1336 Balkan2002 Alien最小圆覆盖
- BZOJ 1336 [Balkan2002]Alien最小圆覆盖
- BZOJ 1336([Balkan2002]Alien最小圆覆盖-最小覆盖圆)
- bzoj 1336/bzoj 1337: [Balkan2002]Alien最小圆覆盖 随机增量法
- BZOJ 1336 Balkan2002 Alien最小圆覆盖 随机增量法
- BZOJ 1336: [Balkan2002]Alien最小圆覆盖 随机增量法
- BZOJ 1336: [Balkan2002]Alien最小圆覆盖 随机增量法
- BZOJ 1337 最小圆覆盖 随机增量法
- BZOJ 1337 最小圆覆盖 随机增量法
- [BZOJ 1336/1337][Balkan2002]Alien最小圆覆盖:随机增量法
- [随机增量法 最小圆覆盖] BZOJ 1366 [Balkan2002]Alien最小圆覆盖 & 1337 最小圆覆盖 & 2823 [AHOI2012]信号塔
- bzoj1336: [Balkan2002]Alien最小圆覆盖&&1337: 最小圆覆盖
- 建立docker私有仓库
- 李炎恢老师PHP第三季视频课程(设计模式+MVC模式+SMARTY+在线商城)
- 启动oracle的三过程
- 《DSAA》 12.2.3 红黑树的自顶向下删除
- Object-C之Bug集中营
- BZOJ 1336&1337最小圆覆盖
- <Principles of fMRI 1>课程笔记5--理解信号、噪音与血流
- 关于java String类的getBytes(String charsetName)和String(byte[] bytes, String charsetName)
- 智能楼控——新风系统控制思路
- 死磕微信!支付宝版“小程序”曝光
- java转行,你应该该做哪些工作
- go 语言 mongodb操作
- C++ 用libcurl库进行http通讯网络编程
- Java Streams,第 1 部分: java.util.stream 库简介