Produza programaticamente todas as soluções de número inteiro positivo da equação x * x + y * y = 1989

#include<stdio.h>
#include<math.h>
int main()
{
    
    
	int x,y;
	for(x=1;x<=sqrt(1989);x++){
    
    
		for(y=1;y<=sqrt(1989);y++){
    
    
			if(x*x+y*y==1989){
    
    
				printf("x=%d,y=%d\n",x,y);
			}
		}
	}
	return 0;
}