結果

問題 No.302 サイコロで確率問題 (2)
コンテスト
ユーザー ciel
提出日時 2015-11-15 01:11:08
言語 C++11
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=gnu++11 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
WA  
実行時間 -
コード長 679 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 353 ms
コンパイル使用メモリ 70,400 KB
実行使用メモリ 6,400 KB
最終ジャッジ日時 2026-04-04 06:30:57
合計ジャッジ時間 3,767 ms
ジャッジサーバーID
（参考情報） 		judge2_0 / judge5_1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 19 WA * 1
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <vector>
#include <algorithm>
#include <cstdio>
#include <cmath>
using namespace std;

long double probably(int n, int s, long long l, long long r){
	vector<long double> a(s*(n+1)+1);
	for(int i=1;i<=s;i++)a[i+s]=powl(1.0/s,n);
	for(int e=0;e<n-1;e++)for(int i=s*n;i>=0;i--){
		long double sum=0;
		for(int j=i;j<i+s;j++)sum+=a[j];
		a[i+s]=sum;
	}
	long double ret=0;
	for(;l<=min(r,(long long)a.size());l++)ret+=a[s+l];
	return ret;
}

int main(){
	long long N,L,R;
	scanf("%lld%lld%lld",&N,&L,&R);
	if(N<10000){
		printf("%Lf\n",probably(N,6,L,R));
	}else{
		double a=N*3.5,v=N*35.0/12,den=sqrt(2*v);
		printf("%f\n",( erf((R+0.5-a)/den)-erf((L-0.5-a)/den) )/2);
	}
}
0