結果
問題 |
No.895 MESE
|
ユーザー |
![]() |
提出日時 | 2020-10-19 17:36:19 |
言語 | Crystal (1.14.0) |
結果 |
AC
|
実行時間 | 53 ms / 2,000 ms |
コード長 | 2,225 bytes |
コンパイル時間 | 11,854 ms |
コンパイル使用メモリ | 295,808 KB |
実行使用メモリ | 11,508 KB |
最終ジャッジ日時 | 2024-06-30 21:23:31 |
合計ジャッジ時間 | 13,265 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 26 |
ソースコード
lib C fun strtoll(s : UInt8*, p : UInt8**, b : Int32) : Int64 end class String def to_i64 C.strtoll(self, nil, 10) end end struct ModInt @@MOD = 1_000_000_007i64 def self.mod @@MOD end def self.zero ModInt.new(0) end def initialize(n) @n = n.to_i64 % @@MOD end getter n : Int64 def + : self self end def - : self ModInt.new(n != 0 ? @@MOD - @n : 0) end def +(m) ModInt.new(@n + m.to_i64 % @@MOD) end def -(m) ModInt.new(@n - m.to_i64 % @@MOD) end def *(m) ModInt.new(@n * m.to_i64 % @@MOD) end def /(m) raise DivisionByZeroError.new if m == 0 a, b, u, v = m.to_i64, @@MOD, 1i64, 0i64 while b != 0 t = a // b a -= t * b a, b = b, a u -= t * v u, v = v, u end ModInt.new(@n * u) end def //(m) self / m end def **(m : Int) t, res = self, ModInt.new(1) while m > 0 res *= t if m.odd? t *= t m >>= 1 end res end def ==(m) @n == m.to_i64 end def !=(m) @n != m.to_i64 end def succ self + 1 end def pred self - 1 end def to_i64 : Int64 @n end delegate to_s, to: @n delegate inspect, to: @n end struct Int def to_mint : ModInt ModInt.new(self) end end class Comb(T) def initialize(@Max = 300009) @fac = Array(T).new(@Max, T.zero) @finv = Array(T).new(@Max, T.zero) @inv = Array(T).new(@Max, T.zero) @fac[0] = @fac[1] = T.new(1) @finv[0] = @finv[1] = T.new(1) @inv[1] = T.new(1) (2...@Max).each do |i| @fac[i] = @fac[i - 1] * i @inv[i] = -@inv[T.mod % i] * (T.mod / i) @finv[i] = @finv[i - 1] * @inv[i] end end def permutation(n, r) (n < r || n < 0 || r < 0) ? 0 : @fac[n] * @finv[n - r] end def combination(n, r) (n < r || n < 0 || r < 0) ? 0 : @fac[n] * @finv[r] * @finv[n - r] end def factorial(n) @fac[n] end end a, b, c = read_line.split.map(&.to_i) n = a + b + c combi = Comb(ModInt).new puts (2..a + 1).sum { |i| left, right = i - 1, n - i x, y, z = a - left, b - 1, c m = x + y + z z_cnt = combi.combination(m, z) * z / m (2.to_mint**m - 1) * z_cnt * combi.combination(x + y, x) }