結果
問題 | No.574 正多面体サイコロ |
ユーザー |
![]() |
提出日時 | 2020-05-24 14:37:51 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 402 ms / 2,000 ms |
コード長 | 1,319 bytes |
コンパイル時間 | 85 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 11,008 KB |
最終ジャッジ日時 | 2024-10-11 07:15:34 |
合計ジャッジ時間 | 2,742 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 20 |
ソースコード
from operator import itemgetterfrom itertools import *from bisect import *from collections import *from heapq import *import mathimport syssys.setrecursionlimit(10 ** 6)def II(): return int(sys.stdin.readline())def MI(): return map(int, sys.stdin.readline().split())def LI(): return list(map(int, sys.stdin.readline().split()))def SI(): return sys.stdin.readline()[:-1]def LLI(rows_number): return [LI() for _ in range(rows_number)]def LLI1(rows_number): return [LI1() for _ in range(rows_number)]int1 = lambda x: int(x) - 1def MI1(): return map(int1, sys.stdin.readline().split())def LI1(): return list(map(int1, sys.stdin.readline().split()))p2D = lambda x: print(*x, sep="\n")dij = [(1, 0), (0, 1), (-1, 0), (0, -1)]def nCr(n,r):if n<r:return 0return math.factorial(n)//math.factorial(r)//math.factorial(n-r)def main():f,n,k=MI()ans=0for d in range(1,f+1):p=0# 自分(d)より大きい数の回数xfor x in range(k):px=(f-d)/f# 自分(d)より小さい数の回数yfor y in range(n-k+1):py=(d-1)/f# 大きい数がx回、小さい数がy回でる確率p+=pow(px,x)*nCr(n,x)*pow(py,y)*nCr(n-x,y)*pow(1/f,n-x-y)ans+=p*dprint(ans)main()