結果

問題 No.58 イカサマなサイコロ
ユーザー SSRS
提出日時 2020-09-02 07:51:06
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 927 bytes
コンパイル時間 756 ms
コンパイル使用メモリ 81,736 KB
最終ジャッジ日時 2025-01-14 03:31:57
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 2 WA * 8
権限があれば一括ダウンロードができます

ソースコード

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 = 0; 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 = 0; 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