NCk mod p dans Ruby

introduction

Lorsque je résolvais la question précédente d'AtCoder (ABC145 --D), j'ai eu du mal à trouver nCk mod p. Par conséquent, Comment trouver le coefficient binomial (nCk mod. P) et l'élément inverse (a ^ -1 mod. P) qui sont souvent effectués Je vais laisser ici la classe ModCombination facile à utiliser, qui est une réécriture de.

code

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

Comment utiliser

--Récrire la valeur de @ MOD si MOD est autre que 10 $ ^ 9 + 7 $ --nCk% MOD peut être trouvé par com (n, k)

max = 10 ** 6
comb = ModCombination.new(max)
puts comb.com(n, k) # nCk mod 10 ** 9 + 7

en conclusion

Je vais l'omettre parce que je pense qu'il sortira si je vérifie ce que fait le code ci-dessus. J'espère que cet article mène à l'AC de quelqu'un.