結果

問題 No.58 イカサマなサイコロ
コンテスト
ユーザー nnasen
提出日時 2017-12-21 18:28:40
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 5,000 ms
コード長 934 bytes
コンパイル時間 1,500 ms
コンパイル使用メモリ 168,176 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-12-17 22:38:58
合計ジャッジ時間 2,184 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 10
権限があれば一括ダウンロードができます

ソースコード

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