結果

問題 No.58 イカサマなサイコロ
コンテスト
ユーザー 37kt_
提出日時 2015-10-30 21:39:06
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 4 ms / 5,000 ms
コード長 1,268 bytes
コンパイル時間 1,944 ms
コンパイル使用メモリ 158,728 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-13 05:32:40
合計ジャッジ時間 2,330 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 10
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

#define lp(i, a, b) for (int i = (int)(a); i < (int)(b); i++)
#define rep(i, n) lp(i, 0, n)
#define rlp(i, a, b) for (int i = (int)(b) - 1; i >= (int)(a); i--)
#define rrep(i, n) rlp(i, 0, n)

#define mp make_pair
#define mt make_tuple
#define pb push_back
#define eb emplace_back
#define all(c) begin(c), end(c)

typedef long long ll;
typedef unsigned long long ull;

template<typename T>
inline bool chmin(T &a, T b){
	if (a <= b) return false;
	a = b; return true;
}

template<typename T>
inline bool chmax(T &a, T b){
	if (a >= b) return false;
	a = b; return true;
}

template<typename T>
inline void unique(T &c){
	sort(all(c));
	c.erase(unique(all(c)), end(c));
}

template<typename T>
inline T sq(T a){
	return a * a;
}

int n, k;
double dp[10][100][100];

double calc(int i, int a, int b)
{
	if (i == n) return a > b ? 1 : 0;
	if (dp[i][a][b] != -1.0) return dp[i][a][b];
	double res = 0;
	if (i < k){
		lp(i1, 4, 7) lp(i2, 1, 7){
			res += 2 * calc(i + 1, a + i1, b + i2);
		}
	}
	else {
		lp(i1, 1, 7) lp(i2, 1, 7){
			res += calc(i + 1, a + i1, b + i2);
		}
	}
	return dp[i][a][b] = res / 36;
}

int main()
{
	cin >> n >> k;
	fill_n(**dp, 100000, -1);
	double res = calc(0, 0, 0);
	printf("%.20f\n", res);
}
0