求两个球的体积并

两个球和两个圆一样，有三种情况。相交，分离（包括两种情况）
对于分离的两种情况，好计算。但是对于相交的情况，不好计算。需要计算球冠的大小。

求两个球的体积并
就是这样。昨天晚上牛客上有一个关于这样的题
代码如下：

#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)/~

求两个球的体积并

相关推荐