結果

問題 No.58 イカサマなサイコロ
コンテスト
ユーザー IL_msta
提出日時 2015-06-29 09:49:24
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 2 ms / 5,000 ms
コード長 1,497 bytes
コンパイル時間 1,126 ms
コンパイル使用メモリ 86,432 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-07 21:08:55
合計ジャッジ時間 1,379 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 10
権限があれば一括ダウンロードができます

ソースコード

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