Ex-circles - UVa 11731 几何

最新推荐文章于 2019-08-19 22:30:00 发布

提比-我有特殊的AC技巧

最新推荐文章于 2019-08-19 22:30:00 发布

阅读量699

点赞数

分类专栏：几何 UVa

本文链接：https://blog.csdn.net/u014733623/article/details/43527599

版权

UVa 同时被 2 个专栏收录

218 篇文章 0 订阅

订阅专栏

几何

39 篇文章 0 订阅

订阅专栏

题意：给你三角形的三条边，求图中DEF的面积和阴影部分的面积。

思路：这是一个旁切圆的问题，三角形关于顶点A、B、C的旁切圆的半径分别是 $\frac{2S}{-a+b+c}$ 、 $\frac{2S}{a-b+c}$ 和 $\frac{2S}{a+b-c}$ ，其中 $S$ 表示三角形面积，a、b、c分别是A、B、C的对边。

另外补充知识点：

若顶点的坐标分别为 $(x_1,y_1)$ 、 $(x_2,y_2)$ 、 $(x_3,y_3)$ ，则三个旁心的座标为：

J_a=(\frac{-ax_1+bx_2+cx_3}{-a+b+c},\frac{-ay_1+by_2+cy_3}{-a+b+c}),J_b=(\frac{ax_1-bx_2+cx_3}{a-b+c},\frac{ay_1-by_2+cy_3}{a-b+c}),J_c=(\frac{ax_1+bx_2-cx_3}{a+b-c},\frac{ay_1+by_2-cy_3}{a+b-c})

AC代码如下：

#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
double angle(double a,double b,double c)
{
    return acos((a*a+b*b-c*c)/(2*a*b));
}
double area(double a,double b,double c)
{
    double p=(a+b+c)/2;
    return sqrt(p*(p-a)*(p-b)*(p-c));
}
double eps=1e-8;
int main()
{
    int t=0,i,j,k;
    double a,b,c,A,B,C,S,ra,rb,rc,ans1,ans2;
    while(~scanf("%lf%lf%lf",&a,&b,&c) && a+b+c>eps)
    {
        S=area(a,b,c);
        A=angle(b,c,a);
        B=angle(a,c,b);
        C=angle(a,b,c);
        ra=2*S/(b+c-a);
        rb=2*S/(a+c-b);
        rc=2*S/(a+b-c);
        ans1=S+(a*ra+b*rb+c*rc)/2;
        ans2=(M_PI/2-A/2)*ra*ra+(M_PI/2-B/2)*rb*rb+(M_PI/2-C/2)*rc*rc;
        ans2/=2;
        printf("Case %d: %.2f %.2f\n",++t,ans1,ans2);
    }
}