結果
問題 | No.2352 Sharpened Knife in Fall |
ユーザー |
|
提出日時 | 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 |
ソースコード
#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'; } }