結果

問題 No.2352 Sharpened Knife in Fall
ユーザー みしあみしあ
提出日時 2023-07-27 18:36:29
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 396 ms / 3,000 ms
コード長 3,796 bytes
コンパイル時間 3,662 ms
コンパイル使用メモリ 263,632 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-04 18:06:17
合計ジャッジ時間 15,079 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
5,248 KB
testcase_01 AC 1 ms
5,248 KB
testcase_02 AC 1 ms
5,248 KB
testcase_03 AC 2 ms
5,248 KB
testcase_04 AC 348 ms
5,248 KB
testcase_05 AC 2 ms
5,248 KB
testcase_06 AC 395 ms
5,248 KB
testcase_07 AC 197 ms
5,248 KB
testcase_08 AC 391 ms
5,248 KB
testcase_09 AC 392 ms
5,248 KB
testcase_10 AC 394 ms
5,248 KB
testcase_11 AC 393 ms
5,248 KB
testcase_12 AC 387 ms
5,248 KB
testcase_13 AC 396 ms
5,248 KB
testcase_14 AC 188 ms
5,248 KB
testcase_15 AC 372 ms
5,248 KB
testcase_16 AC 21 ms
5,248 KB
testcase_17 AC 11 ms
5,248 KB
testcase_18 AC 141 ms
5,248 KB
testcase_19 AC 306 ms
5,248 KB
testcase_20 AC 122 ms
5,248 KB
testcase_21 AC 116 ms
5,248 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#pragma region Macros
#include <bits/stdc++.h>
#include <atcoder/all>
#if defined(LOCAL) || defined(_DEBUG)
#include "debug.hpp"
#else
#define O(...)
#define START()
#define STOP()
#define MEMORY()
#endif
using namespace std;
#define REP(i, n) for(int i=0, i##_len=(n); i<i##_len; ++i)
#define REPR(i, n) for(int i=(n); i>=0; --i)
#define FOR(i, n, m) for(int i=(m), i##_len=(n); i<i##_len; ++i)
#define EACH(i, v) for(const auto& i : v)
#define ALL(x) (x).begin(),(x).end()
#define ALLR(x) (x).rbegin(),(x).rend()
template<class T, class U>bool chmax(T &a, const U &b) { if (a<(T)b) { a=(T)b; return 1; } return 0; }
template<class T, class U>bool chmin(T &a, const U &b) { if (b<(T)a) { a=(T)b; return 1; } return 0; }
#define vec vector
#define umap unordered_map
#define uset unordered_set
using ll = long long;
using ld = long double;
using P = pair<ll, ll>;
using Tup = tuple<ll, ll, ll>;
using vl = vec<ll>;
#define fi first
#define se second
#define el endl
constexpr ll INF = numeric_limits<ll>::max()/2-1;

template<class T> istream &operator>>(istream &stream, vec<T>& o){REP(i, o.size())stream >> o[i];return stream;}
template<class T, class U> istream &operator>>(istream &stream, pair<T, U>& o){cin >> o.fi >> o.se; return stream;}
namespace myinput {
    void input() {}
    template<class T, class... Ts> void input(T&& o, Ts&&... args){cin >> o;input(forward<Ts>(args)...);}
    void assign_vl(size_t siz) {};
    template<class T, class... Ts> void assign_vl(size_t siz, T&& o, Ts&&... args){o.resize(siz);assign_vl(siz, forward<Ts>(args)...);}
    void input_horizon_sub(size_t index) {};
    template<class T, class... Ts> void input_horizon_sub(size_t index, T&& o, Ts&&... args) {cin>>o[index];input_horizon_sub(index, forward<Ts>(args)...);}
    template<class... Ts> void input_horizon(size_t siz, Ts&&... args){REP(i, siz) input_horizon_sub(i, forward<Ts>(args)...);}
}
#define _I(T, ...) ;myinput::input(__VA_ARGS__);
#define I(T, ...) ;T __VA_ARGS__;_I(T, __VA_ARGS__);
#define _Iv(T, siz, ...) ;myinput::assign_vl(siz, __VA_ARGS__);myinput::input(__VA_ARGS__);
#define Iv(T, siz, ...) ;vec<T> __VA_ARGS__;_Iv(T, siz, __VA_ARGS__);
#define _Ih(T, siz, ...) ;myinput::assign_vl(siz, __VA_ARGS__);myinput::input_horizon(siz, __VA_ARGS__);
#define Ih(T, siz, ...) ;vec<T> __VA_ARGS__;_Ih(T, siz, __VA_ARGS__);
#pragma endregion

void Main();

int main(){
    std::cin.tie(nullptr);
    std::cout << std::fixed << std::setprecision(15);
    Main();
    MEMORY();
    return 0;
}

// x∈[l, r] | f(x) = true となる最大のxを返す
template <class Func>
ll binarySearch(ll l, ll r, const Func &f){
    r++;
    while(l < r + 1){
        const ll m = (l<0)==(r<0) ? l+(r-l)/2 : (l<-r) ? -((-l-r-1)/2+1) : (l+r)/2;
        if(f(m)) l = m+1;
        else     r = m-1;
    }
    return r;
}

template <typename T, typename U>
U __typechange(T n) {
    if(n >= 0) return *(U*)(&n);
    n = -n;
    return -*(U*)(&n);
}
template <class Func>
double binarySearchd(double l, double r, const Func &f) {
    return __typechange<ll, double>(
        binarySearch(__typechange<double, ll>(l), __typechange<double, ll>(r), 
            [&](ll mi) { return f(__typechange<ll, double>(mi)); }
        )
    );
}

void Main(){
    /*
    REP(i, 20) {
        ll r = 1;
        double h = r-(i/10.0-1.0)/100.0;
        double theta = 2*acos(1-h/r);
        double s = theta*r*r/2 - (r-h) * sqrt(h*(2*r-h));
        O(s);
    }*/
    I(ll, r, k);
    double piece = M_PI*r*r / (k+1);
    O(piece);
    REP(i, k) {
        double ret = binarySearchd(-r, r, [&](double d) {
            double h = r-d;
            double theta = 2*acos(1-h/r);
            double s = theta*r*r/2 - (r-h) * sqrt(h*(2*r-h));
            return s > piece*(k-i);
        });
        cout << ret << el;
    }
}
0