結果

問題 No.2352 Sharpened Knife in Fall
ユーザー t98slidert98slider
提出日時 2023-06-16 21:41:32
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 266 ms / 3,000 ms
コード長 1,284 bytes
コンパイル時間 1,823 ms
コンパイル使用メモリ 174,244 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-24 13:29:44
合計ジャッジ時間 11,664 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using ll = long long;

template<class T> istream& operator >> (istream& is, vector<T>& vec) {
    for(T& x : vec) is >> x;
    return is;
}

template<class T> ostream& operator << (ostream& os, const vector<T>& vec) {
    if(vec.empty()) return os;
    os << vec[0];
    for(auto it = vec.begin(); ++it != vec.end(); ) os << ' ' << *it;
    return os;
}

int main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
    ll r, k;
    cin >> r >> k;
    double S = r * r * M_PI / (k + 1);
    //cerr << S << '\n';
    //cerr << S * (k + 1) << " " << r * r * M_PI << '\n';
    vector<double> ans;
    if(k & 1){
        ans.emplace_back(0);
    }
    ll k2 = k / 2;
    double S2 = 0;
    for(int i = 0; i < k2; i++){
        double ng = 0, ok = M_PI, mid;
        for(int j = 0; j < 60; j++){
            mid = (ng + ok) / 2;
            double tmp = r * r * mid / 2;
            tmp -= r * r * sin(mid / 2) * cos(mid / 2);
            if(tmp >= S * (i + 1)) ok = mid;
            else ng = mid;
        }
        ans.emplace_back(r * cos(ok / 2));
        ans.emplace_back(-r * cos(ok / 2));
    }
    int c = ans.size();
    sort(ans.begin(), ans.end());
    cout << fixed << setprecision(15);
    for(auto &&v : ans) cout << v << '\n';
}
0