Por que no me funciona y como lo corrijo??

[jairogon]

El siguiente programa compila bien pero al ejecutarlo sale lo siguiente:

sqrt: Domain error
floating point errot: Domain
Abnormal program termination


Como lo corrijo

Gracias

Código:

//Programa que calcula las ternas pitagoricas
// del 1 al 500

#include<iostream.h>
#include<stdlib.h>
#include<conio.h>
#include<math.h>

int main()
{
  double p=0.0;

  system("cls");

	for(int i=1;i<=500;i++)
	{
	    for(int j=1;j<=500;j++)
	     {
		 p=(i*i)+(j*j);

		 for(int k=1;k<=sqrt(p);++k)
                 {
			if((p/(k*k)==1.0))
                        {
			 cout<<"Los numeros "<<i<<" "<<j<<" "<<k<<"  son una terna pitagorica";
                         getche();
                        }
                 }
            }

	}
	return 0;
}