В свое время сделал очень изящное решение для подсчета длинных хвостов распределения, там как всегда факториал на экспорненту. Все считали в лоб и выходили за точность, а я сделал как раз итеративное вычисление с кэшированием и там были конечные числа