結果

問題 No.302 サイコロで確率問題 (2)
ユーザー koba-e964koba-e964
提出日時 2016-03-24 17:17:25
言語 C++11
(gcc 11.4.0)
結果
TLE  
実行時間 -
コード長 1,345 bytes
コンパイル時間 715 ms
コンパイル使用メモリ 86,856 KB
実行使用メモリ 10,112 KB
最終ジャッジ日時 2024-10-02 00:11:14
合計ジャッジ時間 16,343 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

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

ソースコード

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 = 20000;

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 * LIM + 1) {
      dp[1 - t][j] = 0;
    }
    REP(j, 0, 5 * LIM + 1) {
      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