結果

問題 No.2352 Sharpened Knife in Fall
ユーザー srjywrdnprktsrjywrdnprkt
提出日時 2023-06-16 22:33:53
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 845 ms / 3,000 ms
コード長 1,679 bytes
コンパイル時間 1,121 ms
コンパイル使用メモリ 114,672 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-24 15:13:02
合計ジャッジ時間 18,094 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
6,944 KB
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 2 ms
6,944 KB
testcase_04 AC 763 ms
6,944 KB
testcase_05 AC 2 ms
6,940 KB
testcase_06 AC 830 ms
6,940 KB
testcase_07 AC 421 ms
6,940 KB
testcase_08 AC 845 ms
6,940 KB
testcase_09 AC 844 ms
6,940 KB
testcase_10 AC 833 ms
6,940 KB
testcase_11 AC 837 ms
6,940 KB
testcase_12 AC 836 ms
6,944 KB
testcase_13 AC 840 ms
6,940 KB
testcase_14 AC 401 ms
6,940 KB
testcase_15 AC 810 ms
6,940 KB
testcase_16 AC 43 ms
6,944 KB
testcase_17 AC 20 ms
6,944 KB
testcase_18 AC 299 ms
6,944 KB
testcase_19 AC 653 ms
6,940 KB
testcase_20 AC 258 ms
6,940 KB
testcase_21 AC 247 ms
6,940 KB
権限があれば一括ダウンロードができます

ソースコード

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