結果

問題 No.955 ax^2+bx+c=0
ユーザー ianCK
提出日時 2019-12-18 04:21:27
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 3 ms / 2,000 ms
コード長 1,738 bytes
コンパイル時間 1,062 ms
コンパイル使用メモリ 67,356 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-18 22:08:12
合計ジャッジ時間 3,680 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 122
権限があれば一括ダウンロードができます
コンパイルメッセージ
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) {
			if (a >= 100000000) {
				d = sqrtl(d);
				sa = (-tb + d) / ta / 2; 
				sb = (-tb - d) / ta / 2; 
				if (sa < sb) printf("2\n%.20lf\n%.20lf\n", (double) sa, (double) sb);
				else printf("2\n%.20lf\n%.20lf\n", (double) sb, (double) sa);
				return 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