https://ac.nowcoder.com/acm/contest/373/E
题解:
参考文章:https://blog.csdn.net/luyehao1/article/details/86583384
/*
*@Author: STZG
*@Language: C++
*/
#include <bits/stdc++.h>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<string>
#include<vector>
#include<bitset>
#include<queue>
#include<deque>
#include<stack>
#include<cmath>
#include<list>
#include<map>
#include<set>
//#define DEBUG
#define RI register int
using namespace std;
typedef long long ll;
//typedef __int128 lll;
const int N=100000+10;
const int MOD=1e9+7;
const double PI = acos(-1.0);
const double EXP = 1E-8;
const int INF = 0x3f3f3f3f;
double t,n,m,k,q,x,y,z,r,x2,y2,z2,r2;
double ans,cnt,flag,temp;
int a[N];
char str;
int main()
{
#ifdef DEBUG
freopen("input.in", "r", stdin);
//freopen("output.out", "w", stdout);
#endif
scanf("%lf%lf%lf%lf%lf%lf%lf%lf",&x,&y,&z,&r,&x2,&y2,&z2,&r2);
double d=sqrt((x-x2)*(x-x2)+(y-y2)*(y-y2)+(z-z2)*(z-z2));
if(r>r2){
swap(r,r2);
swap(x,x2);
swap(y,y2);
swap(z,z2);
}
if (d >= r + r2) {
ans=(4.0 / 3)*PI*(r*r*r+r2*r2*r2);
}else if (d + r <= r2) {
ans = (4.0 / 3)*PI*r2*r2*r2;
}else {
double co = (r2*r2 + d * d - r*r) / (2.0*d*r2);
double h = r2*(1 - co);
ans += (1.0 / 3)*PI*(3.0*r2 - h)*h*h;
co = (r*r + d * d - r2*r2) / (2.0*d*r);
h = r*(1 - co);
ans += (1.0 / 3)*PI*(3.0*r - h)*h*h;
ans=(4.0 / 3)*PI*(r*r*r+r2*r2*r2)-ans;
}
printf("%.7lf\n", ans);
//cout << "Hello world!" << endl;
return 0;
}