【UVALive 4642】Malfatti Circles（圆，二分）

题

给定三角形，求三个两两相切且与三角形的一条边相切的圆的半径。

二分一个半径，可以得出另外两个半径，需要推一推公式（太久了，我忘记了）

#include<cstdio>
#include<cmath>
#define eps (1e-8)
#define sqr(a) (a)*(a)
#define min(a,b) (a)>(b)?(b):(a)
#define dd double
struct point{dd x,y,v,a;//点x,y，v为角度，a为边长
}q[10];
dd r0,r1,r2,r3;
int read(){int f=0;for(int i=0;i<3;i++){scanf("%lf%lf",&q[i].x,&q[i].y);if(q[i].x||q[i].y) f=1;}return f;
}
dd make(dd r,dd a,dd h,dd n){dd t=r-r/tan(h)/tan(n)+a/tan(n);if(t<=eps)return -1;return sqr((sqrt(t)-sqrt(r))*tan(n));
}
void solve(){dd l=0,r=min(q[0].a*tan(q[0].v),q[2].a*tan(q[2].v));while(r-l>eps){r0=(l+r)/2;r1=make(r0,q[0].a,q[0].v,q[1].v);r2=make(r0,q[2].a,q[0].v,q[2].v);if(r1<=eps||r2<=eps||r2/tan(q[2].v)+r1/tan(q[1].v)+2*sqrt(r2*r1)-q[1].a<eps)r=r0;elsel=r0;}
}
int main(){while(read()){for(int i=0;i<3;i++)q[i].a=sqrt(sqr(q[(i+1)%3].x-q[i].x)+sqr(q[(i+1)%3].y-q[i].y));//计算边长for(int i=0;i<3;i++)q[i].v=acos((sqr(q[i].a)+sqr(q[(i+2)%3].a)-sqr(q[(i+1)%3].a))/2/q[i].a/q[(i+2)%3].a)/2;solve();printf("%lf %lf %lf\n",r0,r1,r2);}
}

【UVALive 4642】Malfatti Circles（圆，二分）相关推荐

[51Nod](1278)相离的圆 ---- 二分查找
平面上有N个圆,他们的圆心都在X轴上,给出所有圆的圆心和半径,求有多少对圆是相离的. 例如:4个圆分别位于1, 2, 3, 4的位置,半径分别为1, 1, 2, 1,那么{1, 2}, {1, 3} ...
51NOD 1278 相离的圆(二分 + 排序好题)
平面上有N个圆,他们的圆心都在X轴上,给出所有圆的圆心和半径,求有多少对圆是相离的. 例如:4个圆分别位于1, 2, 3, 4的位置,半径分别为1, 1, 2, 1,那么{1, 2}, {1, 3} ...
51nod 1278 相离的圆二分+排序
1278 相离的圆基准时间限制:1 秒空间限制:131072 KB 分值: 10 难度:2级算法题平面上有N个圆,他们的圆心都在X轴上,给出所有圆的圆心和半径,求有多少对圆是相离的. 例如:4个 ...
Opencv——霍夫变换以及遇到的一些问题
目录问题1 :颜色空间转换函数参数问题:CV_BGR2GRAY vs CV_GRAY2BGR 问题2:cvRound().cvFloor().cvCeil()函数用法霍夫变换的含义标准霍夫直线变 ...
uni-app 引用地图
插件市场引用为 https://ext.dcloud.net.cn/plugin?id=615 插件引用后进行了部分修改目录 index页面组件页面 [创建nvue文件] 生成高德地图Key i ...
uniapp中进行地图定位
最近在用uniapp接入地图部分,之前用第三方的高德地图接入,内置浏览器正常显示,但是真机调试及打包测试是发现地图页空白 1.先获取appkey 获取高德地图KEY 2.引入map <map : ...
Opencv C++图像处理（全）
文章目录 Opencv官方资料一.入门基础 1.1.头文件说明:#include <opencv2/opencv.hpp> 1.2.头文件说明:#include <opencv2/ ...
数学术语的英汉对照(权威,全面)
abbreviation 简写符号:简写 abscissa 横坐标 absolute complement 绝对补集 absolute error 绝对误差 absolute inequality 绝 ...
数学英语词汇 (1)
abbreviation 简写符号:简写 abscissa 横坐标 absolute complement 绝对补集 absolute error 绝对误差 absolute inequality 绝 ...

【UVALive 4642】Malfatti Circles（圆，二分）

【UVALive 4642】Malfatti Circles（圆，二分）相关推荐

最新文章

热门文章