結果
問題 | No.76 回数の期待値で練習 |
ユーザー | siman |
提出日時 | 2021-08-23 07:27:28 |
言語 | C++17(clang) (17.0.6 + boost 1.83.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,092 bytes |
コンパイル時間 | 2,665 ms |
コンパイル使用メモリ | 142,456 KB |
実行使用メモリ | 11,136 KB |
最終ジャッジ日時 | 2024-11-07 07:59:45 |
合計ジャッジ時間 | 3,176 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ソースコード
#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 i = 1; i <= 6; ++i) { cout << P[i] << endl; } 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(10) << dp[n] << endl; } return 0; }