結果

問題 No.895 MESE
ユーザー aaaaaaaaaa2230
提出日時 2022-05-06 13:58:54
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 104 ms / 2,000 ms
コード長 788 bytes
コンパイル時間 140 ms
コンパイル使用メモリ 82,256 KB
実行使用メモリ 86,912 KB
最終ジャッジ日時 2024-07-05 15:48:25
合計ジャッジ時間 3,174 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

a,b,c = map(int,input().split())
n = a+b+c
mod = 10**9+7

### for bigger prime 
N = n+5
fact = [1]*N
finv = [1]*N
 
for i in range(2,N):
    fact[i] = (fact[i-1]*i)%mod
finv[-1] = pow(fact[-1],mod-2,mod)
for i in range(1,N)[::-1]:
    finv[i-1] = (finv[i]*i)%mod

def nCr(n,r):
    if r > n or r < 0:
        return 0
    else: 
        return fact[n]*finv[r]%mod*finv[n-r]%mod

ans = 0
two = pow(2,c-1,mod)
for i in range(c,n):
    cleft = i-1
    numc = nCr(cleft,c-1)
    numb = nCr(n-1-c,b)
    if i-c >= b:
        numb -= nCr(i-c,b)
        numb %= mod
    base = numc*numb
    ans += two*numc*numb%mod
    each = nCr(cleft-1,c-2)
    # print(i,numb,numc,each)
    ans += (two-1)*numb%mod*each%mod
    ans %= mod
    two = (two*2)%mod
# ans *= pow(2,mod-2,mod)
ans %= mod
print(ans)
0