Submission #1857242

---

Source Code Expand

import numpy as np
import math
for s in range(6):
    v.append([float(a) for a in raw_input().split()])
verts = np.array(v)
#print verts
v1 = verts[1] - verts[0]
v2 = verts[2] - verts[1]
v3 = verts[3] - verts[2]
v4 = verts[4] - verts[3]
v5 = verts[5] - verts[4]
v6 = verts[0] - verts[5]

n1 = np.array([v1[1], -v1[0], 0])
n2 = np.array([v2[1], -v2[0], 0])
n3 = np.array([v3[1], -v3[0], 0])
r12 = n1.dot(n2)
r23 = n2.dot(n3)
r13 = n1.dot(n3)
#print r12, r23, r13
n1[2] = math.sqrt(-r12 * r13 / r23)
n2[2] = math.sqrt(-r12 * r23 / r13)
n3[2] = math.sqrt(-r13 * r23 / r12)
n1 = n1 / np.linalg.norm(n1)
n2 = n2 / np.linalg.norm(n2)
n3 = n3 / np.linalg.norm(n3)
d1 =  np.abs(n1.dot(verts[0]) - n1.dot(verts[3]))
d2 =  np.abs(n2.dot(verts[1]) - n2.dot(verts[4]))
d3 =  np.abs(n3.dot(verts[2]) - n3.dot(verts[5]))
V1= d1 * d2 * d3

n1 = np.array([v1[1], -v1[0], 0])
n2 = np.array([v2[1], -v2[0], 0])
n3 = np.array([v3[1], -v3[0], 0])
r12 = n1.dot(n2)
r23 = n2.dot(n3)
r13 = n1.dot(n3)
#print r12, r23, r13
n1[2] = -math.sqrt(-r12 * r13 / r23)
n2[2] = -math.sqrt(-r12 * r23 / r13)
n3[2] = -math.sqrt(-r13 * r23 / r12)
n1 = n1 / np.linalg.norm(n1)
n2 = n2 / np.linalg.norm(n2)
n3 = n3 / np.linalg.norm(n3)
d1 =  np.abs(n1.dot(verts[0]) - n1.dot(verts[3]))
d2 =  np.abs(n2.dot(verts[1]) - n2.dot(verts[4]))
d3 =  np.abs(n3.dot(verts[2]) - n3.dot(verts[5]))
V2= d1 * d2 * d3
print min(V1, V2)

Submission Info

Judge Result