login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A340903 a(0) = 1; a(n) = Sum_{k=1..n} binomial(n,k) * sigma_0(k) * a(n-k). 7
1, 1, 4, 20, 139, 1192, 12318, 148318, 2041754, 31616757, 544005172, 10296204096, 212589150300, 4755177958104, 114545293676588, 2956316416222300, 81386676426000157, 2380590235918735576, 73729207700492304684, 2410324868012471929670, 82944575892433740648996 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,3
LINKS
FORMULA
E.g.f.: 1 / (1 - Sum_{i>=1} Sum_{j>=1} x^(i*j) / (i*j)!).
E.g.f.: 1 / (1 - Sum_{k>=1} sigma_0(k) * x^k/k!). - Seiichi Manyama, Mar 29 2022
MATHEMATICA
a[0] = 1; a[n_] := a[n] = Sum[Binomial[n, k] DivisorSigma[0, k] a[n - k], {k, 1, n}]; Table[a[n], {n, 0, 20}]
nmax = 20; CoefficientList[Series[1/(1 - Sum[Sum[x^(i j)/(i j)!, {j, 1, nmax}], {i, 1, nmax}]), {x, 0, nmax}], x] Range[0, nmax]!
PROG
(PARI) my(N=40, x='x+O('x^N)); Vec(serlaplace(1/(1-sum(k=1, N, numdiv(k)*x^k/k!)))) \\ Seiichi Manyama, Mar 29 2022
(PARI) a(n) = if(n==0, 1, sum(k=1, n, numdiv(k)*binomial(n, k)*a(n-k))); \\ Seiichi Manyama, Mar 29 2022
CROSSREFS
Sequence in context: A244756 A292932 A187116 * A098541 A004204 A160567
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Jan 26 2021
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 10:55 EDT 2024. Contains 371241 sequences. (Running on oeis4.)