两球体积并

求两个球体积的并,分三种情况,两球有交集,两球相离,包含情况。

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

int main()
{
	double pi=acos(-1.0);
	double x1,y1,z1,r1,x2,y2,z2,r2;
	cin>>x1>>y1>>z1>>r1;
	cin>>x2>>y2>>z2>>r2;
	double d=sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)+(z1-z2)*(z1-z2));
	double v1=4.0/3.0*pi*r1*r1*r1;
	double v2=4.0/3.0*pi*r2*r2*r2;
	if(d>=r1+r2)
	{
		printf("%.12lf\n",v1+v2);
	}
	else if(fabs(r1-r2)>=d)
	{
		printf("%.12lf\n",max(v1,v2));
	}
	else
	{
		if(r1<r2)swap(r1,r2);
		double in=pi*(r1+r2-d)*(r1+r2-d)*(d*d+2*d*r2-3*r2*r2+2*d*r1+6*r2*r1-3*r1*r1)/(12.0*d);
		printf("%.12lf\n",v1+v2-in);
	}
	return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值