Triples de pitagoras C++ (Ayuda sobre optimizacion)

[erest0r]

Buenas noches, he hecho un programa sobre Triples de Pitagoras, ¿Es posible que exista una forma mas optimizada? Dejare mi codigo:

Código (cpp) [Seleccionar] Expandir


/*(Triples de Pitágoras) Un triángulo recto puede tener lados cuyas longitudes sean valores enteros. Un conjunto de tres
valores enteros para los lados de un triángulo recto se conoce como triple de Pitágoras. Estos tres lados deben satisfacer la
relación que establece que la suma de los cuadrados de dos lados es igual al cuadrado de la hipotenusa. Encuentre todos los
triples de Pitágoras para lado1, lado2, y la hipotenusa, que no sean mayores de 500. Use un ciclo for triplemente anidado
para probar todas las posibilidades. Este método es un ejemplo de la computación de fuerza bruta. En cursos de ciencias
computacionales más avanzados aprenderá que existen muchos problemas interesantes para los cuales no hay otra metodología
algorítmica conocida, más que el uso de la fuerza bruta. */

#include <iostream>
#include <cmath>
#include <iomanip>

using namespace std;

int main( int argc, char* args[] )
{
   int lado1 = 1, lado2 = 1, hipotenusa = 1;
   double res_pow_lado2, res_pow_hipotenusa;

   cout << "Lado1 al cuadrado ------ Lado2 al cuadrado ------ Hipotenusa al cuadrado" << endl << endl;

   for( ;; hipotenusa++ )
   {
       res_pow_hipotenusa = pow( hipotenusa, 2 );
       for( ;; lado2++ )
       {
           res_pow_lado2 = pow( lado2, 2 );
           for( ;; lado1++ )
           {
               if( pow( lado1, 2 ) + res_pow_lado2 == res_pow_hipotenusa )
               {
                   cout << lado1 << " = " << pow( lado1, 2 ) << setw(20) << lado2 << " = " << res_pow_lado2 << setw(20) << hipotenusa << " = " << res_pow_hipotenusa << endl;
               }
               if( pow( lado1, 2 ) + res_pow_lado2 > 500 )
                   break;
           }
           lado1 = 1;
           if( lado1 + res_pow_lado2 > 500 )
               break;
       }
       lado2 = 1;
       if( res_pow_hipotenusa > 500 )
           break;
   }

   cin.get();
return 0;
}



Gracias.

EDIT: Me gustaria saber si existe una manera de comprobar que ese es el resultado correcto :S, es decir, si no existen mas combinaciones.

EDIT2: Me di cuenta que las variables res_pow... deben ser double.

[ecfisa]

Hola erest0r.

Proba de este modo:

Código (cpp) [Seleccionar] Expandir

#include <iostream>
#include <iomanip>

#define MAX 500

using namespace std;

int main()
{
 cout << "   h" << "   c1" << "   c2" << "    V" << endl;
 cout << "---------------------" << endl;

  for (int h = 1; h*h <= MAX; h++) {
   for(int a = 1; a < h; a++)
     for(int b = 1; b < h; b++)  // (*)
       if (h*h == a*a + b*b)
         cout << "(" << setw(3) << h << "," << setw(3) << a
              << "," << setw(3) << b << " ) = " << h*h << endl;
 }
 return 0;
}

Pero hay algo del enunciado (comentado en el código) que no me quedó claro y es que si se toman como válidas aquellas ternas que tengan iguales valores en catetos diferentes, como por ejemplo:

Código [Seleccionar] Expandir


h = 5, c1 = 3, c2 = 4
h = 5, c1 = 4, c2 = 3


De no ser así, tendrías que cambiar la línea (*) por:

Código (cpp) [Seleccionar] Expandir

 ...  
   for(int b = 1; b < a; b++)
   ...


Saludos

[erest0r]

Hola ecfisa, no sabria decirte si es necesario o no mostrar ambos catetos con valores invertidos ya que ese ejercicio lo tome de un libro PDF que estoy practicando, me gusto tu solucion, se ve mucho mas rapido de entender que esa maraña de codigo que hice  , gracias.