import numpy as np n = int(input()) Y = [] meany = 0 for x in input().split(): x = int(x) meany += x Y.append(x) meany /= n meanx = (n - 1) / 2 Y = np.array(Y) - meany X = np.array(range(n)) - meanx # print(X) Sxy = 0 Sx2 = 0 for x, y in zip(X,Y): Sxy += x * y Sx2 += x ** 2 d = Sxy / Sx2 b = meany - d * meanx Y += meany ai = b cost = 0 for x in Y: cost += (x - ai) ** 2 ai += d print(b, d) print(cost)