結果

問題 No.58 イカサマなサイコロ
ユーザー IL_mstaIL_msta
提出日時 2015-06-29 09:49:24
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 2 ms / 5,000 ms
コード長 1,497 bytes
コンパイル時間 684 ms
コンパイル使用メモリ 84,620 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-09-22 04:21:54
合計ジャッジ時間 1,436 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

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

ソースコード

diff # 

#define _USE_MATH_DEFINES

#include <iostream>
#include <iomanip>

#include <algorithm>
#include <cmath>

#include <string>
#include <list>
#include <queue>
#include <vector>
#include <complex>
#include <set>

/////////
#define REP(i, x, n) for(int i = x; i < n; i++)
#define rep(i,n) REP(i,0,n)
#define P(p) cout<<(p)<<endl;
/////////
typedef long long LL;
typedef long double LD;
/////////
using namespace::std;
/////////

int main(void){
    std::cin.tie(0); 
    std::ios::sync_with_stdio(false);
    std::cout << std::fixed;//
    //cout << setprecision(6);//
	
	int N,K;//N[1,10]
	cin>>N>>K;
	////
	LL dpA[61];
	LL dpB[61];
	LL dpT[61];
	rep(i,61){
		dpA[i] = 0;
		dpB[i] = 0;
		dpT[i] = 0;
	}
	dpA[0]=1;
	int meA[6]={1,2,3,4,5,6};
	for(int i=0;i<N;++i){//N個サイコロ
		for(int d=0;d<61;++d){
			if(dpA[d]>0){
				for(int k=0;k<6;++k){
					dpT[ d + meA[k] ] += dpA[d];
				}
			}
		}
		rep(zz,61){
			dpA[zz] = dpT[zz];
			dpT[zz] = 0;
		}
	}
	///////
	dpB[0]=1;
	int meB[6]={4,4,5,5,6,6};
	for(int i=0;i<N;++i){//N個サイコロ
		for(int d=0;d<61;++d){
			if(dpB[d]>0){
				for(int k=0;k<6;++k){
					if(i<K){
						dpT[ d + meB[k] ] += dpB[d];
					}else{
						dpT[ d + meA[k] ] += dpB[d];
					}
				}
			}
		}
		rep(zz,61){
			dpB[zz] = dpT[zz];
			dpT[zz] = 0;
		}
	}
	//
	LL ans = 0;
	for(int B=1;B<61;++B){
		if(dpB[B]>0){
			for(int A=B-1;A>0;--A){
				ans += dpA[A]*dpB[B];////ココ
			}
		}
	}
	LD ansLD =1.0*ans;
	rep(i,2*N){
		ansLD /=6;//6^N
	}
	P(ansLD);
	return 0;
}
0