結果
問題 | No.2352 Sharpened Knife in Fall |
ユーザー |
![]() |
提出日時 | 2023-06-16 22:33:53 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 848 ms / 3,000 ms |
コード長 | 1,679 bytes |
コンパイル時間 | 1,147 ms |
コンパイル使用メモリ | 112,228 KB |
最終ジャッジ日時 | 2025-02-14 06:36:14 |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 19 |
ソースコード
#include <iostream> #include <vector> #include <cmath> #include <map> #include <set> #include <iomanip> #include <queue> #include <algorithm> #include <numeric> #include <deque> #include <complex> #include <cassert> using namespace std; using ll = long long; int main(){ long double R, pi=M_PI, S, T, U; int K; cin >> R >> K; S = pi*R*R / (long double) (K+1); vector<long double> v; long double l, r, c; if (K % 2 == 0){ K -= 2; T = 0; l=0, r=pi/2; for (int j=0; j<100; j++){ c = (l+r)/2; U = (c + sin(c)*cos(c)) * R * R; if (U*2 <= S) l=c; else r=c; } T = (l + sin(l)*cos(l)) * R * R; v.push_back(R*sin(l)); v.push_back(-R*sin(l)); for (int i=0; i<K/2; i++){ l=0, r=pi/2; for (int j=0; j<100; j++){ c = (l+r)/2; U = (c + sin(c)*cos(c)) * R * R - T; if (U <= S) l=c; else r=c; } T = (l + sin(l)*cos(l)) * R * R; v.push_back(R*sin(l)); v.push_back(-R*sin(l)); } } else{ v.push_back(0); T = 0; for (int i=0; i<K/2; i++){ l=0, r=pi/2; for (int j=0; j<100; j++){ c = (l+r)/2; U = (c + sin(c)*cos(c)) * R * R - T; if (U <= S) l=c; else r=c; } T = (l + sin(l)*cos(l)) * R * R; v.push_back(R*sin(l)); v.push_back(-R*sin(l)); } } sort(v.begin(), v.end()); for (auto x : v) cout << setprecision(18) << x << endl; return 0; }