%I#29 2024年4月24日13:21:51

%S 0,1,2,3,18,4,5,70,6,32,7182,9929718,81068,43,93785005604,10，

%电话：40058890182776,11682,571744,12135824832118,131181111225770，

%电话：681764132、143141251、15171023156710110684443、166072、821407

%N对于A031396中给出的k值，满足Pell方程x^2-k*y^2=-1的最小整数x。

%H Ray Chandler，n的表，n的a（n）=1..10000</a>

%H K.Lakshmi和R.Someshwari，<a href=“https://www.ijeter.everscience.org/Manuscripts/Volume-4/Issue-7/Vol-4-Issue-7-M-02.pdf“>关于负佩尔方程y^2=72x^2-23，国际工程研究新兴技术期刊（IJETER），第4卷，第7期，2016年7月。

%H R.Suganya和D.Maheswari，<a href=“http://dx.doi.org/10.22457/jmi.v11a9“>关于负Pellian方程y^2=110*x^2-29，《数学与信息学杂志》，第11卷（2017年），第63-71页。

%H S.Vidhyalakshmi、V.Krithika和K.Agalya，<a href=“http://www.ijeter.everscience.org/Manuscripts/Volume-4/Issue-2/Vol-4-Issue-2-M-04.pdf“>关于负佩尔方程y^2=72x^2-8，国际工程研究新兴技术期刊（IJETER）第4卷，第2期，2016年2月。

%e a（5）=18，因为A031396（5）=13，而x^2-13y^2=-1与最小可能x的解具有x=18。

%p A130226：=进程（m）

%p局部xm，x，i，xmo，y2；

%p xm:=[]；#x^2-m*y^2=-1（mod m）需要xm[]中的x

%x从0到m-1的p do

%p如果modp（x^2，m）=modp（-1，m），则

%p xm：=[操作（xm），x]；

%p end if；

%p端do：

%p代表i从0到do

%xm-do中xmo的p

%px：=i*m+xmo；

%p y2:=（x^2+1）/m；

%p如果issqr（y2），则

%p返回x；

%p end if；

%p端do：

%p端do：

%p端程序：

%p L:=文件列表（“b031396.txt”）；

%p n：=1:

%L do中m的p

%p打印f（“%d%d\n”，n，A130226（m））；

%pn:=n+1；

%截止日期：#R.J.Mathar_，2014年10月19日

%t项=1000；

%t a031396=案例[导入[“https://oeis.org/A031396/b031396.txt网站“，”表格“]，{_，_}][[；；术语，2]]；

%t sol[n_]：=求解[x>0&y>0&x^2-ny^2==-1，{x，y}，整数]；

%ta[1]=0；a[n]：=a[n]=x/。溶胶[a031396[[n]]/。C[1]->0//第一//简化//安静；

%t表[打印[n，“”，a031396[[n]]，“”；a[n]，{n，1，术语}]（*Jean-François Alcover_，2020年4月5日*）

%Y参考A094048。

%K nonn公司

%氧1,3

%A Colin Barker_，2007年8月5日