计算几何(三条高线长求三角形面积) - Height to Area - UVA 10522

计算几何(三条高线长求三角形面积) - Height to Area - UVA 10522

题意：

$给定三角形的三条高的长度，计算三角形的面积。$

输入：

$多组测试数据，$

$首行包括一个正整数T，表示不合法数据的组数，$

$当不合法输入的数量达到T时，停止输入。$

$每组包括三个浮点数ha,hb,hc，分别表示三条边a,b,c上对应的高线的长度。$

输出：

$若有解，输出三角形的面积。小数点后保留三位。$

$否则输出：$These are invalid inputs!

Sample Input

1
31.573 22.352 63.448
46.300 50.868 86.683
22.005 24.725 22.914
5.710 25.635 32.805

Sample Output

1517.456
2219.941
311.804
These are invalid inputs!

---

分析：

结论：

$设三角形的三条高线长为:ha,hb,hc，则三角形的面积为:$

$$S_{△ABC}=\frac{1}{\sqrt{(\frac{1}{ha}+\frac{1}{hb}+\frac{1}{hc})(\frac{1}{ha}+\frac{1}{hb}-\frac{1}{hc})(\frac{1}{ha}-\frac{1}{hb}+\frac{1}{hc})(-\frac{1}{ha}+\frac{1}{hb}+\frac{1}{hc})}}$$

简要证明：

$计算三角形面积，我们知道已知三条边长a,b,c时，由海伦公式可求面积S。$

$本题中，已知高线长，我们能够进一步得到边长a=\frac{S}{ha},b=\frac{S}{hb},c=\frac{S}{hc}$

$根据海伦公式：$

$$S=\sqrt{p(p-a)(p-b)(p-c)}，p=\frac{a+b+c}{2}，得$$

$$S=S^2\sqrt{(\frac{1}{ha}+\frac{1}{hb}+\frac{1}{hc})(\frac{1}{ha}+\frac{1}{hb}-\frac{1}{hc})(\frac{1}{ha}-\frac{1}{hb}+\frac{1}{hc})(-\frac{1}{ha}+\frac{1}{hb}+\frac{1}{hc})}$$

$即：$

$$S=\frac{1}{\sqrt{(\frac{1}{ha}+\frac{1}{hb}+\frac{1}{hc})(\frac{1}{ha}+\frac{1}{hb}-\frac{1}{hc})(\frac{1}{ha}-\frac{1}{hb}+\frac{1}{hc})(-\frac{1}{ha}+\frac{1}{hb}+\frac{1}{hc})}}$$

无解情况：

$不满足两边之和大于第三边，即：$

$$b+c\le a或a+b\le c或a+c\le b$$

$对应到高上为：$

$$\frac{S}{hb}+\frac{S}{hc}\le \frac{S}{ha}$$

$由于S>0，我们可以把S约掉$

即：

$$\frac{1}{hb}+\frac{1}{hc}\le \frac{1}{ha}或\frac{1}{ha}+\frac{1}{hb}\le \frac{1}{hc}或\frac{1}{ha}+\frac{1}{hc}\le \frac{1}{hb}$$

代码：

#define _CRT_SECURE_NO_WARNINGS

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

using namespace std;

const double eps = 1e-10;
const double pi = acos(-1.0);

int dcmp(double x)
{
    if (fabs(x) < eps) return 0;
    else return x < 0 ? -1 : 1;
}

int T;
double ha, hb, hc;

double InverseSum(double a, double b, double c)
{
    return 1.0 / a + 1.0 / b + 1.0 / c;
}

double area(double ha, double hb, double hc)
{
    if (dcmp(ha) <= 0 || dcmp(hb) <= 0 || dcmp(hc) <= 0) return -1;
    if (1.0 / ha >= 1.0 / hb + 1.0 / hc 
        || 1.0 / hb >= 1.0 / ha + 1.0 / hc 
        || 1.0 / hc >= 1.0 / ha + 1.0 / hb) return -1;
    double s1 = InverseSum(ha, hb, hc), s2 = InverseSum(ha, hb, -hc),
        s3 = InverseSum(ha, -hb, hc), s4 = InverseSum(-ha, hb, hc);
    double S = s1*s2*s3*s4;
    return sqrt(1.0/S);
}

int main()
{
    scanf("%d", &T);
    while (T)
    {
        scanf("%lf%lf%lf", &ha, &hb, &hc);
        double res = area(ha, hb, hc);
        if (res < 0)
        {
            T--;
            puts("These are invalid inputs!");
        }
        else printf("%.3lf\n", res);
    }
    return 0;
}