#define _USE_MATH_DEFINES #define _CRT_SECURE_NO_WARNINGS #include using namespace std; /*BigInteger #include #include #include namespace xxx = boost::multiprecision; using Bint = xxx::cpp_int; using Real = xxx::number>; */ #define int long long #define pb(x) push_back(x) #define m0(x) memset((x), 0LL, sizeof(x)) #define mm(x) memset((x), -1LL, sizeof(x)) //container #define ALL(x) (x).begin(), (x).end() #define RALL(a) (a).rbegin(), (a).rend() #define EACH(i, c) for (typeof((c).begin()) i = (c).begin(); i != (c).end(); ++i) #define EXIST(s, e) ((s).find(e) != (s).end()) #define UNIQUE(v) (v).erase(unique((v).begin(), (v).end()), (v).end()); #define PERM(c) \ sort(ALL(c)); \ for (bool c##p = 1; c##p; c##p = next_permutation(ALL(c))) // debug #define GET_VAR_NAME(variable) #variable #define test(x) cout << GET_VAR_NAME(x) << " = " << x << endl; // bit_macro #define Bit(n) (1LL << (n)) #define Bitset(a, b) (a) |= (1LL << (b)) #define Bitunset(a, b) (a) &= ~(1LL << (b)) #define Bitcheck(a, b) ((((a) >> (b)) & 1LL) == 1LL) #define Bitcount(a) __builtin_popcountll((a)) //typedef typedef long long lint; typedef unsigned long long ull; typedef complex Complex; typedef pair P; typedef tuple TP; typedef vector vec; typedef vector mat; //constant constexpr int INF = (int)1e18; constexpr int MOD = (int)1e9 + 7; constexpr double PI = (double)acos(-1); constexpr double EPS = (double)1e-10; constexpr int dx[] = {-1, 0, 0, 1, 0, -1, -1, 1, 1}; constexpr int dy[] = {0, -1, 1, 0, 0, 1, -1, 1, -1}; // template void chmax(T &a, T b) { a = max(a, b); } template void chmin(T &a, T b) { a = min(a, b); } // inline int toInt(string s) { int v; istringstream sin(s); sin >> v; return v; } template inline string toString(T x) { ostringstream sout; sout << x; return sout.str(); } // struct Accelerate_Cin { Accelerate_Cin() { cin.tie(0); ios::sync_with_stdio(0); cout << fixed << setprecision(20); }; }; //O(N^2) //最小二乗法 //N個の点を直線の方程式: y=ax+b で近似する。 //N個の点:(v[i].fir,v[i].sec) pair linearApproximation(vector> v) { int n = v.size(); double xy = 0, x = 0, y = 0, xx = 0; for (int i = 0; i < n; i++) xy += v[i].first * v[i].second; for (int i = 0; i < n; i++) x += v[i].first; for (int i = 0; i < n; i++) y += v[i].second; for (int i = 0; i < n; i++) xx += v[i].first * v[i].first; double a = (n * xy - x * y) / (n * xx - x * x); double b = (xx * y - xy * x) / (n * xx - x * x); return {a, b}; } signed main() { int N; cin >> N; vec A(N); for (int i = 0; i < N; i++) cin >> A[i]; vector

v; for (int i = 0; i < N; i++) v.push_back({i, A[i]}); double a, b; tie(a, b) = linearApproximation(v); cout << fixed << setprecision(12) << b << " " << a << endl; double c = 0; for (int i = 0; i < N; i++) { double y = a * i + b; double d = y - A[i]; c += d * d; } cout << fixed << setprecision(12) << c << endl; return 0; }