結果
問題 |
No.731 等差数列がだいすき
|
ユーザー |
|
提出日時 | 2019-01-27 04:11:45 |
言語 | Nim (2.2.0) |
結果 |
AC
|
実行時間 | 2 ms / 1,500 ms |
コード長 | 782 bytes |
コンパイル時間 | 2,718 ms |
コンパイル使用メモリ | 62,192 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-07-01 11:00:10 |
合計ジャッジ時間 | 3,661 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 18 |
ソースコード
import sequtils,math proc getchar_unlocked():char {. importc:"getchar_unlocked",header: "<stdio.h>" .} proc scan(): int = while true: let k = getchar_unlocked() if k < '0': break result = 10 * result + k.ord - '0'.ord let n = scan() let A = newSeqWith(n,scan().float) proc calcCost(b,d:float):float = for i,a in A: let bi = b + d * i.float result += (a - bi)*(a - bi) let S = toSeq(0..<n) let d1 = n.float * S.mapIt(it.float * A[it]).sum() let d2 = S.sum().float * A.sum() let d3 = n.float * S.mapIt(it * it).sum().float let d4 = S.sum().float * S.sum().float let d = (d1 - d2) / (d3 - d4) let b1 = S.mapIt(it * it).sum().float * A.sum() let b2 = S.mapIt(it.float * A[it]).sum() * S.sum().float let b = (b1 - b2) / (d3 - d4) echo b," ",d echo calcCost(b,d)