n = int(input()) h = float(input()) points = [tuple(map(int, input().split())) for _ in range(n)] a1 = 0.0 a2 = 0.0 prev_assigned_f = None while True: assigned_f = [] for x, y in points: residual_f = (y - (a1 * x + h)) ** 2 residual_g = (y - (a2 * x)) ** 2 assigned_f.append(residual_f < residual_g) if prev_assigned_f == assigned_f: break prev_assigned_f = assigned_f.copy() sum_xy_f = 0.0 sum_x2_f = 0.0 sum_xy_g = 0.0 sum_x2_g = 0.0 for i in range(n): x, y = points[i] if assigned_f[i]: sum_xy_f += x * (y - h) sum_x2_f += x * x else: sum_xy_g += x * y sum_x2_g += x * x new_a1 = sum_xy_f / sum_x2_f if sum_x2_f != 0 else 0.0 new_a2 = sum_xy_g / sum_x2_g if sum_x2_g != 0 else 0.0 a1, a2 = new_a1, new_a2 total = 0.0 for x, y in points: residual_f = (y - (a1 * x + h)) ** 2 residual_g = (y - (a2 * x)) ** 2 total += min(residual_f, residual_g) print("{0:.20f}".format(total))