#include #include #include #include #include #include #include #include #include // require sort next_permutation count __gcd reverse etc. #include // require abs exit atof atoi #include // require scanf printf #include #include // require accumulate #include // require fabs #include #include #include #include // require setw #include // require stringstream #include // require memset #include // require tolower, toupper #include // require freopen #include // require srand #define rep(i,n) for(int i=0;i<(n);i++) #define ALL(A) A.begin(), A.end() using namespace std; typedef long long ll; typedef pair P; /* No.58 イカサマなサイコロ ★愚直なDP解 dpT[i][j][k] := 普通のサイコロ i 個の使い イカサマサイコロ j 個使って 合計 k になる確率 dpT[0][0][0] = 1 dpT[i+1][j][k+l] += dpT[i][j][k]/6; (l = 1,2,3,4,5,6 ) dpT[i][j+1][k+l'] += dpT[i][j][k]/3; (l = 4,5,6 ) ★太郎が勝つ確率を求める F := 太郎が勝つ事象 Ei := 二郎のサイコロの合計が i の事象 P(F) = ΣP(F ∩ Ei) (全確率の公式) P(F) = P(F ∩ E1) + P(F ∩ E2) + P(F ∩ E3) + ... + P(F ∩ E60) P(F) = P(F|E1)P(E1) + P(F|E2)P(E2) + P(F|E3)P(E3) + ... + P(F|E60)P(E60) (乗法公式、条件付確率の公式) P(F) = P(F)P(E1) + P(F)P(E2) + P(F)P(E3) + ... + P(F)P(E60) (F と Ei は独立事象) ^ ^ ここは個別の確率 P(E1) の場合の P(F) だから 左辺の P(F) とは異なる! 他の P(F) も同様。 ここは求めたい確率 */ const int MAX_N = 12; const int MAX_M = 12; const int MAX_K = 70; double dpT[MAX_N][MAX_M][MAX_K]; // 太郎のDPテーブル double dpJ[MAX_N][MAX_M][MAX_K]; // 二郎のDPテーブル const double normal[7] = { 0, 1./6., 1./6., 1./6., 1./6., 1./6., 1./6. }; const double fake[7] = { 0, 0., 0., 0., 1./3., 1./3., 1./3. }; int main() { memset (dpT, 0., sizeof (dpT ) ); memset (dpJ, 0., sizeof (dpJ ) ); ios_base::sync_with_stdio(0); int N, K; cin >> N >> K; dpT[0][0][0] = 1.; for (int i = 0; i < N; i++ ){ for (int j = K; j >= 0; j-- ){ for (int k = 6*N; k >= 0; k-- ){ if (dpT[i][j][k] != 0. ){ for (int l = 1; l <= 6; l++ ){ dpT[i+1][j][k+l] += dpT[i][j][k]*normal[l]; dpT[i][j+1][k+l] += dpT[i][j][k]*fake[l]; } // end for } // end if } // end for } // end for } // end for dpJ[0][0][0] = 1.; for (int i = 0; i < N; i++ ){ for (int j = K; j >= 0; j-- ){ for (int k = 6*N; k >= 0; k-- ){ if (dpJ[i][j][k] != 0. ){ for (int l = 1; l <= 6; l++ ){ dpJ[i+1][j][k+l] += dpJ[i][j][k]*normal[l]; } // end for } // end if } // end for } // end for } // end for double res = 0.; for (int i = 1; i <= 6*N; i++ ){ for (int j = 1; j < i; j++ ){ res += dpT[N-K][K][i]*dpJ[N][0][j]; } // end for } // end for printf ("%.6lf\n", res ); return 0; }