Программа для решения A*x^2+B*y^2=C

[GHOST998]

Помогите написать программу на С.
Написать программу, выводящую все целочисленные решения уравнения A*x^2+B*y^2=C, лежащие в диапазоне от -1000 до 1000. Решения отсортировать по возрастанию сначала по х, а потом - по y. Вариант A=B=C=0 не рассматривается. Если решений нет выводит NO. Входной файл содержит три целых числа A,B,C. По модулю не больше миллиона.

[ser70]

Вот набросал кое-что:

#include <iostream>
using namespace std;

int main()
{
const int min=-1000;
const int max=1000;

int a=0, b=0, c=0, k=0, res;
int *masx;
int *masy;
int x=0;

cout << "Введите коэффиценты: "; cin >> a >> b >> c;


for(int i=1; i<=max; i++)
for(int j=1; j<=max; j++)
{
res=a*i*i+b*j*j;
if(res==c)
{
k++;
masx=new int[k];
masy=new int[k];

masx[k-1]=i;
masy[k-1]=j;
}
}

for(int i=min; i>=1; i++)
for(int j=min; j>=1; j++)
{
res=a*i*i+b*j*j;
if(res==c)
{
k++;
masx=new int[k];
masy=new int[k];

masx[k-1]=i;
masy[k-1]=j;
}
}

if(k==0)
{
cout << "NO" << endl;
return 0;
}

int temp;

for(int i=0; i<k; i++)
for(int j=i; j<k; j++)
if(masx[i]>masx[j])
{
temp=masx[i];
masx[i]=masx[j];
masx[j]=temp;
}

for(int i=0; i<k; i++)
cout << masx[i] << ' ';
cout << endl;

for(int i=0; i<k; i++)
for(int j=i; j<k; j++)
if(masy[i]>masy[j])
{
temp=masy[i];
masy[i]=masy[j];
masy[j]=temp;
}

for(int i=0; i<k; i++)
cout << masx[i] << ' ';
cout << endl;

return 0;
}