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