結果

問題 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
権限があれば一括ダウンロードができます

ソースコード

diff #

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)
0