求特殊方程的正整数解

编写程序，对任意给定的正整数N，求方程 X^2 + Y^2 = NX2+Y2=N

正整数N（≤10000）。

输出方程的全部正整数解，其中X≤Y。

每组解占1行，两数字间以1空格分隔，按X的递增顺序输出。

例如 884 输出

10 28
20 22

上代码

#include<stdio.h>
#include<math.h>
int main(int argc,char *argv[])
{   int x,y=0,n,z=0;
    scanf("%d",&n);
    for(x=1;x<=sqrt(n);x++)
    {
          y=sqrt(n-x*x);
          if(y>=x&&x*x+y*y==n)
            {
            
            printf("%d %d\n",x,y);
            z++;}
            
            
            
        
    }
    
    if(z==0) 
printf("No Solution");
}
注意要点，X≤Y。