結果

問題 No.58 イカサマなサイコロ
コンテスト
ユーザー ferin
提出日時 2017-02-12 21:13:16
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 2 ms / 5,000 ms
コード長 1,361 bytes
コンパイル時間 2,182 ms
コンパイル使用メモリ 159,096 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-12-29 16:58:13
合計ジャッジ時間 2,214 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 10
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
typedef vector<int> VI;
typedef vector<VI> VVI;
typedef vector<ll> VL;
typedef vector<VL> VVL;
typedef pair<int, int> PII;

#define FOR(i, a, n) for (int i = (int)a; i < (int)n; ++i)
#define REP(i, n) FOR(i, 0, n)
#define ALL(x) x.begin(), x.end()
#define MOD 1000000007
#define INF 1000000000
#define PI 3.14159265359
#define EPS 1e-12

double dp1[11][65], dp2[11][65];

int main(void)
{
  int n, k;
  cin >> n >> k;

  dp1[0][0] = dp2[0][0] = 1.0;
  for(int i=1; i<=k; ++i) {
    for(int j=6*i; j>=0; --j) {
      for(int k=4; k<=6; ++k) {
        if(j >= k)
          dp1[i][j] += dp1[i-1][j-k] / 3;
      }
    }
  }
  for(int i=k+1; i<=n; ++i) {
    for(int j=6*i; j>=0; --j) {
      for(int k=1; k<=6; k++) {
        if(j >= k)
          dp1[i][j] += dp1[i-1][j-k] / 6;
      }
    }
  }
  for(int i=1; i<=n; ++i) {
    for(int j=6*i; j>=0; --j) {
      for(int k=1; k<=6; ++k) {
        if(j >= k)
          dp2[i][j] += dp2[i-1][j-k] / 6;
      }
    }
  }
  //for(int i=1; i<=12; ++i) cout << dp1[i] << " "; cout << endl;
  //for(int i=1; i<=12; ++i) cout << dp2[i] << " "; cout << endl;

  double ans = 0.0;
  for(int i=1; i<=n*6; ++i) {
    for(int j=0; j < i; ++j) {
      ans += dp1[n][i] * dp2[n][j];
    }
  }

  cout << fixed << setprecision(9) << ans << endl;

  return 0;
}
0