login
Number of ways to write n as an ordered sum of 4 squarefree numbers.
8

%I #6 Feb 04 2021 17:19:52

%S 1,4,10,16,23,32,50,68,83,92,116,148,178,192,224,276,335,360,400,460,

%T 547,580,634,704,821,868,938,1024,1162,1212,1288,1392,1572,1628,1742,

%U 1876,2123,2172,2308,2460,2761,2820,2964,3176,3550,3628,3778,4028,4481,4528,4686,4932,5513,5564

%N Number of ways to write n as an ordered sum of 4 squarefree numbers.

%F G.f.: (Sum_{k>=1} mu(k)^2 * x^k)^4.

%p b:= proc(n, t) option remember;

%p `if`(n=0, `if`(t=0, 1, 0), `if`(t<1, 0, add(

%p `if`(numtheory[issqrfree](j), b(n-j, t-1), 0), j=1..n)))

%p end:

%p a:= n-> b(n, 4):

%p seq(a(n), n=4..57); # _Alois P. Heinz_, Feb 04 2021

%t nmax = 57; CoefficientList[Series[Sum[MoebiusMu[k]^2 x^k, {k, 1, nmax}]^4, {x, 0, nmax}], x] // Drop[#, 4] &

%Y Cf. A005117, A008683, A008966, A098235, A280210, A308767, A341065, A341066, A341067, A341068, A341069, A341070.

%K nonn

%O 4,2

%A _Ilya Gutkovskiy_, Feb 04 2021