結果

問題 No.955 ax^2+bx+c=0
ユーザー ianCKianCK
提出日時 2019-12-18 04:16:17
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,468 bytes
コンパイル時間 776 ms
コンパイル使用メモリ 91,584 KB
最終ジャッジ日時 2025-01-08 12:18:34
ジャッジサーバーID
(参考情報)
judge1 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 121 WA * 1
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:10:14: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   10 |         scanf("%lld%lld%lld", &a, &b, &c);
      |         ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~

ソースコード

diff #

#include <stdio.h>
#include <iostream>
#include <math.h>
#include <assert.h>
#include <iomanip>
using namespace std;
constexpr long double kEps = 1E-13, nEps = 1E-12, kInf = 1E20;
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 {
		if (a < 0) {
			a = -a;
			b = -b;
			c = -c;
		}
		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) {
			cout << 2 << '\n';
			d = sqrtl(d);
			l = -kInf;
			r = -tb / ta / 2 + kEps;
			for (int i = 0; i < 100000; i++) {
				mid = (l + r) / 2;
				if (ta * mid * mid + tb * mid + tc > kEps) l = mid;
				else if(ta * mid * mid + tb * mid + tc < - kEps) r = mid;
			}
			mid = (l + r) / 2;
			cout << setprecision(2000) << mid << '\n';
			sa = a * mid * mid + b * mid + c;
			l = -tb / ta / 2 - kEps;
			r = kInf;
			for (int i = 0; i < 100000; i++) {
				mid = (l + r) / 2;
				if (ta * mid * mid + tb * mid + tc > kEps) r = mid;
				else if(ta * mid * mid + tb * mid + tc < - kEps) l = mid;
			}
			mid = (l + r) / 2;
			cout << setprecision(2000) << mid << '\n';
			sa = a * mid * mid + b * mid + c;
		}
		else if (b * b - 4 * a * c < 0) printf("0\n");
		else printf("1\n%.20lf\n", (double) (-tb / (ta * 2)));
		
	}
	
}
0