昇圧二項係数のログバージョン？

Question

かなり簡単な質問ですが、私はオンラインで何かを見つけるのに困っています。

ログコードboost::math::binomial_coefficientを探していますが、私のコードで生成された大きな値の一部を扱うことができます。 Rで

これを提供lchooseがある、しかし私は、私はあなたがそれを自分で計算する必要が考えるすべてのヘルプははるかに

Answer 1

をいただければ幸いですC++

に相当するものを見つけるように見えることはできません。しかしthis questionはそれをする方法を尋ね、それはあまりにも悪くないことが判明します。最初の近似値は次のとおりです。

n log n − m log m − (n−m) log (n−m) 

より良い近似がログのStirling approximationの別の用語を使用することによって得ることができるあなたが潜在的にある二つの用語を引いているこれらのすべてと

n log n − m log m − (n−m) log (n−m) + ½(log n − log m − log(n−m) − log 2π) 

（N！）同様のサイズであるため、精度の低下に脆弱です。