【两球体积并】

2021-07-19 21:02:58 阅读：181 来源： 互联网

#include <iostream>
#include <cstring>
#include <iomanip>
#include <algorithm>
#include <queue>
#include <stack>
#include <vector>
#include <set>
#include <map>
#include <cmath>
#include <cstdio>
using namespace std;
#define inf 1e18
typedef long long ll;
typedef long double ld;
const ll maxn = 1e5+100;
const ll mod = 1e9+7;
const ld pi = acos(-1.0);

ld xa,ya,za,ra,xb,yb,zb,rb;
ld d,va,vb,ans,ha,hb;

int main() {
	ios::sync_with_stdio(false);

	cin >> xa >> ya >> za >> ra >> xb >> yb >> zb >> rb;

	d = sqrt( (xa-xb)*(xa-xb) + (ya-yb)*(ya-yb) + (za-zb)*(za-zb) ); //圆心距离
	va = (4.0/3.0)*pi*(ra*ra*ra);   // a球体积
	vb = (4.0/3.0)*pi*(rb*rb*rb);   // b球体积

	if(d >= ra+rb) {     //判断相离或相切
		ans = va+vb;
	} else if( d+min(ra,rb) <= max(ra,rb) ) { // 包含
		ans = max(va,vb);
	} else {   //相切
		ha = ra - ( ( (ra*ra) - (rb*rb) + (d*d) ) / (2*d) ); //球冠的高
		hb = rb - ( ( (rb*rb) - (ra*ra) + (d*d) ) / (2*d) );

		ld vra = ( pi*ha*ha*( 3*ra-ha  ) ) / 3.0;     //球冠体积
		ld vrb = ( pi*hb*hb*( 3*rb-hb  ) ) / 3.0;

		ans = va + vb - vra - vrb;
	}

	cout << fixed << setprecision(7) << ans << endl;

	return 0;
}

标签：xb,xa,两球,za,体积,ra,rb,include
来源： https://www.cnblogs.com/Tecode/p/15032183.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

【两球体积并】