結果

問題 No.3672 Volume 3D
コンテスト
ユーザー 👑 cologne
提出日時 2026-09-15 16:27:25
言語 PyPy3
(7.3.23 + ACL)
コンパイル:
pypy3 -mpy_compile _filename_
実行:
pypy3 _filename_
結果
AC  
実行時間 343 ms / 2,000 ms
+ 518µs
コード長 1,487 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 87 ms
コンパイル使用メモリ 82,900 KB
実行使用メモリ 88,512 KB
最終ジャッジ日時 2026-09-15 16:27:31
合計ジャッジ時間 5,520 ms
ジャッジサーバーID
（参考情報） 		judge2_0 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

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

        def cost_m_1(d2, ra, rb):

            D = d2 + ra ** 2 - rb ** 2
            bs = 4 * d2 * ra ** 2

            c2 = (D ** 2) / bs
            if D < 0:
                # same sign
                return -c2 ** 0.5 - 1
            else:
                # cosx - 1 = (cos^2x-1) / (cosx + 1)
                return (c2 - 1) / (c2 ** 0.5 + 1)

        da = cost_m_1(d2, ra, rb)
        db = cost_m_1(d2, rb, ra)

        res = ra ** 3 * da ** 2 * (da + 3) + rb ** 3 * db ** 2 * (db + 3)

        print(float(res) * pi / 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