結果

問題 No.3672 Volume 3D
コンテスト
ユーザー 👑 cologne
提出日時 2026-09-15 16:00:56
言語 PyPy3
(7.3.23 + ACL)
コンパイル:
pypy3 -mpy_compile _filename_
実行:
pypy3 _filename_
結果
WA  
実行時間 -
コード長 1,187 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 70 ms
コンパイル使用メモリ 81,280 KB
実行使用メモリ 87,936 KB
最終ジャッジ日時 2026-09-15 16:01:00
合計ジャッジ時間 3,885 ms
ジャッジサーバーID
(参考情報)
judge1_0 / judge2_0
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 1
other AC * 11 WA * 9
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code

import sys
from math import pi

from typing import List, Tuple

int1 = lambda x: int(x) - 1
input = lambda: sys.stdin.readline().rstrip('\n')
ii = lambda: int(input())
vi = lambda: list(map(int, input().split()))
vi1 = lambda: list(map(int1, input().split()))


def dbg(*args, **kwargs):
    print(*(repr(arg) for arg in args), *(f'{k}: {repr(v)}' for k, v in kwargs.items()),
          sep='; ', file=sys.stderr, flush=True)


def main():
    t = ii()
    for _ in range(t):
        xa, ya, za, ra = vi()
        xb, yb, zb, rb = vi()
        d2 = (xa - xb) ** 2 + (ya - yb) ** 2 + (za - zb) ** 2
        if d2 <= (ra - rb) ** 2:
            print(4 / 3 * min(ra, rb) ** 3 * pi)
            continue
        if d2 >= (ra + rb) ** 2:
            print(0)
            continue
        ca = (d2 + ra ** 2 - rb ** 2) / (2 * d2 ** 0.5 * ra)
        cb = (d2 + rb ** 2 - ra ** 2) / (2 * d2 ** 0.5 * rb)

        print(ra ** 3 * pi * (2 / 3 + (-ca + ca ** 3 / 3)) + rb ** 3 * pi * (2 / 3 + (-cb + cb ** 3 / 3)))


def _start():
    if (ret := main()) is not None:
        print(*ret) if isinstance(ret, List) or isinstance(ret, Tuple) else print(ret)


if __name__ == '__main__':
    _start()
0