結果

問題 No.955 ax^2+bx+c=0
コンテスト
ユーザー ianCK
提出日時 2019-12-18 04:01:20
言語 C++17
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=c++17 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
RE  
実行時間 -
コード長 1,366 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 625 ms
コンパイル使用メモリ 102,564 KB
実行使用メモリ 7,972 KB
最終ジャッジ日時 2026-06-08 09:21:15
合計ジャッジ時間 6,260 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 82 RE * 40
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <stdio.h>
#include <iostream>
#include <math.h>
#include <assert.h>
#include <iomanip>
using namespace std;
constexpr long double kEps = 1E-16, nEps = 1E-12;
int main() {
	long long int a, b, c;
	scanf("%lld%lld%lld", &a, &b, &c);
	if (a == 0) {
		if (b == 0) {
			if (c == 0) printf("-1\n");
			else printf("0\n");
		}
		else printf("1\n%.20lf\n", (double)((long double)c / (long double) -b));
	}
	else {
		long double ta(a), tb(b), tc(c), d, sa, sb, l, r, mid;
		d = b * b - 4 * a * c;
		if (b * b - 4 * a * c > 0) {
			l = -d;
			r = d;
			for (int i = 0; i < 100000; i++) {
				mid = (l + r) / 2;
				if (mid * mid > d + kEps) r = mid;
				else if(mid * mid < d - kEps) l = mid;
				else break;
			}
			mid = (l + r) / 2;
			sa = (-tb + l) / ta / 2; 
			sb = (-tb - l) / ta / 2; 
			if (sa < sb) {
				cout << 2 << '\n';
				cout << setprecision(80) << sa << '\n';
				cout << setprecision(80) << sb << '\n';
			}
			else {
				cout << 2 << '\n';
				cout << setprecision(80) << sb << '\n';
				cout << setprecision(80) << sa << '\n';
			}
			sa = a * sa * sa + b * sa + c;
			assert(-nEps < sa && sa < nEps);
			sb = a * sb * sb + b * sb + c;
			assert(-nEps < sb && sb < nEps);
		}
		else if (b * b - 4 * a * c < 0) printf("0\n");
		else printf("1\n%.20lf\n", (double) (-tb / (ta * 2)));
		
	}
	
}
0