結果
問題 | No.132 点と平面との距離 |
ユーザー |
|
提出日時 | 2024-08-18 22:20:09 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 3,558 ms / 5,000 ms |
コード長 | 1,622 bytes |
コンパイル時間 | 267 ms |
コンパイル使用メモリ | 82,308 KB |
実行使用メモリ | 78,712 KB |
最終ジャッジ日時 | 2024-08-18 22:20:17 |
合計ジャッジ時間 | 6,249 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 3 |
ソースコード
# https://yukicoder.me/problems/no/132import mathdef inner_prod(v1, v2):ans =0for l in range(3):ans += v1[l] * v2[l]return ansdef norm2(v1):ans =0for l in range(3):ans += v1[l] * v1[l]return ansdef main():N = int(input())P = tuple(map(float, input().split()))xyz = []for _ in range(N):p = tuple(map(float, input().split()))xyz.append(p)answer = 0for i in range(N):p0 = xyz[i]v3 = tuple([P[l] - p0[l] for l in range(3)])for j in range(i + 1, N):p1 = xyz[j]v1 = tuple([p1[l] - p0[l] for l in range(3)])for k in range(j + 1, N):p2 = xyz[k]v2 = tuple([p2[l] - p0[l] for l in range(3)])# 係数の計算v1v2 = inner_prod(v1, v2)v3v1 = inner_prod(v3, v1)v3v2 = inner_prod(v3, v2)v1_n2 = norm2(v1)v2_n2 = norm2(v2)v3_n2 = norm2(v3)matrix = [[v2_n2, -v1v2],[-v1v2, v1_n2]]vector = [v3v1, v3v2]new_vector = [0, 0]for i in range(2):for j in range(2):new_vector[i] += vector[j] * matrix[i][j]d = v2_n2 * v1_n2 - v1v2 ** 2ans = v3_n2 - (new_vector[0] * v3v1 / d) - (new_vector[1] * v3v2 / d)answer += math.sqrt(ans)print(answer)if __name__ == "__main__":main()