Si conoces a
y estás probando valores de c
, entonces puedes calcular el valor de b
y evitar probar todos sus valores posibles:
b = (round . sqrt . fromIntegral) (c^2 - a^2)
Veamos si se puede hacer mejor...
Empezando por las matemáticas, tenemos que para números reales ℝ:
∀ a ∈ ℝ , ∃ b,c ∈ ℝ ⇒ a² + b² = c²
En cambio, para números enteros ℤ sólo se cumple para algunos casos que son las llamadas ternas pitagóricas.
Veamos algunas restricciones que podemos deducir para b
y c
:
De a² + b² = c²
, se ve que b
será siempre estrictamente menor que c
: 1 ≤ b < c
que, expresado en haskell sería: b <- [1..c-1]
.
Así mismo, podemos deducir las siguientes restriciones para c
: a < c ≤ (a²+1)/2
que en haskell sería:
c <- [a+1..a^2+1 `div` 2]
Todo junto:
ternasPitagoricas :: Integer -> [(Integer, Integer, Integer)]
ternasPitagoricas a = [(a,b,c) | c <- [a+1..(a^2+1) `div` 2]
, let b = (round . sqrt . fromIntegral) (c^2 - a^2)
, a^2 + b^2 == c^2]
Por último, cuando dices:
Necesito que la función sirva para ternasPitagoricas 20^100000, pero no llego a conseguirlo.
Antes imponerte este objetivonción que necesitas, tienes que analizar mejor qué estás pidiendo, empezando por crear una función que funcione para pequeños casos y estudiar cómo es posible escalar para alcanzar objetivos mayores.
El número 20^100000
es un número gigantesco. La función anterior tendría que recorrer todos los números enteros desde este número hasta su cuadrado, 20^200000
, lo que seguramente serían unos cuantos siglos de cálculo.
Pero hay otro problema aún más irresoluble: este número es imposible de representar como número de coma flotante sin usar librerías especializadas. Si te das cuenta, el cálculo de b
usa sqrt
, que no es usable para números muy muy grandes.
No conozco formas de calcular más rápido las ternas. Realmente es un problema matemático y aquí estamos para resolver problemas de programación.
Control-K
). Tienes ayuda en es.stackoverflow.com/editing-helpn
, dando valores ab
puedes calcularc
reduciendo la complejidad algorítmica.Haskell
lo que me dices. Es decir, no se escribirlo. Muchas gracias, y la proxima vez escribiré todo como un bloque (no sabía cómo se hacía).