結果
問題 | No.731 等差数列がだいすき |
ユーザー |
|
提出日時 | 2018-09-07 21:42:55 |
言語 | D (dmd 2.109.1) |
結果 |
AC
|
実行時間 | 2 ms / 1,500 ms |
コード長 | 923 bytes |
コンパイル時間 | 2,720 ms |
コンパイル使用メモリ | 166,784 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-06-13 01:38:28 |
合計ジャッジ時間 | 3,502 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 18 |
ソースコード
import std.stdio, std.array, std.string, std.conv, std.algorithm;import std.typecons, std.range, std.random, std.math, std.container;import std.numeric, std.bigint, core.bitop, std.bitmanip;immutable long MOD = 10^^9 + 7;immutable long INF = 1L << 59;void main() {auto N = readln.chomp.to!int;auto X = N.iota.map!(to!real).array;auto Y = readln.split.map!(to!real).array;real E(const real[] A) {return A.sum / N;}real S(const real[] A, real e) {return sqrt(A.map!(a => (a - e) * (a - e)).sum / N);}real Ex = E(X);real Ey = E(Y);real Sx = S(X, Ex);real Sy = S(Y, Ey);real cov = E(N.iota.map!(i => (X[i] - Ex) * (Y[i] - Ey)).array);auto A = cov / (Sx * Sx);auto B = Ey - A * Ex;real ans = N.iota.map!(i => (Y[i] - A * X[i] - B) * (Y[i] - A * X[i] - B)).sum;writefln("%.9f %.9f", B, A);writefln("%.9f", ans);}