結果
問題 | No.2352 Sharpened Knife in Fall |
ユーザー |
![]() |
提出日時 | 2023-06-16 22:58:38 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 311 ms / 3,000 ms |
コード長 | 1,418 bytes |
コンパイル時間 | 3,090 ms |
コンパイル使用メモリ | 258,388 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-06-24 15:58:45 |
合計ジャッジ時間 | 15,071 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 19 |
ソースコード
#include <bits/stdc++.h> using namespace std; #define FOR(i,m,n) for(int i=(m);i<(n);++i) #define REP(i,n) FOR(i,0,n) #define ALL(v) (v).begin(),(v).end() using ll = long long; constexpr int INF = 0x3f3f3f3f; constexpr long long LINF = 0x3f3f3f3f3f3f3f3fLL; constexpr double EPS = 1e-8; constexpr int MOD = 998244353; // constexpr int MOD = 1000000007; constexpr int DY4[]{1, 0, -1, 0}, DX4[]{0, -1, 0, 1}; constexpr int DY8[]{1, 1, 0, -1, -1, -1, 0, 1}; constexpr int DX8[]{0, -1, -1, -1, 0, 1, 1, 1}; template <typename T, typename U> inline bool chmax(T& a, U b) { return a < b ? (a = b, true) : false; } template <typename T, typename U> inline bool chmin(T& a, U b) { return a > b ? (a = b, true) : false; } struct IOSetup { IOSetup() { std::cin.tie(nullptr); std::ios_base::sync_with_stdio(false); std::cout << fixed << setprecision(20); } } iosetup; int main() { int r, k; cin >> r >> k; const double area = numbers::pi / (k + 1); vector<double> ans(k, 0); REP(i, k / 2) { const double a = area * (i + 1); double lb = (ans.empty() ? 0 : ans.back()), ub = numbers::pi / 2; REP(_, 80) { const double mid = midpoint(lb, ub); (mid - sin(mid * 2) / 2 > a ? ub : lb) = mid; } ans[i] = lb; } REP(i, k / 2) ans[i] = -cos(ans[i]) * r; FOR(i, k / 2, k / 2 * 2) ans[i] = -ans[i - k / 2]; ranges::sort(ans); REP(i, k) cout << ans[i] << '\n'; return 0; }