結果

問題 No.724 円と円の間の円
ユーザー ats5515ats5515
提出日時 2017-10-04 09:27:46
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,480 bytes
コンパイル時間 693 ms
コンパイル使用メモリ 78,340 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-05-01 13:10:11
合計ジャッジ時間 1,392 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

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

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <map>
#include <set>
#include <string>
#include <algorithm>
#include <cmath>
#include <cassert>
#include <iomanip>
using namespace std;
string tos(int a) {
	string res = "";
	if (a == 0) {
		res = '0';
	}
	while (a > 0) {
		res = char((a % 10) + '0') + res;
		a /= 10;
	}
	return res;
}
int n, a, b;

double calcY(double x) {
	return sqrt(a*b - (4 * a*b*(x - ((double)(a + b) / 2))*(x - ((double)(a + b) / 2)) / ((a + b)*(a + b))));
}
double calcR(double x,double y) {
	return sqrt((x - a)*(x - a) + y*y) - a;
}
double dst(double x1, double y1, double x2, double y2) {
	return sqrt((x1 - x2)*(x1 - x2) + (y1 - y2)*(y1 - y2));
}
int main() {

	cin >> n >> a >> b;
	//cout << setprecision(15) << 4 * a*b*(b - a) / (double)((((double)n*n - 2 * n)* (b - a)*(b - a)) + (a + b)*(a + b)) << endl;
	double u;
	double d;
	double m;
	double x;
	double y;
	double t;
	double r;
	if(n % 2 == 0){
		u = a + b;
		d = 2 * a;
		while(u - d > 1e-9) {
			m = (u + d) * 0.5;
			t = calcY(m);
			if (t > calcR(m, t)) {
				d = m;
			}
			else {
				u = m;
			}
		}
		x = (u + d)*0.5;
		y = calcY(x);
	}
	else {
		x = a + b;
		y = 0;
	}
	n = (n - 1) / 2;
	r = calcR(x, y);
	while (n--) {
		u = x;
		d = 0;
		
		while (u - d > 1e-9) {
			m = (u + d) * 0.5;
			t = calcY(m);
			if (dst(m, t, x, y) > r + calcR(m, t)) {
				d = m;
			}
			else {
				u = m;
			}
		}
		x = (u + d)*0.5;
		y = calcY(x);
		r = calcR(x, y);
	}
	cout << r << endl;
	
	return 0;
}
0