結果

問題 No.731 等差数列がだいすき
ユーザー donkorin_donkorin_
提出日時 2018-09-09 14:55:27
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 2 ms / 1,500 ms
コード長 1,437 bytes
コンパイル時間 1,437 ms
コンパイル使用メモリ 167,080 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-08-24 23:14:30
合計ジャッジ時間 2,686 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

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

ソースコード

diff # 

#include <bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<(b);++i)
#define erep(i,a,b) for(int i=a;i<=(int)(b);++i)
#define per(i,a,b) for(int i=(a);i>(b);--i)
#define eper(i,a,b) for(int i=(a);i>=b;--i)
#define pb push_back
#define mp make_pair
#define INF (1<<30)-1
#define MOD 1000000007
#define all(x) (x).begin(),(x).end()
#define vii vector<int>
#define vll vector<long long>
using namespace std;
typedef long long ll;
typedef pair<int,int> Pii;
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; }
int dy[]={0, 0, 1, -1};
int dx[]={1, -1, 0, 0};
int gcd(int a,int b){return b?gcd(b,a%b):a;}
int lcm(int a,int b){return a/gcd(a, b)*b;}

int n;
int main() {
  cin.tie(0);
  ios::sync_with_stdio(false);
    cin >> n;
    vector<double> a(n), x(n);
    rep(i, 0, n) {
        cin >> a[i];
        x[i] = i;
    }
    double xsum = 0, ysum = 0, xysum = 0, x2sum = 0;
    rep(i, 0, n) {
        xsum += i;
        ysum += a[i];
        xysum += i * a[i];
        x2sum += i * i;
    }
    double d = (n * xysum - xsum * ysum) / (n * x2sum - (xsum * xsum));
    double b1 = (x2sum * ysum - xysum * xsum) / (n * x2sum - (xsum * xsum));
    double cost = 0;
    rep(i, 0, n) {
        cost += (a[i] - (b1 + i * d)) * (a[i] - (b1 + i * d));
    }
    printf("%.15f %.15f\n%.15f\n", b1, d, cost);
    return 0;
}
0