結果
| 問題 |
No.2352 Sharpened Knife in Fall
|
| コンテスト | |
| ユーザー |
srjywrdnprkt
|
| 提出日時 | 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;
}
srjywrdnprkt