結果
問題 |
No.955 ax^2+bx+c=0
|
ユーザー |
![]() |
提出日時 | 2019-12-18 09:43:57 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
RE
|
実行時間 | - |
コード長 | 1,081 bytes |
コンパイル時間 | 633 ms |
コンパイル使用メモリ | 77,908 KB |
実行使用メモリ | 6,948 KB |
最終ジャッジ日時 | 2024-07-07 00:36:20 |
合計ジャッジ時間 | 5,755 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 101 RE * 21 |
ソースコード
#include <cassert> #include <cmath> #include <iomanip> #include <iostream> using namespace std; double bs(double pos, double neg, double a, double b, double c, double eps) { assert(a*pos*pos+b*pos+c>0); assert(a*neg*neg+b*neg+c<0); for (int i = 0; i < 1000; i++) { double m = (pos + neg) / 2; if (a*m*m+b*m+c < 0) neg = m; else pos = m; } return (pos+neg)/2; } int main() { long long a, b, c; cin >> a >> b >> c; cout << setprecision(16); if (a == 0) { if (b == 0) { if (c == 0) cout << "-1\n"; else cout << "0\n"; } else { cout << "1\n" << -(double)c/b << endl; } } else { long long d = b * b - 4 * a * c; if (d == 0) { cout << "1\n" << -b/2.0/a << endl; } else if (d < 0) { cout << "0\n"; } else { double v = -b/2.0/a, x = bs(v-d+1, v, a, b, c, 1e-11), y = bs(v+d+1, v, a, b, c, 1e-11); cout << "2\n" << x << endl << y << endl; } } }