結果

問題 No.2352 Sharpened Knife in Fall
ユーザー norikamenorikame
提出日時 2023-06-16 23:19:51
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 355 ms / 3,000 ms
コード長 1,501 bytes
コンパイル時間 5,059 ms
コンパイル使用メモリ 311,840 KB
実行使用メモリ 4,872 KB
最終ジャッジ日時 2023-09-06 22:19:04
合計ジャッジ時間 15,361 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,376 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 2 ms
4,376 KB
testcase_03 AC 1 ms
4,376 KB
testcase_04 AC 233 ms
4,624 KB
testcase_05 AC 1 ms
4,380 KB
testcase_06 AC 354 ms
4,636 KB
testcase_07 AC 177 ms
4,376 KB
testcase_08 AC 354 ms
4,588 KB
testcase_09 AC 355 ms
4,580 KB
testcase_10 AC 352 ms
4,872 KB
testcase_11 AC 345 ms
4,624 KB
testcase_12 AC 350 ms
4,620 KB
testcase_13 AC 349 ms
4,572 KB
testcase_14 AC 159 ms
4,380 KB
testcase_15 AC 316 ms
4,624 KB
testcase_16 AC 18 ms
4,376 KB
testcase_17 AC 9 ms
4,376 KB
testcase_18 AC 125 ms
4,380 KB
testcase_19 AC 278 ms
4,568 KB
testcase_20 AC 107 ms
4,380 KB
testcase_21 AC 105 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

// 

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;

using ll = long long;

#define rep(i, n) for (int i=0; i<(int)(n); ++(i))
#define rep3(i, m, n) for (int i=(m); (i)<(int)(n); ++(i))
#define repr(i, n) for (int i=(int)(n)-1; (i)>=0; --(i))
#define rep3r(i, m, n) for (int i=(int)(n)-1; (i)>=(int)(m); --(i))
#define all(x) (x).begin(), (x).end()

const double PI = acos(-1);

int main() {
	int r, k;
	cin >> r >> k;
	if (k%2 == 1) {
		vector<double> res;
		double tar = ((double)r * r * PI / 2.0) / (k/2+1);
		res.push_back(0.0);
		rep3(i, 1, k/2+1) {
			double li = 0.0, ri = r;
			while (ri-li > 1e-5) {
				double ci = (li + ri) / 2.0, th = asin(ci/r);
				double area = ((double)r * r * PI / 2.0) * (th / PI) * 2.0 + r * cos(th) * 2.0 * ci / 2.0;
				if (area < tar*i) li = ci;
				else ri = ci;
			}
			res.push_back(ri);
			res.push_back(-ri);
		}
		sort(all(res));
		rep(i, k) cout << fixed << setprecision(10) << res[i] << endl;
	}
	else {
		vector<double> res;
		double tar = ((double)r * r * PI / 2.0) / (k+1);
		for (int i=1; i<k; i+=2) {
			double li = 0.0, ri = r;
			while (ri-li > 1e-5) {
				double ci = (li + ri) / 2.0, th = asin(ci/r);
				double area = ((double)r * r * PI / 2.0) * (th / PI) * 2.0 + r * cos(th) * 2.0 * ci / 2.0;
				if (area < tar*i) li = ci;
				else ri = ci;
			}
			res.push_back(ri);
			res.push_back(-ri);
		}
		sort(all(res));
		rep(i, k) cout << fixed << setprecision(10) << res[i] << endl;
	}
	return 0;
}
0