// includes #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include // macros #define ll long long int #define pb emplace_back #define mk make_pair #define pq priority_queue #define FOR(i, a, b) for(int i=(a);i<(b);++i) #define rep(i, n) FOR(i, 0, n) #define rrep(i, n) for(int i=((int)(n)-1);i>=0;i--) #define all(x) (x).begin(),(x).end() #define sz(x) ((int)(x).size()) #define UNIQUE(v) v.erase(unique(v.begin(), v.end()), v.end()) using namespace std; // types typedef pair P; typedef pair Pl; typedef pair Pll; typedef pair Pd; // constants const int inf = 1e9; const ll linf = 1LL << 50; const double EPS = 1e-10; const int mod = 1e9 + 7; // solve template bool chmax(T &a, const T &b){if(a < b){a = b; return 1;} return 0;} template bool chmin(T &a, const T &b){if(a > b){a = b; return 1;} return 0;} int main(int argc, char const* argv[]) { ios_base::sync_with_stdio(false); cin.tie(0); int n; cin >> n; double N = n; vector a(n); rep(i, n)cin >> a[i]; double sum = 0., ssum = 0., psum = 0.; rep(i, n){ sum += a[i]; ssum += a[i] * a[i]; psum += i * a[i]; } cout << setprecision(15); double d = - (sum * (N - 1.) - 2. * psum) / (N * (N * N - 1.) / 12.) / 2.; double b = - (- 2. * sum + N * (N - 1.) * d) / (2. * N); cout << b << " " << d << endl; double res = 0.; rep(i, n){ res += (a[i] - b - i * d) * (a[i] - b - i * d); } cout << res << endl; return 0; }