結果
問題 | No.731 等差数列がだいすき |
ユーザー |
![]() |
提出日時 | 2018-09-07 21:51:20 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 34 ms / 1,500 ms |
コード長 | 693 bytes |
コンパイル時間 | 162 ms |
コンパイル使用メモリ | 12,800 KB |
実行使用メモリ | 11,008 KB |
最終ジャッジ日時 | 2024-11-29 12:50:02 |
合計ジャッジ時間 | 1,666 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 18 |
ソースコード
#!/usr/bin/env python3 def lsm(n, qs): d = n * sum(i * qs[i] for i in range(n)) - sum(range(n)) * sum(qs) d /= n * sum(i ** 2 for i in range(n)) - (sum(range(n))) ** 2 b0 = sum(i ** 2 for i in range(n)) * sum(qs) b0 -= sum(i * qs[i] for i in range(n)) * sum(range(n)) b0 /= n * sum(i ** 2 for i in range(n)) - (sum(range(n))) ** 2 bs = (b0 + i * d for i in range(n)) cost = sum((q - b) ** 2 for q, b in zip(qs, bs)) return b0, d, cost def main(): n = int(input()) qs = [int(q) for q in input().split()] b0, d, cost = lsm(n, qs) print("{:.12f} {:.12f}".format(b0, d)) print("{:.12f}".format(cost)) if __name__ == '__main__': main()