結果
問題 |
No.731 等差数列がだいすき
|
ユーザー |
![]() |
提出日時 | 2017-04-05 00:35:19 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 2 ms / 1,500 ms |
コード長 | 724 bytes |
コンパイル時間 | 680 ms |
コンパイル使用メモリ | 64,348 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-07-08 09:09:59 |
合計ジャッジ時間 | 1,313 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 18 |
ソースコード
#include<iostream> #include<vector> #include<cmath> #include<cstdio> #include<assert.h> #define loop(i,a,b) for(int i=a;i<b;i++) #define rep(i,a) loop(i,0,a) using namespace std; typedef long long ll; typedef vector<ll> vi; int main(){ ll n; cin>>n; assert(2<=n&&n<=1000); vi in(n); rep(i,n)cin>>in[i]; rep(i,n)assert(0<=in[i]&&in[i]<=100000); ll n1=n*(n-1)/2,n2=n*(n-1)*(2*n-1)/6; ll sum1=0,sum2=0; rep(i,n){ sum1+=in[i]; sum2+=in[i]*i; } double b=(double)(n2*sum1-n1*sum2)/(n*n2-n1*n1); double d=(double)(n*sum2-n1*sum1)/(n*n2-n1*n1); double c=0; rep(i,n)c+=pow(in[i]-(d*i+b),2); printf("%.7f %.7f\n%.7f\n",b,d,c); } /* sum(d*i+b-a_i)が最小な時のb, dを求める問題 自分は偏微分した */