%I#23 2024年2月26日19:17:09

%S 0,1,1,2,2,2,2,2,2,3,3,33,3,1,3,4,4,4，4,4,1,4,5,5,5，5,5,1,5,4,5，

%T第5、5、6、6、6,6、6和6、5、6,6,6和6,6，6、7、7、7,7、6、7，

%U 8,8,8,10,10,10,9,9,10,10,9，9,9,9

%N将N表示为四个平方之和，x^2+y^2+z^2+w^2，其中x>=y>=z>=w>=0，最大化x的值，那么a（N）就是x。

%拉格朗日定理告诉我们，每个正整数都可以写成四个平方和。

%H David Consiglio，Jr.，<a href=“/A178786/b178786.txt”>n表，n=0..10000的a（n）</a>

%H David Consiglio，Jr.，<a href=“/A178786/A178786.txt”>Python程序</a>

%o（Python）

%o从数学导入*

%o对于范围（0.500）内的nbre：#或大于500！

%o最大值4=0

%o对于范围（0，int（sqrt（nbre/4））+1）内的c1：

%o对于范围（c1，int（sqrt（nbre/3））+1）内的c2：

%o对于范围内的c3（c2，int（sqrt（nbre/2））+1）：

%o s3=c3**2+c2**2+c 1**2

%o如果s3<=nbre：

%o c4=平方英尺（nbre-s3）

%o如果int（c4）==c4且c4>=c3：

%o如果c4>最大c4：

%o最大c4=整数（c4）

%o打印（maxc4，end='，'）

%Y参见A122922、A122923、A122824、A122225、A122966、A122627、A002330、A129921。

%三个正方形的Y模拟：A261904和A261915。

%K nonn公司

%0、5

%A·Sébastien Dumortier，2011年6月24日