結果

問題 No.132 点と平面との距離
コンテスト
ユーザー aru aru
提出日時 2020-08-14 15:58:09
言語 Go
(1.23.4)
結果
WA  
実行時間 -
コード長 2,066 bytes
コンパイル時間 11,146 ms
コンパイル使用メモリ 230,472 KB
実行使用メモリ 8,216 KB
最終ジャッジ日時 2024-10-10 12:00:55
合計ジャッジ時間 21,552 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other WA * 1 OLE * 2
権限があれば一括ダウンロードができます

ソースコード

diff # 

package main

import (
	"bufio"
	"fmt"
	"math"
	"os"
	"sort"
	"strconv"
)

func out(x ...interface{}) {
	fmt.Println(x...)
}

var sc = bufio.NewScanner(os.Stdin)

func getInt() int {
	sc.Scan()
	i, e := strconv.Atoi(sc.Text())
	if e != nil {
		panic(e)
	}
	return i
}

func getInts(N int) []int {
	ret := make([]int, N)
	for i := 0; i < N; i++ {
		ret[i] = getInt()
	}
	return ret
}

func getString() string {
	sc.Scan()
	return sc.Text()
}

// min, max, asub, absなど基本関数
func max(a, b int) int {
	if a > b {
		return a
	}
	return b
}

func min(a, b int) int {
	if a < b {
		return a
	}
	return b
}

func asub(a, b int) int {
	if a > b {
		return a - b
	}
	return b - a
}

func abs(a int) int {
	if a >= 0 {
		return a
	}
	return -a
}

func lowerBound(a []int, x int) int {
	idx := sort.Search(len(a), func(i int) bool {
		return a[i] >= x
	})
	return idx
}

func upperBound(a []int, x int) int {
	idx := sort.Search(len(a), func(i int) bool {
		return a[i] > x
	})
	return idx
}

func calcL(x0, y0, z0, x1, y1, z1, x2, y2, z2, x3, y3, z3 float64) float64 {
	ax := x2 - x1
	ay := y2 - y1
	az := z2 - z1
	bx := x3 - x1
	by := y3 - y1
	bz := z3 - z1

	// out(x1, y1, z1)
	// out(x2, y2, z2)
	// out(x3, y3, z3)
	// out(ax, ay, az, bx, by, bz)
	a := ay*bz - az*by
	b := az*bx - ax*bz
	c := ax*by - ay*bx
	d := -(a*x1 + b*y1 + c*z1)

	// out(a, b, c, d)

	ret := math.Abs(a*x0 + b*y0 + c*z0 + d)
	ret /= math.Sqrt(a*a + b*b + c*c)

	return ret
}

func getf() float64 {
	sc.Scan()
	f, _ := strconv.ParseFloat(sc.Text(), 64)
	return f
}

func main() {
	sc.Split(bufio.ScanWords)
	sc.Buffer([]byte{}, 1000000)
	N := getInt()
	x, y, z := getf(), getf(), getf()
	f := make([][3]float64, N)
	for i := 0; i < N; i++ {
		f[i][0], f[i][1], f[i][2] = getf(), getf(), getf()
	}

	ans := 0.0

	for i := 0; i < N; i++ {
		for j := i + 1; j < N; j++ {
			for k := j + 1; k < N; k++ {
				ret := calcL(x, y, z,
					f[i][0], f[i][1], f[i][2],
					f[j][0], f[j][1], f[j][2],
					f[k][0], f[k][1], f[k][2],
				)
				out(i, j, k, ret)
				ans += ret
			}
		}
	}
	out(ans)
}
0