%I#36 2019年9月13日03:58:58

%S 5,3,2,17,13,7,37,41,29,11,73,61,53,19,23109,89101,97113,31149，

%电话157137,43,47193197181173,5927724128126925723229,67,71，

%电话：7931333734935317293,83409421433389401371397103107

%N素数按二次无理连分式项排序。

%C对于每个素数p，如果p与1 mod 4全等，则计算（1+sqrt（p））/2，否则计算sqrt（p）。将其表示为周期性的连续分数。按周期部分中最大的项进行排序；在那些具有相同最大项的项中，按项的几何平均值对其进行排序。

%这些二次无理数用于Richtmyer低偏差序列生成器。按这种方式排序，黄金比率在二次无理数列表中居首位，因为（frac（n*phi））在形式序列（frac，n*a））之间的差异最小。

%H Pierre Abbat，n的表格，n=1..4228的a（n）</a>

%H Pierre Abbat，<a href=“https://github.com/phma/quadlods网站“>四边形</a>

%e 17==1（mod 4），因此计算（sqrt（17）+1）/2=2.561552812808830……其连分式展开为[2；（1,1,3）]。最大项为3。

%e 13==1（mod 4），因此计算（sqrt（13）+1）/2=2.30277563773199……其连分式展开为[2；（3）]。最大项再次为3，但平均项大于（sqrt（17）+1）/2中的平均项，因此13在17之后。

%e 7==3（mod 4），因此计算sqrt（7）=2.645751311064590……其连分式展开为[2；（1,1,4）]。最大项是4，所以7在13之后。

%o（C++）请参阅Quadlods链接。程序生成序列的6542个项，但在4228个项之后，有大于65536的项，而它没有生成这些项。

%A000040的Y排列。参见A001622（φ）。

%K nonn公司

%O 1，1

%A_Pierre Abbat，2019年9月9日