結果
問題 | No.58 イカサマなサイコロ |
ユーザー | Mister |
提出日時 | 2020-04-12 03:30:23 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 2 ms / 5,000 ms |
コード長 | 1,254 bytes |
コンパイル時間 | 836 ms |
コンパイル使用メモリ | 81,224 KB |
実行使用メモリ | 4,348 KB |
最終ジャッジ日時 | 2023-10-21 14:44:30 |
合計ジャッジ時間 | 1,397 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge13 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
4,348 KB |
testcase_01 | AC | 2 ms
4,348 KB |
testcase_02 | AC | 2 ms
4,348 KB |
testcase_03 | AC | 2 ms
4,348 KB |
testcase_04 | AC | 2 ms
4,348 KB |
testcase_05 | AC | 2 ms
4,348 KB |
testcase_06 | AC | 2 ms
4,348 KB |
testcase_07 | AC | 1 ms
4,348 KB |
testcase_08 | AC | 2 ms
4,348 KB |
testcase_09 | AC | 2 ms
4,348 KB |
ソースコード
#include <iostream> #include <iomanip> #include <vector> using ldouble = long double; std::vector<ldouble> calc(int n, int k) { int m = n * 6; std::vector<ldouble> dp(m + 1, 0); auto ndp = dp; dp[0] = 1; for (int i = 0; i < n; ++i) { std::fill(ndp.begin(), ndp.end(), 0); for (int d = 0; d <= m; ++d) { if (i < k) { for (int f = 4; f <= 6; ++f) { if (d + f > m) continue; ndp[d + f] += dp[d] / 3; } } else { for (int f = 1; f <= 6; ++f) { if (d + f > m) continue; ndp[d + f] += dp[d] / 6; } } } std::swap(dp, ndp); } return dp; } void solve() { int n, k; std::cin >> n >> k; int m = n * 6; auto adp = calc(n, 0), bdp = calc(n, k); ldouble ans = 0; for (int b = 0; b <= m; ++b) { for (int a = 0; a < b; ++a) { ans += adp[a] * bdp[b]; } } std::cout << std::fixed << std::setprecision(10) << ans << std::endl; } int main() { std::cin.tie(nullptr); std::ios::sync_with_stdio(false); solve(); return 0; }