結果

問題 No.3672 Volume 3D
コンテスト
ユーザー daiota
提出日時 2026-09-05 00:35:07
言語 C++17
(gcc 15.3.0 + boost 1.92.0 + ACL)
コンパイル:
g++-15 -O2 -lm -std=c++17 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
WA  
実行時間 -
コード長 974 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,231 ms
コンパイル使用メモリ 211,476 KB
実行使用メモリ 6,400 KB
最終ジャッジ日時 2026-09-05 00:35:13
合計ジャッジ時間 3,572 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge4_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 11 WA * 9
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef pair<ll,ll> P;
#define REP(i,n) for(ll i=0;i<ll(n);i++)




const double pi=acos(-1.0);





int main(void){
	cin.tie(nullptr);  ios_base::sync_with_stdio(false);
	ll i,j;




	int T;
	cin >> T;
	while(T--){

		ll a,b,c,r;
		cin >> a >> b >> c >> r;

		ll a1,b1,c1,r1;
		cin >> a1 >> b1 >> c1 >> r1;


		ll d=(a-a1)*(a-a1)+(b-b1)*(b-b1)+(c-c1)*(c-c1);

		if((r+r1)*(r+r1)<=d){
			cout << 0 << endl;
		}

		else if(d<=(r-r1)*(r-r1)){

			if(r<=r1) cout << fixed << setprecision(12) << (double)(4.0/3.0*pi*r*r*r) << endl;
			else cout << fixed << setprecision(12) << (double)(4.0/3.0*pi*r1*r1*r1) << endl;


		}
		else{

			double x=(r*r-r1*r1+d)/(2.0*sqrt(1.0*d));

			double v=pi*(2.0*r*r*r/3.0-1.0*r*r*x+x*x*x/3.0);
			double q=sqrt(1.0*d)-x;
			double u=pi/3.0*r1*r1*r1*2.0-pi*1.0*r1*r1*q+pi*q/3.0*q*q;

			cout << fixed << setprecision(12) << v+u << endl;

		}









	}






	return 0;

}
0