結果

問題 No.2352 Sharpened Knife in Fall
ユーザー Jeroen Op de Beek
提出日時 2023-06-16 21:34:47
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 277 ms / 3,000 ms
コード長 1,588 bytes
コンパイル時間 1,973 ms
コンパイル使用メモリ 194,512 KB
最終ジャッジ日時 2025-02-14 04:36:37
ジャッジサーバーID
（参考情報） 		judge5 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include "bits/stdc++.h"
using namespace std;
#define all(x) begin(x),end(x)
template<typename A, typename B> ostream& operator<<(ostream &os, const pair<A, B> &p) { return os << '(' << p.first << ", " << p.second << ')'; }
template<typename T_container, typename T = typename enable_if<!is_same<T_container, string>::value, typename T_container::value_type>::type> ostream& operator<<(ostream &os, const T_container &v) { string sep; for (const T &x : v) os << sep << x, sep = " "; return os; }
#define debug(a) cerr << "(" << #a << ": " << a << ")\n";
typedef long long ll;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef pair<int,int> pi;
const int mxN = 1e5+1, oo = 1e9;
int main() {
    ios_base::sync_with_stdio(false);
    cin.tie(NULL);
    double r,k; cin >> r >> k;
    double vol = r*r*acos(-1);
    for(int i=1;i<k+1;++i) {
        double need = vol*i/(k+1);
        double lo = -r, hi = r;
        for(int rep=0;rep<50;++rep) {
            auto mid = (lo+hi)/2;
            // need to 
            auto calc = [&]() {
                if(mid<0) {
                    auto ang = acos(abs(mid)/r)*2;
                    double have = ang*r*r/2-(abs(mid)*sqrt(r*r-mid*mid));
                    return have;
                } else {
                    auto ang = acos(abs(mid)/r)*2;
                    double have = ang*r*r/2-(abs(mid)*sqrt(r*r-mid*mid));
                    return vol-have;
                }
            };
            auto have = calc();

            if(have>=need) hi = mid;
            else lo = mid;

        }
        cout << lo << '\n';
    }
}
0