結果

問題 No.2352 Sharpened Knife in Fall
ユーザー coindarw
提出日時 2023-06-16 21:44:59
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 278 ms / 3,000 ms
コード長 2,882 bytes
コンパイル時間 4,095 ms
コンパイル使用メモリ 255,204 KB
最終ジャッジ日時 2025-02-14 05:02:42
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

#include <atcoder/all>
using ll = long long;
#define rep(i, n) for (int i = 0, i##_len = (n); i < i##_len; ++i)
#define reps(i, n) for (int i = 1, i##_len = (n); i <= i##_len; ++i)
#define rrep(i, n) for (int i = ((int)(n)-1); i >= 0; --i)
#define rreps(i, n) for (int i = ((int)(n)); i > 0; --i)
#define rep2(i, s, n) for (int i = (s); i < (int)(n); ++i)
#define repc2(i, s, n) for (int i = (s); i <= (int)(n); ++i)
constexpr int inf = 2000'000'000;
constexpr ll linf = 4000000000000000000ll;
constexpr ll M7 = 1000000007ll;
constexpr ll M09 = 1000000009ll;
constexpr ll M9 = 998244353ll;
#define all(v) begin(v), end(v)
#define rall(v) rbegin(v), rend(v)
using namespace std;
using namespace atcoder;
template <typename T>
inline ostream& operator<<(ostream& os, vector<T>& v) {
    for (auto& e : v) os << e << " ";
    return os;
}
template <typename T, typename U>
std::ostream& operator<<(std::ostream& os, const std::pair<T, U>& p) noexcept {
    return os << "(" << p.first << ", " << p.second << ")";
}
#ifdef ONLINE_JUDGE
#define debug(...)
#else
#define debug(...) cerr << "<" << #__VA_ARGS__ << ">: ", debug_out(__VA_ARGS__)
template <typename T>
void debug_out(T t) {
    cerr << t << endl;
}
template <typename T, typename... Args>
void debug_out(T t, Args... args) {
    cerr << t << ", ";
    debug_out(args...);
}
#endif

int main() {
    ios_base::sync_with_stdio(false);
    cin.tie(NULL);
    int r, k;
    cin >> r >> k;
    auto f = [&](double t) { return (t - sin(t)) / 2; };
    vector<double> v;
    if (k % 2 == 0) {
        reps(i, k / 2) {
            auto cmp = [&](double t) { return f(t) < i * M_PI / (k + 1); };
            auto bsearch = [&]() {
                double ok = 0, ng = M_PI;
                rep(_, 100) {
                    double mid = (ng + ok) / 2;
                    if (cmp(mid))
                        ok = mid;
                    else
                        ng = mid;
                }
                return ok;
            };
            double t = bsearch();
            v.push_back(cos(t / 2) * r);
            v.push_back(-cos(t / 2) * r);
        }
    } else {
        v.push_back(0);
        reps(i, k / 2) {
            auto cmp = [&](double t) { return f(t) < i * M_PI / (k + 1); };
            auto bsearch = [&]() {
                double ok = 0, ng = M_PI;
                rep(_, 100) {
                    double mid = (ng + ok) / 2;
                    if (cmp(mid))
                        ok = mid;
                    else
                        ng = mid;
                }
                return ok;
            };
            double t = bsearch();
            v.push_back(cos(t / 2) * r);
            v.push_back(-cos(t / 2) * r);
        }
    }
    sort(all(v));
    cout << fixed << setprecision(10);
    for (double x : v) {
        cout << x << "\n";
    }
    return 0;
}
0