結果
| 問題 | 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)