結果

問題 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
権限があれば一括ダウンロードができます

ソースコード

diff #

#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;
}
0