计算几何(三角形内切圆 + 海伦公式) - InCircle - UVA 11524

本文链接：https://blog.csdn.net/njuptACMcxk/article/details/108759780

计算几何(三角形内切圆 + 海伦公式) - InCircle - UVA 11524

题意：

$对于\triangle ABC，$

在这里插入图片描述
$假设其内切圆与其相切于 P 、 Q 、 R 三点，$

$给出比例：m_1:n_1，m_2:n_2，m_3:n_3，以及内切圆的半径，$

$求出\triangle ABC的面积。$

输入：

$T 组测试数据，$

$每组首行包括一个浮点数 R ，表示内切圆的半径。$

$接着三行，每行两个浮点数，m_i和n_i$

输出：

$一个浮点数，表示S_{\triangle ABC}，小数点后保留4位数字。$

Sample Input

2
140.9500536497
15.3010457320 550.3704847907
464.9681681852 65.9737378230
55.0132446384 10.7791711946
208.2835101182
145.7725891419 8.8264176452
7.6610997600 436.1911036207
483.6031801012 140.2797089713

Sample Output

400156.4075
908824.1322

数据范围：

$0 < T < 50001，1 < R < 5000，1 < m_1, n_1, m_2, n_2, m_3, n_3 < 50000$

分析：

在这里插入图片描述
$根据三角形内切圆的性质，有： A P = A R ， B P = B Q ， C Q = C R$

$设 A P = x = A R$

$则：BP=\frac{n_1}{m_1}x=BQ，CR=\frac{m_3}{n_3}x=CQ$

$\qquad AB=c=(1+\frac{n_1}{m_1})x，AC=b=(1+\frac{m_3}{n_3})x，BC=a=(\frac{n_1}{m_1}+\frac{m_3}{n_3})x=b+c-2x$

$一个未知数，可列一个方程来解出，$

$我们知道三条边 a, b, c ，以及内切圆的半径 R ，容易想到通过三角形面积来建立等量关系，$

$由海伦公式：$

$S_{\triangle{ABC}}=\sqrt{p(p-a)(p-b)(p-c)}，其中p=\frac{a+b+c}{2}$

$又：$

$S_{\triangle{ABC}}=\frac{(a+b+c)R}{2}=pR$

$联立得：$

$p-a)(p-b)(p-c)=pR^2$

$而：$

$p=\frac{(a+b+c)}{2}=b+c-1=(1+\frac{n_1}{m_1}+\frac{m_3}{n_3})x$

$带入解得：$

$x=\sqrt{\frac{(1+\frac{n_1}{m_1}+\frac{m_3}{n_3})R^2}{\frac{n_1m_3}{m_1n_3}}}$

$带入能够解得三角形面积。$

代码：

#include<cstring>
#include<iostream>
#include<cmath>
#include<cstdio>
#include<algorithm>

using namespace std;

int main()
{
    int T;
    double R;
    double n[4],m[4];
    scanf("%d",&T);
    while(T--)
    {
        scanf("%lf",&R);
        for(int i=1;i<=3;i++) scanf("%lf%lf",&m[i],&n[i]);
        double x=sqrt((1+n[1]/m[1]+m[3]/n[3])*R*R/(n[1]*m[3]/m[1]/n[3]));
        double S=(1+n[1]/m[1]+m[3]/n[3])*x*R;
        printf("%.4lf\n",S);
    }
    return 0;
}