結果

問題 No.2352 Sharpened Knife in Fall
ユーザー norikame
提出日時 2023-06-16 23:19:51
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 344 ms / 3,000 ms
コード長 1,501 bytes
コンパイル時間 5,471 ms
コンパイル使用メモリ 313,316 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-24 16:29:54
合計ジャッジ時間 15,002 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

// 
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
using ll = long long;
#define rep(i, n) for (int i=0; i<(int)(n); ++(i))
#define rep3(i, m, n) for (int i=(m); (i)<(int)(n); ++(i))
#define repr(i, n) for (int i=(int)(n)-1; (i)>=0; --(i))
#define rep3r(i, m, n) for (int i=(int)(n)-1; (i)>=(int)(m); --(i))
#define all(x) (x).begin(), (x).end()
const double PI = acos(-1);
int main() {
    int r, k;
    cin >> r >> k;
    if (k%2 == 1) {
        vector<double> res;
        double tar = ((double)r * r * PI / 2.0) / (k/2+1);
        res.push_back(0.0);
        rep3(i, 1, k/2+1) {
            double li = 0.0, ri = r;
            while (ri-li > 1e-5) {
                double ci = (li + ri) / 2.0, th = asin(ci/r);
                double area = ((double)r * r * PI / 2.0) * (th / PI) * 2.0 + r * cos(th) * 2.0 * ci / 2.0;
                if (area < tar*i) li = ci;
                else ri = ci;
            }
            res.push_back(ri);
            res.push_back(-ri);
        }
        sort(all(res));
        rep(i, k) cout << fixed << setprecision(10) << res[i] << endl;
    }
    else {
        vector<double> res;
        double tar = ((double)r * r * PI / 2.0) / (k+1);
        for (int i=1; i<k; i+=2) {
            double li = 0.0, ri = r;
            while (ri-li > 1e-5) {
                double ci = (li + ri) / 2.0, th = asin(ci/r);
                double area = ((double)r * r * PI / 2.0) * (th / PI) * 2.0 + r * cos(th) * 2.0 * ci / 2.0;
                if (area < tar*i) li = ci;
                else ri = ci;
            }
            res.push_back(ri);
            res.push_back(-ri);
        }
        sort(all(res));
        rep(i, k) cout << fixed << setprecision(10) << res[i] << endl;
    }
    return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0