、

xy = [x[i]*y[i] for i in range(len(x))]
x2 = [x[i]**2 for i in range(len(x))]

x_ = sum(x)/n
y_ = sum(y)/n

a = (sum(xy) - n*x_*y_) / (sum(x2)-n*(x_**2))
b = y_ - a*x_

print("%f\t%f" % (a,b))