import numpy as np from math import pi, cos, sin, atan2 import sys ZERO = np.array([0.0, 0.0]) get = lambda: np.array(list(map(float, input().split()))) def norm(a) -> float: x = abs(a[0]) y = abs(a[1]) return x ** p + y ** p def dist(a, b): return norm(a - b) def cross(a, b): return a[0] * b[1] - a[1] * b[0] p = float(input()) a = get() b = get() c = get() if a[0] > b[0]: a, b = b, a if a[0] > c[0]: a, c = c, a if cross(b - a, c - a) < 0: b, c = c, b A = a.copy() a -= A b -= A c -= A B = -pi C = pi def atan(diff) -> float: return atan2(diff[1], diff[0]) def near(diff) -> int: da = dist(a, diff) db = dist(b, diff) dc = dist(c, diff) mn = min(da, db, dc) if da == mn: return 0 if db == mn: return 1 return 2 def bisect_a(di): ok = 1e7 ng = 0.0 for _ in range(100): mid = (ok + ng) / 2 if near(di * mid) == 0: ng = mid else: ok = mid return di * ok def direction(theta): return np.array([cos(theta), sin(theta)]) def check(theta): di = direction(theta) x = bisect_a(di) n = near(x) if n: return n if theta < atan(b): return 1 return 2 for _ in range(100): theta = (B + C) / 2 if check(theta) == 1: B = theta else: C = theta ans = bisect_a(direction(B)) if max(abs(ans[0]), abs(ans[1])) > 1e6: ans = bisect_a(direction(C)) print(dist(a, ans), dist(b, ans), dist(c, ans), file=sys.stderr) print(*(ans + A))