結果

問題 No.731 等差数列がだいすき
ユーザー 0w10w1
提出日時 2018-09-09 00:05:35
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 2 ms / 1,500 ms
コード長 1,227 bytes
コンパイル時間 2,075 ms
コンパイル使用メモリ 205,476 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-08-23 10:32:43
合計ジャッジ時間 3,189 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 2 ms
4,380 KB
testcase_02 AC 1 ms
4,376 KB
testcase_03 AC 2 ms
4,376 KB
testcase_04 AC 2 ms
4,380 KB
testcase_05 AC 2 ms
4,376 KB
testcase_06 AC 2 ms
4,376 KB
testcase_07 AC 2 ms
4,376 KB
testcase_08 AC 2 ms
4,380 KB
testcase_09 AC 2 ms
4,380 KB
testcase_10 AC 2 ms
4,376 KB
testcase_11 AC 1 ms
4,380 KB
testcase_12 AC 2 ms
4,380 KB
testcase_13 AC 2 ms
4,376 KB
testcase_14 AC 1 ms
4,376 KB
testcase_15 AC 1 ms
4,376 KB
testcase_16 AC 1 ms
4,380 KB
testcase_17 AC 2 ms
4,376 KB
testcase_18 AC 2 ms
4,376 KB
testcase_19 AC 2 ms
4,376 KB
testcase_20 AC 2 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

pair<double, double> least_squares(const vector<int> &x, const vector<int> &y) {
  int n = x.size();
  vector<double> xy(n);
  for (int i = 0; i < n; ++i) xy[i] = 1.0 * x[i] * y[i];
  vector<double> x2(n);
  for (int i = 0; i < n; ++i) x2[i] = 1.0 * x[i] * x[i];
  double xsum = accumulate(x.begin(), x.end(), 0.0);
  double ysum = accumulate(y.begin(), y.end(), 0.0);
  double xysum = accumulate(xy.begin(), xy.end(), 0.0);
  double x2sum = accumulate(x2.begin(), x2.end(), 0.0);
  double au = n * xysum - xsum * ysum;
  double av = n * x2sum - xsum * xsum;
  double bu = x2sum * ysum - xysum * xsum;
  double bv = n * x2sum - xsum * xsum;
  return {au / av, bu / bv};
}

signed main() {
  ios::sync_with_stdio(false);

  int N;
  cin >> N;

  vector<int> A(N);
  for (int i = 0; i < N; ++i) cin >> A[i];

  vector<int> X(N);
  iota(X.begin(), X.end(), 0);
  double a, b;
  tie(a, b) = least_squares(X, A);

  vector<double> v(N);
  v[0] = b;
  for (int i = 0; i + 1 < N; ++i) v[i + 1] = v[i] + a;

  double c = 0;
  for (int i = 0; i < N; ++i) {
    double x = v[i] - A[i];
    c += x * x;
  }

  cout << fixed << setprecision(8) << b << ' ' << a << '\n' << ' ' << c << endl;
}
0