結果
| 問題 |
No.731 等差数列がだいすき
|
| ユーザー |
ixmel
|
| 提出日時 | 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を求める問題
自分は偏微分した
*/
ixmel