%I#16 2023年2月24日11:53:08

%S 10,26,50,58,65,74,82,85106122125130145170185202218226250，

%电话：265274290298314325338346362365370394425442458481，

%电话：48549353053338554565586606266296363468568746

%N复合数字N，这样x^2-N*y^2表示-1。

%C对于k=0..，小于或等于10^k的项数：0, 1, 8, 71, 712, 6702, 63485, 602870, ... . - _Robert G.Wilson v_，2008年7月20日

%H Robert G.Wilson v，n表，n=1..63485的a（n）</a>

%H J.P.Robertson和K.R.Matthews，<a href=“http://www.jstor.org/stable/27642477“>对Feit结果的连续分数法</A>，《美国数学月刊》，115（2008年第4期），346-349。

%H Eric Weisstein的数学世界，<a href=“http://mathworld.wolfram.com/PellEquation.html“>佩尔方程式</a>。

%e3^2-10*1^2=-1，所以10是一个成员。

%e 4005^2-106*389^2=-1，所以106是一个成员。

%t lst={}；做[If[！PrimeQ@n&&FindInstance[x^2-n*y^2==-1，{x，y}，Integers]！={}，附加到[lst，n]]，{n，2，1000}]

%Y有关具有此属性的素数，请参见A002313。A134406是一个子集。

%K nonn公司

%O 1,1号机组

%A _N.J.A.Sloane，2008年4月8日

%E更多条款摘自2008年7月20日_Robert G.Wilson v_