結果

問題 No.2971 無理積分
コンテスト
ユーザー haihamabossu
提出日時 2024-11-30 00:18:54
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 30 ms / 1,000 ms
コード長 795 bytes
コンパイル時間 2,400 ms
コンパイル使用メモリ 193,492 KB
最終ジャッジ日時 2025-02-26 09:32:30
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 36
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
#define rep(i, n) for (int i = 0; i < (int)(n); i++)

double f(double a, ll n) {
  double res = 1;
  rep(i, n) res *= a + i;
  return res;
}

void solve() {
  ll n;
  cin >> n;
  if (n == 0) {
    cout << "0.4\n";
    return;
  }
  double ans = 0, c = 1, d = 1;
  rep(i, 50) {
    double now = 1;
    now *= f(-1.0 / 2.0, i);
    now *= f(1.0 / 3.0, i);
    now /= f(4.0 / 3.0, i);
    now *= c;
    now /= d;
    c *= -1.0 / double(n * n * n);
    d *= i + 1;
    ans += now;
  }
  ans *= n;
  ans *= sqrt(double(n));
  cout << fixed << setprecision(6) << ans << '\n';
}

int main() {
  std::cin.tie(nullptr);
  std::ios_base::sync_with_stdio(false);
  int T = 1;
  for (int t = 0; t < T; t++) {
    solve();
  }
  return 0;
}
0