結果

問題 No.302 サイコロで確率問題 (2)
ユーザー koba-e964koba-e964
提出日時 2016-03-24 17:26:41
言語 C++11
(gcc 11.4.0)
結果
RE  
実行時間 -
コード長 1,339 bytes
コンパイル時間 651 ms
コンパイル使用メモリ 88,688 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-10-02 00:11:30
合計ジャッジ時間 4,225 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 58 ms
5,248 KB
testcase_01 AC 1,128 ms
5,248 KB
testcase_02 AC 2 ms
5,248 KB
testcase_03 AC 102 ms
5,248 KB
testcase_04 AC 1,123 ms
5,248 KB
testcase_05 AC 1 ms
5,248 KB
testcase_06 AC 2 ms
5,248 KB
testcase_07 AC 2 ms
5,248 KB
testcase_08 AC 114 ms
5,248 KB
testcase_09 AC 2 ms
5,248 KB
testcase_10 AC 2 ms
5,248 KB
testcase_11 AC 152 ms
5,248 KB
testcase_12 AC 2 ms
5,248 KB
testcase_13 AC 2 ms
5,248 KB
testcase_14 AC 2 ms
5,248 KB
testcase_15 AC 2 ms
5,248 KB
testcase_16 RE -
testcase_17 AC 2 ms
6,816 KB
testcase_18 AC 2 ms
6,820 KB
testcase_19 AC 2 ms
6,816 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <algorithm>
#include <bitset>
#include <cassert>
#include <cctype>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <ctime>
#include <deque>
#include <functional>
#include <iomanip>
#include <iostream>
#include <list>
#include <map>
#include <numeric>
#include <queue>
#include <set>
#include <sstream>
#include <stack>
#include <string>
#include <utility>
#include <vector>

#define REP(i,s,n) for(int i=(int)(s);i<(int)(n);i++)

using namespace std;
typedef long long int ll;
typedef vector<int> VI;
typedef pair<int, int> PI;
const double EPS=1e-9;
const int LIM = 8000;

double dp[2][5 * LIM + 1];


double calc(int n, int l, int r) {
  dp[0][0] = 1.0;
  REP(i, 0, n) {
    int t = i % 2;
    REP(j, 0, 5 * i + 6) {
      dp[1 - t][j] = 0;
    }
    REP(j, 0, 5 * i + 6) {
      REP(k, 0, min(5, j) + 1) {
	dp[1 - t][j] += dp[t][j - k] / 6.0;
      }
    }
  }
  double sum = 0;
  REP(i, l, r + 1) {
    sum += dp[n % 2][i];
  }
  return sum;
}

int main(void){
  ll n, l, r;
  cin >> n >> l >> r;
  if (n < LIM) {
    printf("%.15f\n", calc(n, max(l - n, 0LL) , min(r - n, 5 * n)));
  } else {
    double lb = (l - 0.5 - 3.5 * n) / sqrt(35.0 / 6.0 * n);
    double rb = (r + 0.5 - 3.5 * n) / sqrt(35.0 / 6.0 * n);
    double res = 0.5 * (erf(rb) - erf(lb));
    printf("%.15f\n", res);
  }
}
0