結果
問題 | 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 python3def 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))) ** 2b0 = 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))) ** 2bs = (b0 + i * d for i in range(n))cost = sum((q - b) ** 2 for q, b in zip(qs, bs))return b0, d, costdef 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()