結果

問題 No.132 点と平面との距離
ユーザー LyricalMaestro
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

# https://yukicoder.me/problems/no/132
import math
def inner_prod(v1, v2):
ans =0
for l in range(3):
ans += v1[l] * v2[l]
return ans
def norm2(v1):
ans =0
for l in range(3):
ans += v1[l] * v1[l]
return ans
def 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 = 0
for 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 ** 2
ans = v3_n2 - (new_vector[0] * v3v1 / d) - (new_vector[1] * v3v2 / d)
answer += math.sqrt(ans)
print(answer)
if __name__ == "__main__":
main()
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0