結果

問題 No.2352 Sharpened Knife in Fall
コンテスト
ユーザー kumakuma
提出日時 2023-05-13 22:47:00
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 345 ms / 3,000 ms
コード長 1,205 bytes
コンパイル時間 2,314 ms
コンパイル使用メモリ 201,208 KB
最終ジャッジ日時 2025-02-13 00:12:57
ジャッジサーバーID
（参考情報） 		judge3 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
// #include <atcoder/modint>
#define rng(a) a.begin(),a.end()
#define rrng(a) a.rbegin(),a.rend()
#define INF 2000000000000000000
#define ll long long
#define ull unsigned long long
#define ld long double
#define pll pair<ll, ll>
using namespace std;
template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; }
template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return 1; } return 0; }
const double PI = 3.141592653589793238462643383279;

int main() {
  ios::sync_with_stdio(false);
  cin.tie(nullptr);
  ll R, K;
  cin >> R >> K;
  ld S = R * R * PI;
  vector<ld> ans;
  for (ll i = 0; i < K / 2; ++i) {
    ld want_s = (i + 1) * S / (ld)(K + 1);
    ld l = 0, r = R;
    while (l + 0.0000001 <= r) {
      ld m = (l + r) / 2;
      ld half_sita = acosl(m / R);
      ld s = R * R * half_sita - m * sqrtl(R * R - m * m);
      if (s <= want_s) {
        r = m;
      }
      else {
        l = m;
      }
    }
    ans.push_back(r);
    ans.push_back(-r);
  }
  if (K % 2 == 1) {
    ans.push_back(0);
  }
  sort(rng(ans));
  cout << fixed << setprecision(15);
  for (ll i = 0; i < ans.size(); ++i) {
    cout << ans.at(i) << "\n";
  }
}
0