%I#15 2015年11月11日05:39:17

%S 0,0,0,1,0,0 0,0，0,0 0，0,1,0，1,0，1,0,0，

%T 1,1,1,1,0,1,2,1,0,0,2,1,1,0,2,0,1,2,0,1,0,1,1,0,1,1,0,0,1,1,1,0，

%U 1,2,1,0,1,3,1,0,2,1,0,1,1,1,1,1,2,0,0,3,1,0,1,2,0,0,1,2,0,0,1,2,0,0,0,2,0,1,4,0,2,2,2,1

%N将N划分为3个不同正方形的数量。

%H Robert Israel，n的表，a（n）表示n=0..10000</a>

%F G.F.：（1+theta_3（x））^3/48-（1+theta_3_罗伯特·伊斯雷尔（Robert Israel），2015年11月11日

%p N：=1000；#获取a（0）到a（N）

%p A:=矢量（N+1）；

%p表示从2到地板（sqrt（N））do

%p表示b从1到min（a-1，楼层（sqrt（N-a^2））do

%p表示c从0到min（b-1，楼层（平方（N-a^2-b^2））do

%px:=a^2+b^2+c^2；

%p A[x+1]：=A[x+1]+1

%p od od od：

%p转换（A，列表）；#_罗伯特·伊斯雷尔（Robert Israel），2015年11月11日

%t n=100；清除[A]；A[_]=0；对于[a=2，a<=Floor[Sqrt[n]]，a++，对于[b=1，b<=Min[a-1，Floor[Sqrt[n-a^2]]，b++，对于[c=0，c<=Min[b-1，Floor[Sqrt[n-a|2-b^2]]，c++，x=a^2+b^2+c^2；A[x+1]=A[x+1]+1]]；数组[A，n]（*_Jean-François Alcover_，2015年11月11日，改编自_Robert Israel_的枫叶脚本*）

%K nonn公司

%O 0,27号

%A·热心的W·威尔逊_