結果

問題 No.2352 Sharpened Knife in Fall
ユーザー norikame
提出日時 2023-06-16 23:19:51
言語 C++23
(gcc 13.3.0 + boost 1.87.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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

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