結果

問題 No.731 等差数列がだいすき
ユーザー 0w10w1
提出日時 2018-09-09 00:04:27
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 1,219 bytes
コンパイル時間 2,229 ms
コンパイル使用メモリ 204,532 KB
実行使用メモリ 4,504 KB
最終ジャッジ日時 2023-08-23 10:31:50
合計ジャッジ時間 3,370 ms
ジャッジサーバーID
(参考情報)
judge13 / judge15
このコードへのチャレンジ
(要ログイン)

テストケース

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

ソースコード

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 << ' ' << c << endl;
}
0