結果

問題 No.302 サイコロで確率問題 (2)
コンテスト
ユーザー koba-e964
提出日時 2016-03-24 17:18:03
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 1,344 bytes
コンパイル時間 953 ms
コンパイル使用メモリ 86,880 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-02 00:11:16
合計ジャッジ時間 1,737 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 8 WA * 11 RE * 1
権限があれば一括ダウンロードができます

ソースコード

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

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