結果
| 問題 | No.1408 Nice Dice Game | 
| コンテスト | |
| ユーザー |  QCFium | 
| 提出日時 | 2021-02-26 22:12:05 | 
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) | 
| 結果 | 
                                AC
                                 
                             | 
| 実行時間 | 12 ms / 2,000 ms | 
| コード長 | 978 bytes | 
| コンパイル時間 | 1,571 ms | 
| コンパイル使用メモリ | 165,768 KB | 
| 実行使用メモリ | 5,248 KB | 
| 最終ジャッジ日時 | 2024-10-02 15:00:28 | 
| 合計ジャッジ時間 | 2,659 ms | 
| ジャッジサーバーID (参考情報) | judge5 / judge3 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 1 | 
| other | AC * 36 | 
ソースコード
#include <bits/stdc++.h>
int ri() {
	int n;
	scanf("%d", &n);
	return n;
}
int main() {
	int n = ri();
	if (n == 1) {
		puts("1.20205690315959428539973816151144");
	} else if (n == 2) {
		puts("1.0823232337111381915160036965411679027747");
	} else if (n == 3) {
		puts("1.036927755");
	} else if (n > 50) puts("1");
	else {
		long double res = 0;
		for (int b = 2; b <= 1000; b++) {
			long double x = 1;
			for (int i = 0; i < n; i++) x *= (long double) 1 / b;
			res += (1 - x) / ((long double) b * (b - 1));
		}
		res += 0.001;
		
		long double r0 = 0;
		for (int a = 2; a <= 400; a++) {
			long double cur = 0;
			for (int b = a + 1; b <= 1000; b++) {
				long double x = 1;
				for (int i = 0; i < n; i++) x *= (long double) a / b;
				cur += (1 - x) / ((long double) b * (b - a));
			}
			for (int i = 1; i <= a; i++) cur += 1.0L / (1000 - a + i) / a;
			for (int i = 1; i <= n; i++) cur /= a;
			r0 += cur;
		}
		
		printf("%.11Lf\n", r0 + res);
		
	}
	
	return 0;
}
            
            
            
        