OFFSET
1,2
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..1000
Eric Weisstein's World of Mathematics, Root-Mean-Square
EXAMPLE
a(14) = 6 subsets: {1}, {2}, {7}, {14}, {1, 7} and {2, 14}.
MAPLE
a:= proc(n) option remember; uses numtheory; local l, b;
l, b:= sort([divisors(n)[]]),
proc(i, s, c) option remember;
`if`(i=0, `if`(c>0 and issqr(s/c), 1, 0),
b(i-1, s, c)+b(i-1, s+l[i]^2, c+1))
end; forget(b); b(nops(l), 0$2)
end:
seq(a(n), n=1..100); # Alois P. Heinz, Sep 30 2022
MATHEMATICA
a[n_] := a[n] = Module[{b, l = Divisors[n]}, b[i_, s_, c_] := b[i, s, c] = If[i == 0, If[c > 0 && IntegerQ @ Sqrt[s/c], 1, 0], b[i-1, s, c]+b[i-1, s+l[[i]]^2, c+1]]; b[Length[l], 0, 0]];
Table[a[n], {n, 1, 100}] (* Jean-François Alcover, Oct 13 2022, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Dec 11 2020
STATUS
approved