結果

問題 No.76 回数の期待値で練習
ユーザー simansiman
提出日時 2021-08-23 07:28:59
言語 C++17(clang)
(17.0.6 + boost 1.83.0)
結果
AC  
実行時間 19 ms / 5,000 ms
コード長 1,029 bytes
コンパイル時間 3,918 ms
コンパイル使用メモリ 142,928 KB
実行使用メモリ 11,008 KB
最終ジャッジ日時 2024-11-07 08:01:41
合計ジャッジ時間 1,848 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 19 ms
11,008 KB
testcase_01 AC 19 ms
11,008 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cassert>
#include <cmath>
#include <algorithm>
#include <iostream>
#include <iomanip>
#include <limits.h>
#include <map>
#include <queue>
#include <set>
#include <string.h>
#include <vector>

using namespace std;
typedef long long ll;

int main() {
  int T;
  cin >> T;
  vector<int> N(T);

  for (int i = 0; i < T; ++i) {
    cin >> N[i];
  }

  vector<double> P(7, 0.0);
  vector<double> dp(1000010, 0.0);
  dp[1] = 1.0;
  dp[2] = 1.0833333333333333;
  dp[3] = 1.2569444444444444;
  dp[4] = 1.5353009259259260;
  dp[5] = 1.6915991512345676;
  dp[6] = 2.0513639724794235;

  for (int n = 1; n <= 5; ++n) {
    double p = dp[n + 1] - 1;

    for (int i = 1; i <= n - 1; ++i) {
      p -= dp[n - i + 1] * P[i];
    }

    P[n] = p;
  }

  P[6] = 1.0 - (P[1] + P[2] + P[3] + P[4] + P[5]);

  for (int n = 7; n <= 1000000; ++n) {
    dp[n] = 1.0;

    for (int i = 1; i <= 6; ++i) {
      dp[n] += P[i] * dp[n - i];
    }
  }

  for (int n : N) {
    cout << fixed << setprecision(15) << dp[n] << endl;
  }

  return 0;
}
0