結果
問題 | No.731 等差数列がだいすき |
ユーザー |
![]() |
提出日時 | 2018-09-07 22:01:05 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 1,500 ms |
コード長 | 2,286 bytes |
コンパイル時間 | 1,470 ms |
コンパイル使用メモリ | 172,884 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-29 13:29:38 |
合計ジャッジ時間 | 2,139 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 18 |
ソースコード
#include<bits/stdc++.h>#define rep(i,a,b) for(int i=a;i<b;i++)#define rrep(i,a,b) for(int i=a;i>=b;i--)#define fore(i,a) for(auto &i:a)#define all(x) (x).begin(),(x).end()#pragma GCC optimize ("-O3")using namespace std; void _main(); int main() { cin.tie(0); ios::sync_with_stdio(false); _main(); }typedef long long ll; const int inf = INT_MAX / 2; const ll infl = 1LL << 60;template<class T>bool chmax(T &a, const T &b) { if (a<b) { a = b; return 1; } return 0; }template<class T>bool chmin(T &a, const T &b) { if (b<a) { a = b; return 1; } return 0; }//---------------------------------------------------------------------------------------------------/*---------------------------------------------------------------------------------------------------∧_∧∧_∧ (´<_` ) Welcome to My Coding Space!( ´_ゝ`) / ⌒i/ \ | |/ / ̄ ̄ ̄ ̄/ |__(__ニつ/ _/ .| .|____\/____/ (u ⊃---------------------------------------------------------------------------------------------------*/int N, A[1010];//---------------------------------------------------------------------------------------------------pair<double, double> linearApproximation(vector<pair<int, int>> v) {int n = v.size();double xy = 0;rep(i, 0, n) xy += v[i].first * v[i].second;double x = 0;rep(i, 0, n) x += v[i].first;double y = 0;rep(i, 0, n) y += v[i].second;double xx = 0;rep(i, 0, n) 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 };}//---------------------------------------------------------------------------------------------------void _main() {cin >> N;rep(i, 0, N) cin >> A[i];vector<pair<int, int>> v;rep(i, 0, N) v.push_back({ i, A[i] });double a, b;tie(a, b) = linearApproximation(v);printf("%.10f %.10f\n", b, a);double c = 0;rep(i, 0, N) {double y = a * i + b;double d = y - A[i];c += d * d;}printf("%.10f\n", c);}