factorial - How to implement Stirling's formula with BigDecimal and BigInteger in Kotlin? -


i'm trying factorial program. if input number below 250 000 use tail recursive function find factorial of number. if input number exceed value of 250 000 try use stirling's formula (enter image description here). want able work bigintegers , bigdecimals, whenever try calculate 250 102 nan (not number ... error). can please me?

here code in kotlin:

import java.io.file import java.math.bigdecimal import java.math.biginteger  tailrec fun tail_recursion_factorial(n: biginteger, factorialofn: biginteger = biginteger.valueof(2)): biginteger {     return when(n){         biginteger.zero -> biginteger.one         biginteger.one ->  biginteger.one         biginteger.valueof(2) ->  factorialofn         else -> tail_recursion_factorial(n.minus(biginteger.one), n.times(factorialofn))     } }  // calculate approximate value of n! //  using stirling's approximation:  fun stirling_factorial(n: biginteger): bigdecimal {     return bigdecimal.valueof(math.sqrt((2*n.todouble()*1.0/3.0)*math.pi))*bigdecimal.valueof(math.pow(n.todouble(),n             .todouble()))*bigdecimal.valueof(math.pow(math.e,-1.0*n.todouble())) }   fun main(args: array<string>){     print("n == ")     var n = readline()     try {         when {             biginteger(n) < biginteger.zero -> println("sorry bro! can't factorial negative number.")             biginteger(n) >= biginteger.zero -> {                 when{                     biginteger(n) <= biginteger.valueof(250000) -> {                         file("factorial.txt").writetext("\"$n! ${tail_recursion_factorial(biginteger(n))}\"")                         println("check factorial.txt in project directory!")                     }                     else -> {                         println("since number bigger 250 000\nwe're calculating using stirling's formula wich approximation of factorial.\n")                         file("factorial.txt").writetext("the aproximation of $n! is\n${stirling_factorial                         (biginteger(n))}")                         println("check factorial.txt in project directory!")                     }                 }              }         }     }catch (e: numberformatexception){         println("sorry bro! can't ...")     } }

edit: problem solved

import ch.obermuhlner.math.big.bigdecimalmath import java.io.file import java.math.bigdecimal import java.math.mathcontext

and .... 

fun stirling_formula(n: bigdecimal) : bigdecimal {     val mathcontext = mathcontext(121)     return (bigdecimal.one + bigdecimal.one.divide(bigdecimal.valueof(12).times(n),mathcontext))             .times(bigdecimalmath.sqrt((n.times(bigdecimal.valueof(2.0))).times(bigdecimalmath.pi(mathcontext)), mathcontext))             .times(bigdecimalmath.pow(n,n,mathcontext))             .times(bigdecimalmath.pow(bigdecimalmath.e(mathcontext),-n,mathcontext))

all goes alexander romanov

ps: i've updated formula more accurate approximation one: enter image description here

you need proper library such large computations. take @ kotlin-big-math.

val = math.pow(250102.0, 250102.0) println(a) // >> infinity  val b = bigdecimal(a) // exception in thread "main" java.lang.numberformatexception:                       // infinite or nan

Comments

Post a Comment

Popular posts from this blog

javascript - WinJS appendTextAsync producing scheduler errors -

im trying implement logger in uwp app (windows 10) following code in both this , this answers. i'm afraid in both cases end getting scheduler ( winjs.utilities.scheduler ) errors @ random times (couldn't establish pattern). main problem the application freezes sometimes , suspect due log implementation because final line reads transitioning job (1) from: running to: blocked when happens. i messages like transitioning job (1) from: blocked_waiting to: cooperative_yield transitioning job (1) from: cooperative_yield to: scheduled in middle of execution, varying states. initialization method this.init = function () { let logger = new winjs.promise(function (complete) { let = new date(); let logfilename = `log-${now.getdate()}-${now.getmonth() + 1}-${now.getfullyear()}.log`; localfolder .createfolderasync("logs", windows.storage.creationcollisionoption.openifexists) .then(fun...
Read more

minify - Minimizing css files -

i want know how css files can minimized. currently, using external sites css minifier compress css files. however, after removing whitespaces, comments, etc, file sizes reduced 10-15%. how frameworks bootstrap manage make minimized versions kilobytes. as example, css is: body { background-color: #ff0; } the reduction goes 28 bytes 22 bytes other css files have lot more characters still smaller this. actually minifying code means it's remove thing code white space characters new line characters comments block delimiters it reduces amount of code has transferred on web , speed site load fastly.but unreadable format.
Read more

Sockets with kotlin -

i looking develop application in kotlin uses sockets, unfortunately kotlinx sockets broken , not far connect current kotlin 1.1 co-routines. (i see noted unanswered issue on github @sundwarf on july 16 already) specifically, sample code kotlinx.sockets uses connect in form: asocket().tcp().connect(inetsocketaddress(inetaddress.getbyname("google.com"), 80)) which fails noted trace here . i looking either sample of connect work, or recommendation of alternate package kotlix.sockets. there many kotlin apps making use of sockets, either there alternatives kotlinx.sockets or there way around connect problem.
Read more