OFFSET
0,79
LINKS
FORMULA
a(n) = [x^n y^4] Product_{k>=1} (1 + y*x^(k^2)). - Ilya Gutkovskiy, Apr 22 2019
MAPLE
b:= proc(n, i, t) option remember; `if`(n=0, `if`(t=0, 1, 0),
`if`(t*i^2<n, 0, `if`(i=1, 0, b(n, i-1, t))+
`if`(i^2>n, 0, b(n-i^2, i-1, t-1))))
end:
a:= n-> b(n, isqrt(n), 4):
seq(a(n), n=0..150); # Alois P. Heinz, Feb 07 2013
MATHEMATICA
b[n_, i_, t_] := b[n, i, t] = If[n==0, If[t==0, 1, 0], If[t*i^2<n, 0, If[i == 1, 0, b[n, i-1, t]] + If[i^2>n, 0, b[n-i^2, i-1, t-1]]]]; a[n_] := b[n, Sqrt[n] // Floor, 4]; Table[a[n], {n, 0, 150}] (* Jean-François Alcover, Feb 29 2016, after Alois P. Heinz*)
dnzs[n_]:=Length[Select[IntegerPartitions[n, {4}], Length[Union[#]]==4&&AllTrue[ Sqrt[ #], IntegerQ] && FreeQ[#, 0]&]]; Array[dnzs, 110, 0] (* Harvey P. Dale, Jun 09 2024 *)
CROSSREFS
KEYWORD
nonn,look
AUTHOR
STATUS
approved