結果
問題 | No.2352 Sharpened Knife in Fall |
ユーザー |
![]() |
提出日時 | 2023-06-16 21:40:28 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 536 ms / 3,000 ms |
コード長 | 1,245 bytes |
コンパイル時間 | 4,016 ms |
コンパイル使用メモリ | 257,752 KB |
最終ジャッジ日時 | 2025-02-14 04:54:06 |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 19 |
ソースコード
#include <stdio.h> #include <atcoder/all> #include <bits/stdc++.h> using namespace std; using namespace atcoder; using mint = modint998244353; #define rep(i,n) for (int i = 0; i < (n); ++i) #define Inf32 1000000001 #define Inf64 4000000000000000001 long long R,K; long double get(long double h){ long double r = R; return r * r * acosl(1.0 - h/r) - (r-h) * sqrtl(2.0 * r * h - h*h); } int main(){ long double pi = acosl(-1.0); cin>>R>>K; long double S = R; S *= R; S *= pi; S /= K+1; vector<long double> ans; long double last = 0.0; rep(i,1000000){ long double ok = last,ng = 2.0*R; rep(_,50){ long double mid = (ok+ng)/2.0; if(get(mid)<=S*(i+1))ok = mid; else ng = mid; } if(ok >= R-1e-10)break; ans.push_back(-R + ok); last = ok; } int t = ans.size(); rep(i,t)ans.push_back(-ans[i]); if(K%2==1)ans.push_back(0.0); sort(ans.begin(),ans.end()); vector<long double> ta; rep(i,ans.size()){ if(ta.size()>0 && abs(ans[i]-ta.back())<1e-11)continue; ta.push_back(ans[i]); } swap(ans,ta); //cout<<ans.size()<<endl; //cout<<fixed<<setprecision(10)<<ans[K/2]<<','<<ans[K/2+1]<<endl; rep(i,ans.size()){ cout<<fixed<<setprecision(10)<<ans[i]<<endl; } //cout<<ans.size()<<endl; return 0; }