#include<bits/stdc++.h>
using namespace std;

int main()
{
    int i,j,t;
    double r1,r2,h,p,r ,ans,v,v1,x,v2,pi=acos(-1);
    cin>>t;
    for(i=1;i<=t;i++)
    {
        cin>>r1>>r2>>h>>p;
        r=((p*(r1-r2))/h)+r2;
        x=(h*r2)/(r1-r2);
        v1=(1/3.0)*pi*(r2*r2)*x;
        v2=(1/3.0)*pi*(r*r)*(p+x);
        v=v2-v1;
        printf("Case %d: %.10lf\n",i,v);
    }
    return 0;
}