%I#31 2023年8月2日07:11:56

%第30、39、46、50、51、54、57、62、63、65、66、70、71、74、75、78、79、81、84、85、86、87、90页，

%电话：91,93,94,95,98,99102105106107109110111113114116117118，

%电话：119120121122123125126127129130131133134135137

%N是4个不同非零平方和的数字：形式为w^2+x^2+y^2+z^2，其中0<w<x<y<z。

%H T.D.Noe，n的表格，n=1..1000的a（n）</a>

%H<a href=“/index/Su#ssq”>与平方和相关的序列的索引项</a>

%F{n:A025443（n）>=1}。A025386和A025376的联合_R.J.Mathar，2018年6月15日

%e 30=1^2+2^2+3^2+4^2。

%t data=Flatten[DeleteCase[FindInstance[w^2+x^2+y^2+z^2==#&0<w<x<y<z<#，{w，x，y，z}，Integers]&/@Range[137]，{}]，1]；w ^2+x ^2+y ^2+z ^2/。数据（*_Ant King，2010年10月17日*）

%t选择[Union[Total[#^2]和/@Subsets[Range[10]，{4}]]，#<=137&]（*哈维·P·戴尔，2011年7月3日*）

%o（哈斯克尔）

%o a004433 n=a004433_列表！！（n-1）

%o a004433_list=过滤器（p4$tail a000290_list）[1..]其中

%o p k（q:qs）m=k==0&&m==0||

%o q<=m&&k>=0&&（p（k-1）qs（m-q）| | p k qs m）

%o——_ Inhard Zumkeller_，2013年4月22日

%o（PARI）列表（lim）=我的（v=列表（））；lim=1；对于（z=4，平方（lim\4），对于（y=3，min（平方（（lim-z^2）\3），z-1），对于；集（v）\\_Charles R Greathouse IV_，2017年2月7日

%Y参见A001944、A001995、A003995、P004431、A004432、A00443、A224981、A224992、A224883、A000290。

%不，简单，好

%O 1,1号机组

%A _N.J.A.斯隆_