結果

問題 No.132 点と平面との距離
コンテスト
ユーザー obakyan
提出日時 2019-05-10 20:24:54
言語 Lua
(LuaJit 2.1.1734355927)
結果
AC  
実行時間 2,835 ms / 5,000 ms
コード長 1,489 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 70 ms
コンパイル使用メモリ 6,816 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-02 01:01:43
合計ジャッジ時間 4,643 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 3
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

local ior = io.input()
local n = ior:read("*n")
local px, py, pz = ior:read("*n", "*n", "*n")
local tx, ty, tz = {}, {}, {}
for i = 1, n do
  tx[i], ty[i], tz[i] = ior:read("*n", "*n", "*n")
end

local function getlen(idxs)
  local x, y, z = {}, {}, {}
  for i = 1, 3 do
    x[i], y[i], z[i] = tx[idxs[i]], ty[idxs[i]], tz[idxs[i]]
  end
  x[4], y[4], z[4] = px, py, pz
  for i = 2, 4 do
    x[i], y[i], z[i] = x[i] - x[1], y[i] - y[1], z[i] - z[1]
  end
  local e1sz = math.sqrt(x[2] * x[2] + y[2] * y[2] + z[2] * z[2])
  local e1 = {x[2] / e1sz, y[2] / e1sz, z[2] / e1sz}

  local innerpr = e1[1] * x[3] + e1[2] * y[3] + e1[3] * z[3]
  local e2 = {
    x[3] - innerpr * e1[1],
    y[3] - innerpr * e1[2],
    z[3] - innerpr * e1[3]
    }
  local e2sz = math.sqrt(e2[1] * e2[1] + e2[2] * e2[2] + e2[3] * e2[3])
  for i = 1, 3 do e2[i] = e2[i] / e2sz end

  local e3 = {x[4], y[4], z[4]}
  innerpr = e3[1] * e1[1] + e3[2] * e1[2] + e3[3] * e1[3]
  for i = 1, 3 do
    e3[i] = e3[i] - e1[i] * innerpr
  end
  innerpr = e3[1] * e2[1] + e3[2] * e2[2] + e3[3] * e2[3]
  for i = 1, 3 do
    e3[i] = e3[i] - e2[i] * innerpr
  end
  local e3sz = math.sqrt(e3[1] * e3[1] + e3[2] * e3[2] + e3[3] * e3[3])
  for i = 1, 3 do e3[i] = e3[i] / e3sz end
  innerpr = e3[1] * x[4] + e3[2] * y[4] + e3[3] * z[4]
  return innerpr
end

local sum = 0
for i = 1, n - 2 do
  for j = i + 1, n - 1 do
    for k = j + 1, n do
      sum = sum + getlen({i, j, k})
    end
  end
end
print(string.format("%.10f", sum))
0