#882195 (C++23) No.2352 Sharpened Knife in Fall

• 問題: No.2352 Sharpened Knife in Fall
• ユーザー: norikame
• 提出日時: 2023-06-16 23:19:51
• 言語: C++23 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 344 ms / 3,000 ms
• コード長: 1,501 bytes
• コンパイル時間: 5,471 ms
• コンパイル使用メモリ: 313,316 KB
• 実行使用メモリ: 5,376 KB
• 最終ジャッジ日時: 2024-06-24 16:29:54
• 合計ジャッジ時間: 15,002 ms
• ジャッジサーバーID （参考情報）: judge1 / judge5

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 5,248 KB
testcase_01	AC	2 ms; 5,376 KB
testcase_02	AC	2 ms; 5,376 KB
testcase_03	AC	1 ms; 5,376 KB
testcase_04	AC	219 ms; 5,376 KB
testcase_05	AC	2 ms; 5,376 KB
testcase_06	AC	344 ms; 5,376 KB
testcase_07	AC	168 ms; 5,376 KB
testcase_08	AC	332 ms; 5,376 KB
testcase_09	AC	328 ms; 5,376 KB
testcase_10	AC	332 ms; 5,376 KB
testcase_11	AC	321 ms; 5,376 KB
testcase_12	AC	330 ms; 5,376 KB
testcase_13	AC	328 ms; 5,376 KB
testcase_14	AC	148 ms; 5,376 KB
testcase_15	AC	319 ms; 5,376 KB
testcase_16	AC	18 ms; 5,376 KB
testcase_17	AC	8 ms; 5,376 KB
testcase_18	AC	120 ms; 5,376 KB
testcase_19	AC	265 ms; 5,376 KB
testcase_20	AC	100 ms; 5,376 KB
testcase_21	AC	97 ms; 5,376 KB

ソースコード

// 

#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;
}