5-21 求特殊方程的正整数解 (15分)

#include <iostream>
#include <sstream>
#include <string.h>
#include <map>
#include <stdio.h>
#include <algorithm>
#include <math.h>
using namespace std;

int b[10050]= {0};
int main()
{
memset(b,0,sizeof(b));
int a[120]= {0},N;
scanf("%d",&N);
for(int i=1; i<=sqrt(N); i++)
{
a[i]=i*i;
b[i*i]=1;
}
int c=0;
for(int i=1; i<=sqrt(N); i++)
{
if(b[N-i*i])
{
int m=N-i*i;
m=sqrt(m);
printf("%d %d\n",i,m);
b[N-i*i]=0;
b[i*i]=0;
c=1;
}
}
if(c==0) printf("No Solution\n");
return 0;
}