結果

問題 No.2352 Sharpened Knife in Fall
ユーザー hari64hari64
提出日時 2023-06-16 22:13:55
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 592 ms / 3,000 ms
コード長 3,320 bytes
コンパイル時間 2,055 ms
コンパイル使用メモリ 203,892 KB
実行使用メモリ 5,680 KB
最終ジャッジ日時 2023-09-06 20:14:57
合計ジャッジ時間 16,244 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,384 KB
testcase_01 AC 1 ms
4,384 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 1 ms
4,384 KB
testcase_04 AC 548 ms
5,616 KB
testcase_05 AC 2 ms
4,380 KB
testcase_06 AC 575 ms
5,512 KB
testcase_07 AC 287 ms
4,568 KB
testcase_08 AC 592 ms
5,492 KB
testcase_09 AC 582 ms
5,512 KB
testcase_10 AC 569 ms
5,680 KB
testcase_11 AC 563 ms
5,564 KB
testcase_12 AC 566 ms
5,560 KB
testcase_13 AC 565 ms
5,552 KB
testcase_14 AC 273 ms
4,484 KB
testcase_15 AC 539 ms
5,544 KB
testcase_16 AC 27 ms
4,388 KB
testcase_17 AC 15 ms
4,384 KB
testcase_18 AC 201 ms
4,460 KB
testcase_19 AC 457 ms
5,516 KB
testcase_20 AC 172 ms
4,384 KB
testcase_21 AC 161 ms
4,388 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef __int128_t lll;
typedef long double ld;
typedef pair<ll, ll> pl;
typedef tuple<ll, ll, ll> tt;
typedef vector<ll> vl;
typedef vector<vl> vvl;
typedef vector<vvl> vvvl;
typedef vector<double> vd;
typedef vector<vd> vvd;
typedef vector<vvd> vvvd;
typedef vector<bool> vb;
typedef vector<vb> vvb;
typedef vector<vvb> vvvb;
typedef vector<char> vc;
typedef vector<vc> vvc;
typedef vector<vvc> vvvc;
typedef vector<pl> vp;
typedef vector<tt> vt;
typedef vector<string> vs;

const ll INF = 1000000010;
const ll INFL = INF * INF;
const double MYPI = acos(-1);
const ld epsilon = 1e-10;

#define ovl(a, b, c, d, e, ...) e
#define pb push_back
#define eb emplace_back
#define MP make_pair
#define DEBUG(...) DEBUG_(#__VA_ARGS__, __VA_ARGS__)
#define REP1(i, n) for (int i = 0; i < (n); i++)
#define REP2(i, l, r) for (int i = (l); i < (r); i++)
#define REP3(i, l, r, d) for (int i = (l); i < (r); i += (d))
#define REP(...) ovl(__VA_ARGS__, REP3, REP2, REP1)(__VA_ARGS__)
#define RREP1(i, n) for (int i = (n)-1; i >= 0; i--)
#define RREP2(i, l, r) for (int i = (r)-1; i >= (l); i--)
#define RREP3(i, l, r, d) for (int i = (r)-1; i >= (l); i -= (d))
#define RREP(...) ovl(__VA_ARGS__, RREP3, RREP2, RREP1)(__VA_ARGS__)
#define EACH1(e, a) for (auto &e : a)
#define EACH2(x, y, a) for (auto &[x, y] : a)
#define EACH3(x, y, z, a) for (auto &[x, y, z] : a)
#define EACH(...) ovl(__VA_ARGS__, EACH3, EACH2, EACH1)(__VA_ARGS__)
#define ALL(x) begin(x), end(x)
#define RALL(x) (x).rbegin(), (x).rend()
#define sz(x) (ll) x.size()
#define LB(a, x) (lower_bound(ALL(a), x) - a.begin())
#define UB(a, x) (upper_bound(ALL(a), x) - a.begin())
#define FLG(x, i) (((x) >> (i)) & 1)
#define CNTBIT(x) __builtin_popcountll(x)
#define TOPBIT(t) (t == 0 ? -1 : 63 - __builtin_clzll(t))
#define IN(x, a, b) ((a) <= (x) && (x) < (b))

template <typename T, typename... Ts>
void OUT(const T &a, const Ts &...b) {
    cout << a;
    (cout << ... << (cout << " ", b));
    cout << endl;
}

template <typename T>
void DEBUG_(string_view name, const T &a) {
    cout << name << ": " << a << endl;
}

template <typename T>
void DEBUG_(string_view name, const vector<T> &a) {
    cout << name << ": ";
    REP(i, sz(a)) cout << a[i] << " ";
    cout << endl;
}

template <typename T1, typename T2>
bool chmax(T1 &x, const T2 &y) {
    return x < y ? x = y, 1 : 0;
}
template <typename T1, typename T2>
bool chmin(T1 &x, const T2 &y) {
    return x > y ? x = y, 1 : 0;
}

int main() {
    cin.tie(0);
    ios::sync_with_stdio(false);

    long double PI = acos(-1);

    int R, K;
    cin >> R >> K;

    auto is_ok = [&](long double theta, int k) -> bool {
        return (PI * (k + 1) / (K + 1)) > (theta - sin(theta) * cos(theta));
    };

    vector<long double> ans;
    if (K % 2 == 1) ans.pb(0);

    REP(k, K / 2) {
        long double ok, ng;
        ok = 0.0;
        ng = PI / 2;
        while (abs(ok - ng) > 1e-10) {
            long double mid = (ok + ng) / 2;
            if (is_ok(mid, k)) {
                ok = mid;
            } else {
                ng = mid;
            }
        }
        ans.pb(R * cos(ok));
        ans.pb(-R * cos(ok));
    }

    sort(ALL(ans));

    EACH1(a, ans) { cout << fixed << setprecision(16) << a << endl; }

    return 0;
}
0