Ural 1550 Dean's Pyramid 3

Sol：V-ans  = V-正四面体  -  V-截柱

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

using namespace std;

const double pi=acos(-1.0);

double h,w,x,y,r;

int main()
{
    scanf("%lf%lf%lf%lf%lf",&h,&w,&x,&y,&r);
    double H=max(fabs(x),fabs(y));
    double v1=h*w*w/3.0;
    double v2=pi*r*r*(h-H/w*h*2);
    printf("%.6f\n",v1-v2);
    return 0;
}