%I#48 2019年5月23日08:38:21
%S 2,50325110581255525105625276257182513812581250160225,
%电话122102524420501795625801125462656252082925412597656254005625,
%电话:44890625305256256105125051250512502832530346062553955950781253540972510140625836775781250
%N最小正整数,即正整数的两个平方之和,精确到N次。
%D A.Beiler,《数字理论中的娱乐》,多佛,第140-141页。
%H T.D.Noe和Ray Chandler,<a href=“/A016032/b016032.txt”>n的表,a(n)对于n=1..2178</a>(a(2179)超过1000位)。
%H C.里维拉,<a href=“http://www.primepuzzles.net/puzzles/puzz_062.htm“>拼图62</a>
%H Eric Weisstein的数学世界,<a href=“http://mathworld.wolfram.com/SquareNumber.html“>平方数</a>
%H G.Xiao,<a href=“http://wims.unice.fr/~wims/en_tool~number~twosquares.en.html“>两个方块</a>
%H<a href=“/index/Su#ssq”>与平方和相关的序列的索引项</a>
%F a(n)=最小值(2*A018782(2n-1),A018782。
%e a(0)=1,因为1是不能表示为两个正平方和的最小正整数。
%e a(1)=2从2=1^2+1^2。
%e a(2)=50从50=1^2+7^2=5^2+5^2开始。
%t数组[Block[{k=1},While[Length@DeleteCase[PowersRepresentations[k,2,2],_?(!FreeQ[#,0]&)]!=#,k++];k] &,6](*米歇尔·德弗里格,2019年3月31日*)
%o(PARI)b(k)=我的(c=0);对于(i=1,平方(k\2),如果(平方(k-i^2),c+=1));抄送\\A025426
%o表示(n=1,10,k=1;而(k,如果(b(k)==n,则打印1(k,“,”);断裂);k+=1))\\_Derek Orr_,2019年3月20日
%Y参见A018825、A048610、A025284-A025293(第一个条目)。
%Y其他版本见A000446、A124980和A093195。
%K nonn很好
%O 1,1号机组
%A _罗伯特·G·威尔逊v_
%E由_Jud McCranie更正和扩展_
%2007年11月12日,几位记者改进了E定义
|