結果
問題 | No.129 お年玉(2) |
ユーザー |
|
提出日時 | 2016-09-08 11:24:16 |
言語 | Python2 (2.7.18) |
結果 |
AC
|
実行時間 | 22 ms / 5,000 ms |
コード長 | 782 bytes |
コンパイル時間 | 331 ms |
コンパイル使用メモリ | 6,912 KB |
実行使用メモリ | 6,784 KB |
最終ジャッジ日時 | 2024-11-28 00:34:23 |
合計ジャッジ時間 | 2,333 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 46 |
ソースコード
MOD = 10**9 def ncr(n, r): if n - r < r: r = n-r if r == 0: return 1 if r == 1: return n numerator = [0 for i in range(r)] denominator = [0 for i in range(r)] for k in range(r): numerator[k] = n - r + k + 1 denominator[k] = k + 1 for p in range(2, r+1): pivot = denominator[p-1] if pivot > 1: offset = (n-r) % p for k in range(p-1, r, p): numerator[k-offset] /= pivot denominator[k] /= pivot result = 1 for k in range(r): if numerator[k] > 1: result = (result*numerator[k]) % MOD return result N,M = input(), input() N -= N/(M*1000)*M*1000 N /= 1000 if N >= M: print pow(2, M) % MOD else: print ncr(M, N)