結果
問題 |
No.731 等差数列がだいすき
|
ユーザー |
![]() |
提出日時 | 2018-09-07 22:36:53 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,019 bytes |
コンパイル時間 | 1,552 ms |
コンパイル使用メモリ | 168,976 KB |
実行使用メモリ | 6,824 KB |
最終ジャッジ日時 | 2024-11-29 16:08:35 |
合計ジャッジ時間 | 7,953 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 WA * 1 |
other | WA * 18 |
ソースコード
#include <bits/stdc++.h> using namespace std; typedef unsigned long long ul; typedef signed long long ll; ul over = 1000000007; mt19937 mt(time(nullptr)); // uniform_int_distribution< int > rand_int(0, 10) // rand[0, 10] inclusive int main(void) { cin.tie(0); ios::sync_with_stdio(false); cout << fixed; int n; cin >> n; vector< double > a(n); for (int i=0; i<n; ++i) { cin >> a[i]; } //double b=a[0], d=(a[n-1]-a[0])/(n-1); double b=0.01, d=0.01; double pdb=0, pdd=0; ul max_ep = 100000000/n; double rate = 1.0; for (ul ep=0; ep<max_ep; ++ep) { double db=0, dd=0; for (int i=0; i<n; ++i) { db += 2*(b+d*i-a[i]); dd += 2*(b+d*i-a[i]) * i; } //db += pdb/2; //dd += pdd/2; pdb=db; pdd=dd; b -= rate * db/n * 0.1; d -= rate * dd/n * 0.1; rate -= 0.5/max_ep; } cout << b << " " << d << endl; double cost = 0; for (ul i=0; i<n; ++i) { double diff = b+d*i-a[i]; cost += diff*diff; } cout << cost << endl; return 0; }