import numpy as np
n = int(input())
a = np.array(list(map(float, input().split())))
b, d = np.polyfit(np.array(range(n)), a, 1)
print(d, b)
print(sum([ (d + b * i - a_i) ** 2 for i, a_i in enumerate(a) ]))