結果

問題 No.58 イカサマなサイコロ
ユーザー SSRSSSRS
提出日時 2020-09-02 07:53:24
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 2 ms / 5,000 ms
コード長 930 bytes
コンパイル時間 845 ms
コンパイル使用メモリ 83,232 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-05-01 10:23:53
合計ジャッジ時間 1,322 ms
ジャッジサーバーID
(参考情報)
judge4 / judge1
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
6,940 KB
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 2 ms
6,940 KB
testcase_04 AC 2 ms
6,940 KB
testcase_05 AC 2 ms
6,940 KB
testcase_06 AC 2 ms
6,944 KB
testcase_07 AC 2 ms
6,940 KB
testcase_08 AC 2 ms
6,944 KB
testcase_09 AC 2 ms
6,944 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <iostream>
#include <iomanip>
#include <vector>
using namespace std;
int main(){
	cout << fixed << setprecision(5);
	int N;
	cin >> N;
	int K;
	cin >> K;
	vector<vector<double>> dp(N + 1, vector<double>(N * 6 + 1, 0));
	dp[0][0] = 1;
	for (int i = 0; i < N; i++){
		for (int j = 0; j <= i * 6; j++){
			for (int k = 1; k <= 6; k++){
				dp[i + 1][j + k] += dp[i][j] / 6;
			}
		}
	}
	vector<vector<double>> dp2(N + 1, vector<double>(N * 6 + 1, 0));
	dp2[0][0] = 1;
	for (int i = 0; i < K; i++){
		for (int j = 0; j <= i * 6; j++){
			for (int k = 4; k <= 6; k++){
				dp2[i + 1][j + k] += dp2[i][j] / 3;
			}
		}
	}
	for (int i = K; i < N; i++){
		for (int j = 0; j <= i * 6; j++){
			for (int k = 1; k <= 6; k++){
				dp2[i + 1][j + k] += dp2[i][j] / 6;
			}
		}
	}
	double ans = 0;
	for (int i = 0; i <= N * 6; i++){
		for (int j = i + 1; j <= N * 6; j++){
			ans += dp[N][i] * dp2[N][j];
		}
	}
	cout << ans << endl;
}
0