結果
問題 |
No.1126 SUM
|
ユーザー |
![]() |
提出日時 | 2020-07-25 13:23:54 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 158 ms / 1,000 ms |
コード長 | 882 bytes |
コンパイル時間 | 141 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 22,716 KB |
最終ジャッジ日時 | 2024-06-27 04:01:07 |
合計ジャッジ時間 | 5,373 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 26 |
ソースコード
import sys #import collections as col input=sys.stdin.readline def linput(ty=int, cvt=list): return cvt(map(ty,input().rstrip().split())) def vinput(rep=1, ty=int, cvt=list): return cvt(ty(input().rstrip()) for _ in "*"*rep) def facto(x, y, mod): re = 1 for i in range(y, x+1): re *= i re %= mod #print(re, file=sys.stderr) return re def cmb(n, r, mod): if ( r<0 or r>n ): return 0 r = min(r, n-r) #return facto(n, n-r+1, mod) * g2[r] % mod return g1[n] * g2[r] * g2[n-r] % mod mod = 10**9+7 NN = 10**5+10 g1 = [1, 1] # factorial g2 = [1, 1] # factorial^(-1) inverse = [0, 1] # tmp for i in range( 2, NN + 1 ): g1.append( ( g1[-1] * i ) % mod ) inverse.append( ( -inverse[mod % i] * (mod//i) ) % mod ) g2.append( (g2[-1] * inverse[-1]) % mod ) def main(): A,B = linput() M = mod res = cmb(B+1, A+1, M) print(res) main()