結果

問題 No.764 浮動点
ユーザー square1001square1001
提出日時 2019-03-27 22:04:36
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 8 ms / 1,500 ms
コード長 1,945 bytes
コンパイル時間 769 ms
コンパイル使用メモリ 83,416 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-10-11 06:13:21
合計ジャッジ時間 1,971 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
6,820 KB
testcase_02 AC 2 ms
6,816 KB
testcase_03 AC 2 ms
6,820 KB
testcase_04 AC 6 ms
6,816 KB
testcase_05 AC 3 ms
6,816 KB
testcase_06 AC 6 ms
6,816 KB
testcase_07 AC 3 ms
6,816 KB
testcase_08 AC 7 ms
6,816 KB
testcase_09 AC 3 ms
6,820 KB
testcase_10 AC 7 ms
6,816 KB
testcase_11 AC 3 ms
6,816 KB
testcase_12 AC 7 ms
6,816 KB
testcase_13 AC 7 ms
6,816 KB
testcase_14 AC 7 ms
6,820 KB
testcase_15 AC 7 ms
6,820 KB
testcase_16 AC 7 ms
6,820 KB
testcase_17 AC 7 ms
6,816 KB
testcase_18 AC 2 ms
6,816 KB
testcase_19 AC 6 ms
6,816 KB
testcase_20 AC 8 ms
6,816 KB
testcase_21 AC 8 ms
6,816 KB
testcase_22 AC 7 ms
6,816 KB
testcase_23 AC 8 ms
6,820 KB
testcase_24 AC 7 ms
6,820 KB
testcase_25 AC 2 ms
6,820 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cmath>
#include <vector>
#include <iostream>
#include <algorithm>
using namespace std;
const long double pi = acos(-1);
pair<int, int> calc(vector<int> v) {
	int sum = 0, mx = 0;
	for (int i = 0; i < v.size(); ++i) {
		sum += v[i];
		mx = max(mx, v[i]);
	}
	int L = 0, R = sum;
	if (mx * 2 > sum) L = mx * 2 - sum;
	return make_pair(L, R);
}
long double area(long double a, long double b, long double c) {
	if (a + b <= c) return 0;
	if (a < b) swap(a, b);
	if (a > b + c) return b * b * pi;
	long double l = 0, r = min(a, b);
	if (a * a <= b * b + c * c) {
		for (int i = 0; i < 60; ++i) {
			long double m = (l + r) * 0.5;
			long double x = sqrt(a * a - m * m) + sqrt(b * b - m * m);
			if (x >= c) l = m;
			else r = m;
		}
	}
	else {
		for (int i = 0; i < 60; ++i) {
			long double m = (l + r) * 0.5;
			long double x = sqrt(a * a - m * m) - sqrt(b * b - m * m);
			if (x >= c) r = m;
			else l = m;
		}
	}
	long double rad1 = asin(l / a);
	long double rad2 = (a * a <= b * b + c * c ? asin(l / b) : pi - asin(l / b));
	long double circle1 = rad1 * a * a;
	long double circle2 = rad2 * b * b;
	long double ans = circle1 + circle2 - c * l;
	return ans;
}
int main() {
	int N;
	cin >> N;
	vector<int> A(N + 2);
	int sum = 0, mx = 0;
	for (int i = 0; i < N + 2; ++i) {
		cin >> A[i];
	}
	cout.precision(15);
	pair<int, int> resbase = calc(A);
	if (resbase.first > 0) {
		for (int i = 0; i < N; ++i) {
			cout << 0 << '\n';
		}
	}
	else {
		for (int i = 1; i <= N; ++i) {
			vector<int> vl, vr;
			for (int j = 1; j <= i; ++j) {
				vl.push_back(A[j]);
			}
			for (int j = i + 1; j <= N + 1; ++j) {
				vr.push_back(A[j]);
			}
			pair<int, int> wl = calc(vl);
			pair<int, int> wr = calc(vr);
			long double ans = 0;
			ans += area(wl.second, wr.second, A[0]);
			ans += area(wl.first, wr.first, A[0]);
			ans -= area(wl.second, wr.first, A[0]);
			ans -= area(wl.first, wr.second, A[0]);
			cout << ans << '\n';
		}
	}
	return 0;
}
0