import numpy as np
N = int(input())
y = list(map(int, input().split()))
x = [i for i in range(N)]
var_x = np.var(x)
mu_x = np.mean(x)
mu_y = np.mean(y)
cov = sum((x-mu_x)*(y-mu_y)) / N
a = cov / var_x
b = mu_y - a*mu_x
cost = 0
for i in range(N):
    cost += (y[i] - a*x[i] - b) ** 2
print(b, a)
print(cost)