結果
問題 | No.895 MESE |
ユーザー | convexineq |
提出日時 | 2021-04-08 10:04:16 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 100 ms / 2,000 ms |
コード長 | 888 bytes |
コンパイル時間 | 349 ms |
コンパイル使用メモリ | 82,224 KB |
実行使用メモリ | 86,288 KB |
最終ジャッジ日時 | 2024-06-23 14:22:25 |
合計ジャッジ時間 | 3,509 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 26 |
ソースコード
SIZE=3*10**5+1; MOD=10**9+7 #998244353 #ここを変更する inv = [0]*SIZE # inv[j] = j^{-1} mod MOD fac = [0]*SIZE # fac[j] = j! mod MOD finv = [0]*SIZE # finv[j] = (j!)^{-1} mod MOD fac[0] = fac[1] = 1 finv[0] = finv[1] = 1 for i in range(2,SIZE): fac[i] = fac[i-1]*i%MOD finv[-1] = pow(fac[-1],MOD-2,MOD) for i in range(SIZE-1,0,-1): finv[i-1] = finv[i]*i%MOD inv[i] = finv[i]*fac[i-1]%MOD def choose(n,r): # nCk mod MOD の計算 if 0 <= r <= n: return (fac[n]*finv[r]%MOD)*finv[n-r]%MOD else: return 0 a,b,c = map(int,input().split()) n = a+b+c r = 0 p = 1 for i in range(1,a+1): r += choose(n-i-1-1,c-1)*choose(a+b-i-1,b-1)%MOD ans = 0 for i in range(n): ans += r*p%MOD if n-1-i-1 <= a: j = n-1-i-1 r -= choose(n-j-1-1,c-1)*choose(a+b-j-1,b-1) r %= MOD if r==0: break p = p*2%MOD print(ans%MOD)