enum = range(len(minT))

c1 = (rand.random_integers(0, 14000), rand.random_integers(0, 100))
c2 = (rand.random_integers(0, 14000), rand.random_integers(0, 100))
c3 = (rand.random_integers(0, 14000), rand.random_integers(0, 100))

c1p = []
c2p = []
c3p = []

for point in enum:
    dist = []
    x1 = point - c1[0]
    y1 = minT[point] - c1[1]
    dist1 = sqrt(x1**2 + y1**2)
    x2 = point - c2[0]
    y2 = minT[point] - c2[1]
    dist2 = sqrt(x2**2 + y2**2)
    x3 = point - c3[0]
    y3 = minT[point] - c3[1]
    dist3 = sqrt(x3**2 + y3**2)
    c1p.append((point, minT[point])) if dist1 < dist2 and dist1 < dist3 else c2p.append((point, minT[point])) if dist2 < dist3 else c3p.append((point, minT[point]))

scatter(*zip(*c1p), marker='x', c='red')
scatter(*zip(*c2p), marker='x', c='blue')
scatter(*zip(*c3p), marker='x', c='green')
scatter(*c1, marker='o', c='red', s=100)
scatter(*c2, marker='o', c='blue', s=100)
scatter(*c3, marker='o', c='green', s=100)
show()