连分数与佩尔方程特解（最小整数解）

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll a[20000];
bool pell_minimum_solution(ll n,ll &x0,ll &y0){//求PELL方程最小整数解
	ll m=(ll)sqrt((double)n);//M是N的平方根向下取整
	if(m*m==n)return false;//当n是完全平方数则佩尔方程无解（不讨论正负一，零）
	//下面是把N用连分数形式存，B,C,TMP（即AI）的递推见解释以22为例参考
	int i=0;    //连分数的数位
	a[i++]=m;   //A0位整数部分m=4
	ll b=m,c=1; //B=4即整数部分位,C=1即求RN时的分母
	double sq=sqrt(n);//SQ是N的高精度根,相当于r0
	double tmp;//tp在下面的循环就是r