Als ich AtCoders frühere Frage löste (ABC145 - D), fiel es mir schwer, nCk mod p zu finden. Daher So finden Sie den Binomialkoeffizienten (nCk mod. P) und das inverse Element (a ^ -1 mod. P), die häufig ausgeführt werden Ich werde die benutzerfreundliche ModCombination-Klasse hier lassen, die neu geschrieben wurde.
class ModCombination
def initialize(max)
@MOD = 10 ** 9 + 7
@fact = [1, 1]
@factinv = [1, 1]
@inv = [0, 1]
2.upto(max) do |i|
@fact << @fact[i-1] * i % @MOD
@inv << @MOD - @inv[@MOD % i] * (@MOD / i) % @MOD
@factinv << @factinv[i-1] * @inv[i] % @MOD
end
end
def com(n, k)
return 0 if n < k
return 0 if n < 0 || k < 0
return @fact[n] * (@factinv[k] * @factinv[n-k] % @MOD) % @MOD
end
end
nCk% MOD
kann durch com (n, k)
gefunden werdenmax = 10 ** 6
comb = ModCombination.new(max)
puts comb.com(n, k) # nCk mod 10 ** 9 + 7
Ich werde es weglassen, weil ich denke, dass es herauskommt, wenn ich überprüfe, was der obige Code tut. Ich hoffe, dieser Artikel führt zu jemandes AC.
Recommended Posts