Стабильный расчет больших биномиальных коэффициентов в C++ ⇐ C++

[Anonymous]

Я ищу способ надежного вычисления больших биномиальных коэффициентов в C++. Моя текущая реализация может обрабатывать числа примерно до n=60, но после этого f * (n - i) превышает предел long long. Я хочу дойти до n=100, после чего результирующие биномиальные коэффициенты могут достигать 10E29.
Мне бы хотелось, чтобы точность результата была такой же высокой. насколько это возможно, но было бы нормально передать менее точный результат, если биномиальные коэффициенты очень велики (например, 15 значащих цифр). Но я не знаю, как реализовать это с каким-либо фиксированным типом возвращаемого значения — long long или другим.
long long comb(long long n, long long r)
{
long long f = 1;
for(auto i = 0; i < r; i++){
f = (f * (n - i)) / (i + 1);
}
return f ;
}


Подробнее здесь: https://stackoverflow.com/questions/787 ... ients-in-c