#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;
}