求解佩尔方程的基本解

 佩尔方程是关于x,y的二元不定式方程，表述为：x*x-n*y*y=1;(n为非平方整数)

当x=1或x=-1,y=0显然满足方程。常把x,y为满足方程的最小正数的解又称基本解。

#include<iostream>
using namespace std;
#include<math.h>
void main(){
double a,m,n;
double x,y;
cout<<"解佩尔方程：x*x-n*y*y=1"<<endl;
cout<<"请输入非平方整数n"<<endl;
cin>>n;
m=floor(sqrt(n+1));
if(m*m==n){cout<<"方程无解"<<endl;  return;}
y=1;
while(y<=10000000){
y++;
a=n*y*y;
x=floor(sqrt(a+1));
if(x*x==a+1){
cout<<"方程：x*x-n*y*y=1的基本解为"<<endl;
cout<<"x="<<(int)x<<endl;
cout<<"y="<<y<<endl;
break;
}
}
if(y>10000000) cout<<"未求出方程解"<<endl;
}