%I#16 2022年6月5日10:42:50
%S 7,17,26,31,49,71,97,99127161199241248287337362391449485,
%电话:5115776477217998468819671057115112491351145715671681,
%电话:179919212024204721772311244925912737288728893041319933613363
%N在Pellian方程x^2-D*y^2=1中,具有x不是最小可能值的性质的数x,其中D=(x^2-1)的平方自由部分。
%C或者将k编号为A033314(k)<>A068310(k)。
%H Eric Weisstein的数学世界,<a href=“http://mathworld.wolfram.com/PellEquation.html“>Pell方程式</a>。
%e a(2)=17。17^2-1=288的平方自由部分是D=2。但是x^2-2*y^2=1的最小可能解不是x=17,而是x=3(y=2)。
%e15不是一个项:15^2-1=224的无平方部分是D=14,x^2-14*y^2=1确实有最小解x=15(和y=4)。
%t平方自由部分[n]:=
%t倍@@Power@@@({#[[1]],Mod[#[2]],2]}&/@FactorInteger[n]);
%t a={};NMAX=3400;dict//清除;
%t对于[n=2,n<=NMAX,n++,s=平方自由部分[n^2-1];
%t如果[!IntegerQ[dict[s]],dict[s]=1,AppendTo[a,n]]];一
%K nonn公司
%O 1,1号机组
%A _ Herbert Kociemba,2022年6月2日
|