cp's OEIS Frontend

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.

A318290 Expansion of Sum_{k>=1} (-1 + Product_{j>=1} (1 + j*x^(k*j))).

This page as a plain text file.
%I A318290 #9 Apr 03 2019 02:59:03
%S A318290 1,3,6,10,16,33,44,74,126,204,289,503,696,1151,1749,2599,3742,5928,
%T A318290 8245,12658,18351,26715,37828,55296,78346,111882,159664,226782,315416,
%U A318290 446670,618667,860764,1199995,1649820,2289020,3157349,4303996,5878786,8033272,10894516,14749052
%N A318290 Expansion of Sum_{k>=1} (-1 + Product_{j>=1} (1 + j*x^(k*j))).
%C A318290 Inverse Moebius transform of A022629.
%H A318290 N. J. A. Sloane, <a href="/transforms.txt">Transforms</a>
%F A318290 G.f.: Sum_{k>=1} A022629(k)*x^k/(1 - x^k).
%F A318290 a(n) = Sum_{d|n} A022629(d).
%p A318290 a:=series(add(-1+mul(1+j*x^(k*j),j=1..100),k=1..100),x=0,42): seq(coeff(a,x,n),n=1..41); # _Paolo P. Lava_, Apr 02 2019
%t A318290 nmax = 41; Rest[CoefficientList[Series[Sum[-1 + Product[(1 + j x^(k j)), {j, 1, nmax}], {k, 1, nmax}], {x, 0, nmax}], x]]
%t A318290 b[n_] := b[n] = SeriesCoefficient[Product[(1 + k x^k), {k, 1, n}], {x, 0, n}]; a[n_] := a[n] = SeriesCoefficient[Sum[b[k] x^k/(1 - x^k), {k, 1, n}], {x, 0, n}]; Table[a[n], {n, 41}]
%t A318290 b[0] = 1; b[n_] := b[n] = Sum[Sum[(-d)^(k/d + 1), {d, Divisors[k]}] b[n - k], {k, 1, n}]/n; a[n_] := a[n] = Sum[b[d], {d, Divisors[n]}]; Table[a[n], {n, 41}]
%Y A318290 Cf. A022629, A047966, A300278, A302550, A318025.
%K A318290 nonn
%O A318290 1,2
%A A318290 _Ilya Gutkovskiy_, Aug 23 2018