HDU 1451（Area in Triangle）

#include <cstdio>
#include <cmath>
using namespace std;
const double eps = 1e-9;
const double PI = acos(-1.0);

int main()
{
    int Case = 1;
    double a, b, c, len;
    while (scanf("%lf%lf%lf%lf", &a, &b, &c, &len) != EOF)
    {
        if (a == 0 && b == 0 && c == 0 && len == 0)
            break;

        double p0 = a + b + c;
        double p1 = p0 / 2;
        double S = sqrt(p1 * (p1 - a) * (p1 - b) * (p1 - c));

        printf("Case %d: ", Case++);
        if (len >= p0)
        { 
            printf("%.2lf\n", S); 
            continue; 
        }

        double R = 2 * S / p0;
        if (2 * PI * R - len > eps)
        {
            R = len / PI / 2;
            S = PI * R * R;
            printf("%.2lf\n", S);
            continue;
        }
        double r = (p0 - len) / (p0 / R - 2 * PI);
        a = a / R * (R - r);
        b = b / R * (R - r);
        c = c / R * (R - r);
        double p = (a + b + c) / 2;
        S = PI * r * r + sqrt(p * (p - a) * (p - b) * (p - c)) + r * (a + b + c);
        printf("%.2lf\n", S);
    }
    return 0;
}