|
%I #16 Apr 06 2022 11:46:26
%S 1,1,0,-2,-3,-2,3,8,8,-2,-16,-24,-10,24,59,54,-11,-117,-174,-90,162,
%T 431,449,-20,-835,-1393,-848,1062,3352,3748,317,-6257,-11134,-7583,
%U 7294,25956,30786,5217,-46545,-88132,-65062,48534,199234,249263,63034,-342174,-691679,-554002
%N Expansion of 1/(1 - Sum_{k>=1} mu(k)*x^k), where mu() is the Moebius function (A008683).
%C Invert transform of A008683.
%H Seiichi Manyama, <a href="/A300663/b300663.txt">Table of n, a(n) for n = 0..5000</a>
%H N. J. A. Sloane, <a href="/transforms.txt">Transforms</a>
%F G.f.: 1/(1 - Sum_{k>=1} A008683(k)*x^k).
%F a(0) = 1; a(n) = Sum_{k=1..n} mu(k) * a(n-k). - _Seiichi Manyama_, Apr 06 2022
%p a:= proc(n) option remember; `if`(n=0, 1, add(
%p numtheory[mobius](j)*a(n-j), j=1..n))
%p end:
%p seq(a(n), n=0..50); # _Alois P. Heinz_, Mar 10 2018
%t nmax = 47; CoefficientList[Series[1/(1 - Sum[MoebiusMu[k] x^k, {k, 1, nmax}]), {x, 0, nmax}], x]
%t a[0] = 1; a[n_] := a[n] = Sum[MoebiusMu[k] a[n - k], {k, 1, n}]; Table[a[n], {n, 0, 47}]
%o (PARI) my(N=66, x='x+O('x^N)); Vec(1/(1-sum(k=1, N, moebius(k)*x^k))) \\ _Seiichi Manyama_, Apr 06 2022
%o (PARI) a(n) = if(n==0, 1, sum(k=1, n, moebius(k)*a(n-k))); \\ _Seiichi Manyama_, Apr 06 2022
%Y Cf. A008683, A050385, A068341, A073776, A073777, A117209, A117210, A185694, A280194.
%K sign
%O 0,4
%A _Ilya Gutkovskiy_, Mar 10 2018
|
