%I#18 2021年5月30日11:11:12

%S 0,0,00,0,1,0,0，0,0'，0,0，

%温度2,2,0,3,2,3,1,2,2,3,3,4,4,0,4,3,3,1,3,6,3,0,6,4,2,4,4,1,4,6,6,6，

%U 3,3,6,4,3,7,5,5,6,1,6,2,10,4,5,7,5,7,5

%N将8*N+4表示为四个不同奇数平方和的方法的数量。

%每个奇数平方都是8*k+1形式的数，所以四个奇数平方的和都是8*k+4形式的数。

%C A316489列出了8*k+4形式的所有正数，这些正数不能用四个不同的奇数平方和表示；对于每个这样的数字，a（k）=0。

%C A316834列出了所有只能用一种方式表示为四个不同奇数平方和的数字；每个这样的数字的形式为8*k+4，对于每个这样的数，a（k）=1。

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

%e n=1:8*1+4=12不能表示为四个不同的奇数平方和，因此a（1）=0。

%en=10:8*10+4=84可以用一种方式表示为四个不同的奇数平方和（84=1^2+3^2+5^2+7^2），因此a（10）=1。

%e n=19:8*19+4=156可以用两种方式表示为四个不同的奇数平方之和（156=1^2+3^2+5^2+11^2=1^2+5 ^2+7 ^2+9 ^2），因此a（19）=2。

%pb:=proc（n，i，t）选项记忆`如果`（n=0，`如果`（t=0，1，0），

%p`if`（最小（i，t）<1，0，b（n，i-2，t）+

%p`if`（i^2>n，0，b（n-i^2，i-2，t-1））

%p端：

%pa:=n->（m->b（m，（r->r+1-irem（r，2））（isqrt（m）），4））（8*n+4）：

%p序列（a（n），n=0..100）；#_Alois P.Heinz，2018年8月5日

%t a[n_]：=计数[整数分区[8 n+4，{4}，范围[1，Sqrt[8 n%4]，2]^2]，w_/；最大@差异@w＜0]；阵列[a，87，0]（*Giovanni Resta_，2018年8月12日*）

%tb[n_，i_，t_]：=b[n，i，t]=如果[n==0，如果[t==0、1、0]，

%t如果[Min[i，t]<1，0，b[n，i-2，t]+

%t如果[i^2>n，0，b[n-i^2，i-2，t-1]]]；

%t a[n]：=函数[m，b[m，函数[r，r+1-Mod[r，2]][楼层@平方米[m] ]，4]][8n+4]；

%t a/@Range[0100]（*_Jean-François Alcover_，2021年5月30日，以_Alois P.Heinz_*命名）

%Y参见A316834、A316489、A316490。

%K nonn公司

%0、20

%A _乔恩·肖恩菲尔德，2018年7月29日