結果

問題 No.58 イカサマなサイコロ
ユーザー nnasennnasen
提出日時 2017-12-21 18:28:40
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 2 ms / 5,000 ms
コード長 934 bytes
コンパイル時間 2,472 ms
コンパイル使用メモリ 165,088 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-08-22 20:51:57
合計ジャッジ時間 2,520 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge11
このコードへのチャレンジ（β）

テストケース

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

ソースコード

diff # 

#include <bits/stdc++.h>
#define FOR(i,a,b) for(int i = (a); i < (b); ++i)
#define REP(i,n) FOR(i,0,n)
#define SZ(n) (int)(n).size()
#define ALL(n) (n).begin(), (n).end()
#define MOD 1000003
#define INF 100000000
using namespace std;
typedef long long LL;
typedef vector<int> VI;
typedef pair<int, int> PI;


int main() {
	int n, k;
	cin >> n >> k;
	double dp1[101] = {};
	double dp2[101] = {};
	dp1[0] = 1; 
	dp2[0] = 1;
	REP(i, n) {
		for (int j = 60; j >= 0; --j) {
			FOR(k, 1, 7) {
				dp1[j + k] += dp1[j] / 6;
			}
			dp1[j] = 0;
		}
	}
	REP(i, k) {
		for (int j = 60; j >= 0; --j) {
			FOR(k, 4, 7) {
				dp2[j + k] += dp2[j] / 3;
			}
			dp2[j] = 0;
		}
	}
	REP(i, n - k) {
		for (int j = 60; j >= 0; --j) {
			FOR(k, 1, 7) {
				dp2[j + k] += dp2[j] / 6;
			}
			dp2[j] = 0;
		}
	}
	double ans = 0;
	REP(i, 100) {
		REP(j, i) {
			ans += dp1[j] * dp2[i];
		}
	}
	cout << setprecision(6) << fixed << ans << endl;
	return 0;
}
0