两个球和两个圆一样,有三种情况。相交,分离(包括两种情况)
对于分离的两种情况,好计算。但是对于相交的情况,不好计算。需要计算球冠的大小。
就是这样。昨天晚上牛客上有一个关于这样的题
代码如下:
#include<bits/stdc++.h>
using namespace std;
const double PI = acos(-1.0);
double x1,y11,z1,r1;
double x2,y2,z2,r2;
//球冠的体积:PI*h*h/(3*r-h)其中h是圆的半径r-圆心到截圆的平面的距离dis。h1=r1-(r1*r1-r2*r2+d*d)/(2*d)其中d为两圆心之间的距离。
int main()
{
scanf("%lf%lf%lf%lf",&x1,&y11,&z1,&r1);
scanf("%lf%lf%lf%lf",&x2,&y2,&z2,&r2);
double dis=sqrt((x1-x2)*(x1-x2)+(y11-y2)*(y11-y2)+(z1-z2)*(z1-z2));
if(dis>=r1+r2)
{
double s=4.0/3.0*PI*(r1*r1*r1+r2*r2*r2);
printf("%.7lf\n",s);
}
else if(dis<=fabs(r1-r2))
{
double maxn=max(r1,r2);
double s=4.0/3.0*PI*(maxn*maxn*maxn);
printf("%.7lf\n",s);
}
else
{
double s=3.0/4.0*PI*(r1*r1*r1+r2*r2*r2);
double h1=r1-(r1*r1-r2*r2+dis*dis)/(2*dis);
s-=(PI*h1*h1*(3*r1-h1)/3.0);
double h2=r2-(r2*r2-r1*r1+dis*dis)/(2.0*dis);
s-=(PI*h2*h2*(3*r2-h2)/3.0);
printf("%.7lf\n",s);
}
}
还是把公式记住,要不比赛的时候自己推太麻烦了。
球冠的体积:PIhh/(3r-h)其中h是圆的半径r-圆心到截圆的平面的距离dis。h1=r1-(r1r1-r2r2+dd)/(2*d)其中d为两圆心之间的距离
努力加油a啊(o)/~