data Vector3D = Vector3D{x :: Double, y :: Double, z :: Double}deriving Show data Plane = Plane{norm :: Vector3D, point::Vector3D} deriving Show vcross :: Vector3D -> Vector3D -> Vector3D vcross (Vector3D ax ay az) (Vector3D bx by bz) = Vector3D (ay*bz-az*by) (az*bx-ax*bz) (ax*by-ay*bx) vadd :: Vector3D -> Vector3D -> Vector3D vadd (Vector3D ax ay az) (Vector3D bx by bz) = Vector3D (ax+bx) (ay+by) (az+bz) vinv :: Vector3D -> Vector3D vinv (Vector3D x y z) = Vector3D (-x) (-y) (-z) vsub :: Vector3D -> Vector3D -> Vector3D vsub a b = vadd a (vinv b) main = do getLine p <- return . (\(x:y:z:_) -> Vector3D x y z) . map read . words =<< getLine q <- return . map ((\(x:y:z:_) -> Vector3D x y z) . map read . words) . lines =<< getContents let pl = [Plane (vcross (vsub (q!!i) (q!!j)) (vsub (q!!k) (q!!j))) (q!!i) | i<-[0..length q-1], j<-[i+1..length q-1], k<-[j+1..length q-1]] print $ sum $ map (\(Plane (Vector3D a b c) (Vector3D d e f)) -> abs (a*(x p) + b*(y p) + c*(z p) - a*d - b*e - c*f) / sqrt(a*a+b*b+c*c) ) pl