結果

問題 No.58 イカサマなサイコロ
ユーザー ferinferin
提出日時 2017-02-12 21:13:16
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 2 ms / 5,000 ms
コード長 1,361 bytes
コンパイル時間 1,066 ms
コンパイル使用メモリ 145,220 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-08-28 14:54:58
合計ジャッジ時間 1,643 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

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

ソースコード

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