This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A300672 #8 Feb 11 2021 11:21:40 %S A300672 1,0,1,1,3,3,8,10,23,32,64,98,187,296,543,891,1595,2660,4694,7924, %T A300672 13854,23556,40940,69939,121122,207490,358517,615292,1061635,1824013, %U A300672 3144404,5406257,9314645,16021922,27595176,47478950,81757104,140691461,242232918,416890645,717712748,1235289624 %N A300672 Expansion of 1/(1 - Sum_{p prime, k>=1} x^(p^k)/(1 - x^(p^k))). %C A300672 Invert transform of A001222. %H A300672 Alois P. Heinz, <a href="/A300672/b300672.txt">Table of n, a(n) for n = 0..2000</a> %H A300672 N. J. A. Sloane, <a href="/transforms.txt">Transforms</a> %F A300672 G.f.: 1/(1 - Sum_{k>=2} A001222(k)*x^k). %p A300672 a:= proc(n) option remember; `if`(n=0, 1, %p A300672 add(a(n-i)*numtheory[bigomega](i), i=1..n)) %p A300672 end: %p A300672 seq(a(n), n=0..42); # _Alois P. Heinz_, Feb 11 2021 %t A300672 nmax = 41; CoefficientList[Series[1/(1 - Sum[Boole[PrimePowerQ[k]] x^k/(1 - x^k), {k, 1, nmax}]), {x, 0, nmax}], x] %t A300672 nmax = 41; CoefficientList[Series[1/(1 - Sum[PrimeOmega[k] x^k, {k, 2, nmax}]), {x, 0, nmax}], x] %t A300672 a[0] = 1; a[n_] := a[n] = Sum[PrimeOmega[k] a[n - k], {k, 1, n}]; Table[a[n], {n, 0, 41}] %Y A300672 Cf. A001222, A022559, A112967, A129921, A180305, A293549, A300671. %K A300672 nonn %O A300672 0,5 %A A300672 _Ilya Gutkovskiy_, Mar 11 2018