結果

問題 No.731 等差数列がだいすき
コンテスト
ユーザー むらため
提出日時 2019-01-27 04:11:45
言語 Nim
(2.2.8)
コンパイル:
nim --nimcache=~ --hints:off -o:a.out -d:release cpp _filename_
実行:
./a.out
結果
AC  
実行時間 2 ms / 1,500 ms
コード長 782 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,232 ms
コンパイル使用メモリ 70,516 KB
実行使用メモリ 7,844 KB
最終ジャッジ日時 2026-03-22 10:13:46
合計ジャッジ時間 3,039 ms
ジャッジサーバーID
(参考情報)
judge2_1 / judge1_0
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 18
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code

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