🚀 go-pugleaf

Saludos, Jes�s Carrete Monta�a.

El 19 Jun 2001 00:42:31 +0200 nos contabas:


>        Me gustar�a si es posible calcular en tiempo razonable (menos
>de una semana en un ordenador dom�stico) la expresi�n arriba
>indicada. He probado a hacerlo en bc y apcalc, pero ambos dicen que el
>exponente es demasiado grande (a pesar de ser lenguajes de precisi�n
>arbitraria). 

Se hace en mucho menos de una semana planteando unos logaritmos decimales:

 log(2^(2^64)) = log(2)*2^64 = .30102999566398119521 * 18446744073709551616 =
 5553023288523357132.21137434038486695936 (valores obtenidos con bc)
 
Si esto lo separamos en sus partes entera y fraccionaria, podemos hallar el
valor deseado:

 alog(5553023288523357132.21137434038486695936) =
 alog(0.21137434038486695936) * alog(5553023288523357132) = 
 1.62695050249084300164 * 10^5553023288523357132 

Como ves, es un n�mero tan exageradamente enorme que calculandolo
directamente no hay forma de manejarlo.

Ahora viene el truco: �como se hace esto con bc, si s�lo tiene el logaritmo
neperiano (l) y la exponencial con base e (e)? 

F�cil:

 log(x) = l(x)/l(10)
 alog(x) = e(x*l(10))
 
>A�n no he probado ninguna librer�a como gmp, pero me
>gustar�a antes conocer las posibilidades.

No future.

>Un saludo.

-- 
Outside of a dog, a book is a man's best friend.  Inside a dog it's too
dark to read.
                -- Groucho Marx