結果
問題 | No.2352 Sharpened Knife in Fall |
ユーザー |
![]() |
提出日時 | 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 |
ソースコード
// #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; }