結果
問題 |
No.731 等差数列がだいすき
|
ユーザー |
![]() |
提出日時 | 2018-09-22 18:41:32 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 2 ms / 1,500 ms |
コード長 | 1,156 bytes |
コンパイル時間 | 1,128 ms |
コンパイル使用メモリ | 161,144 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-07-18 10:18:25 |
合計ジャッジ時間 | 1,784 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 18 |
ソースコード
#include <bits/stdc++.h> using namespace std; typedef long long ll; typedef vector<int> vint; typedef pair<int,int> pint; typedef vector<pint> vpint; #define rep(i,n) for(int i=0;i<(n);i++) #define REP(i,n) for(int i=n-1;i>=(0);i--) #define reps(i,f,n) for(int i=(f);i<(n);i++) #define each(it,v) for(__typeof((v).begin()) it=(v).begin();it!=(v).end();it++) #define all(v) (v).begin(),(v).end() #define eall(v) unique(all(v), v.end()) #define pb push_back #define mp make_pair #define fi first #define se second #define chmax(a, b) a = (((a)<(b)) ? (b) : (a)) #define chmin(a, b) a = (((a)>(b)) ? (b) : (a)) const int MOD = 1e9 + 7; const int INF = 1e9; const ll INFF = 1e18; int main() { int n; cin >> n; vector<double> a(n); { rep(i, n) cin >> a[i]; } double x = 0, x2 = 0, y = 0, y2 = 0, xy = 0; rep(i, n) { x += i; x2 += i * i; y += a[i]; xy += i * a[i]; } double ans_a = (n * xy - x * y) / (n * x2 - x * x); double ans_b = (x2 * y - xy * x) / (n * x2 - x * x); double ans = 0; { rep(i, n) { double act = ans_a * i + ans_b; ans += pow(act - a[i], 2); } } printf("%.9f %9f\n", ans_b, ans_a); printf("%.9f\n", ans); }