結果
問題 | No.731 等差数列がだいすき |
ユーザー |
![]() |
提出日時 | 2018-09-07 21:34:22 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 6 ms / 1,500 ms |
コード長 | 801 bytes |
コンパイル時間 | 1,800 ms |
コンパイル使用メモリ | 170,316 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-29 05:45:23 |
合計ジャッジ時間 | 2,696 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 18 |
ソースコード
#include<bits/stdc++.h>using namespace std;using Int = long long;using Double = long double;struct Precision{Precision(){cout<<fixed<<setprecision(12);}}precision_beet;//INSERT ABOVE HEREsigned main(){int n;cin>>n;vector<Double> a(n);for(int i=0;i<n;i++) cin>>a[i];auto calc=[&](Double x,bool f=0)->Double{auto v(a);for(int i=0;i<n;i++) v[i]-=i*x;Double b=0;for(int i=0;i<n;i++) b+=v[i];b/=n;Double c=0;for(int i=0;i<n;i++) c+=(v[i]-b)*(v[i]-b);if(f) cout<<b<<" "<<x<<endl;return c;};Double l=-1e9,r=1e9;for(int k=0;k<200;k++){Double m1=l+(r-l)/3;Double m2=l+(r-l)*2/3;if(calc(m1)<calc(m2)) r=m2;else l=m1;}Double c=calc(l,1);cout<<c<<endl;return 0;}