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.

Showing 1-3 of 3 results.

A274804 The exponential transform of sigma(n).

Original entry on oeis.org

1, 1, 4, 14, 69, 367, 2284, 15430, 115146, 924555, 7991892, 73547322, 718621516, 7410375897, 80405501540, 914492881330, 10873902417225, 134808633318271, 1738734267608613, 23282225008741565, 323082222240744379, 4638440974576329923, 68794595993688306903
Offset: 0

Views

Author

Johannes W. Meijer, Jul 27 2016

Keywords

Comments

The exponential transform [EXP] transforms an input sequence b(n) into the output sequence a(n). The EXP transform is the inverse of the logarithmic transform [LOG], see the Weisstein link and the Sloane and Plouffe reference. This relation goes by the name of Riddell's formula. For information about the logarithmic transform see A274805. The EXP transform is related to the multinomial transform, see A274760 and the second formula.
The definition of the EXP transform, see the second formula, shows that n >= 1. To preserve the identity LOG[EXP[b(n)]] = b(n) for n >= 0 for a sequence b(n) with offset 0 the shifted sequence b(n-1) with offset 1 has to be used as input for the exponential transform, otherwise information about b(0) will be lost in transformation.
In the a(n) formulas, see the examples, the multinomial coefficients A178867 appear.
We observe that a(0) = 1 and provides no information about any value of b(n), this notwithstanding it is customary to start the a(n) sequence with a(0) = 1.
The Maple programs can be used to generate the exponential transform of a sequence. The first program uses a formula found by Alois P. Heinz, see A007446 and the first formula. The second program uses the definition of the exponential transform, see the Weisstein link and the second formula. The third program uses information about the inverse of the exponential transform, see A274805.
Some EXP transform pairs are, n >= 1: A000435(n) and A065440(n-1); 1/A000027(n) and A177208(n-1)/A177209(n-1); A000670(n) and A075729(n-1); A000670(n-1) and A014304(n-1); A000045(n) and A256180(n-1); A000290(n) and A033462(n-1); A006125(n) and A197505(n-1); A053549(n) and A198046(n-1); A000311(n) and A006351(n); A030019(n) and A134954(n-1); A038048(n) and A053529(n-1); A193356(n) and A003727(n-1).

Examples

			Some a(n) formulas, see A178867:
a(0) = 1
a(1) = x(1)
a(2) = x(1)^2 + x(2)
a(3) = x(1)^3 + 3*x(1)*x(2) + x(3)
a(4) = x(1)^4 + 6*x(1)^2*x(2) + 4*x(1)*x(3) + 3*x(2)^2 + x(4)
a(5) = x(1)^5 + 10*x(1)^3*x(2) + 10*x(1)^2*x(3) + 15*x(1)*x(2)^2 + 5*x(1)*x(4) + 10*x(2)*x(3) + x(5)

References

  • Frank Harary and Edgar M. Palmer, Graphical Enumeration, 1973.
  • Robert James Riddell, Contributions to the theory of condensation, Dissertation, University of Michigan, Ann Arbor, 1951.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, 1995, pp. 18-23.

Links

Crossrefs

Cf. A274805, A274760, A178867, A000203, A007446.
Cf. A177208, A177209, A006351, A197505, A144180, A256180, A033462, A198046, A134954, A145460, A188489, A005432, A029725, A124213, A002801.
Cf. A036040, another version.

Programs

  • Maple
    nmax:=21: with(numtheory): b := proc(n): sigma(n) end: a:= proc(n) option remember; if n=0 then 1 else add(binomial(n-1, j-1) * b(j) *a(n-j), j=1..n) fi: end: seq(a(n), n=0..nmax); # End first EXP program.
nmax:= 21: with(numtheory): b := proc(n): sigma(n) end: t1 := exp(add(b(n)*x^n/n!, n=1..nmax+1)): t2 := series(t1, x, nmax+1): a := proc(n): n!*coeff(t2, x, n) end: seq(a(n), n=0..nmax); # End second EXP program.
nmax:=21: with(numtheory): b := proc(n): sigma(n) end: f := series(log(1+add(q(n)*x^n/n!, n=1..nmax+1)), x, nmax+1): d := proc(n): n!*coeff(f, x, n) end: a(0):=1: q(0):=1: a(1):=b(1): q(1):=b(1): for n from 2 to nmax+1 do q(n) := solve(d(n)-b(n), q(n)): a(n):=q(n): od: seq(a(n), n=0..nmax); # End third EXP program.
  • Mathematica
    a[0] = 1; a[n_] := a[n] = Sum[Binomial[n-1, j-1]*DivisorSigma[1, j]*a[n-j], {j, 1, n}]; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Feb 22 2017 *)
nmax = 20; CoefficientList[Series[Exp[Sum[DivisorSigma[1, k]*x^k/k!, {k, 1, nmax}]], {x, 0, nmax}], x] * Range[0, nmax]! (* Vaclav Kotesovec, Jun 08 2021 *)

Formula

a(n) = Sum_{j=1..n} (binomial(n-1,j-1) * b(j) * a(n-j)), n >= 1 and a(0) = 1, with b(n) = A000203(n) = sigma(n).
E.g.f.: exp(Sum_{n >= 1} b(n)*x^n/n!) with b(n) = sigma(n) = A000203(n).

A274844 The inverse multinomial transform of A001818(n) = ((2*n-1)!!)^2.

Original entry on oeis.org

1, 8, 100, 1664, 34336, 843776, 24046912, 779780096, 28357004800, 1143189536768, 50612287301632, 2441525866790912, 127479926768287744, 7163315850315825152, 431046122080208896000, 27655699473265974050816, 1884658377677216933085184
Offset: 1

Views

Author

Johannes W. Meijer, Jul 27 2016

Keywords

Comments

The inverse multinomial transform [IML] transforms an input sequence b(n) into the output sequence a(n). The IML transform inverses the effect of the multinomial transform [MNL], see A274760, and is related to the logarithmic transform, see A274805 and the first formula.
To preserve the identity MNL[IML[b(n)]] = b(n) for n >= 0 for a sequence b(n) with offset 0 the shifted sequence b(n-1) with offset 1 has to be used as input for the MNL.
In the a(n) formulas, see the examples, the cumulant expansion numbers A127671 appear.
We observe that the inverse multinomial transform leaves the value of a(0) undefined.
The Maple programs can be used to generate the inverse multinomial transform of a sequence. The first program is derived from a formula given by Alois P. Heinz for the logarithmic transform, see the first formula and A001187. The second program uses the e.g.f. for multivariate row polynomials, see A127671 and the examples. The third program uses information about the inverse of the inverse of the multinomial transform, see A274760.
The IML transform of A001818(n) = ((2*n-1)!!)^2 leads quite unexpectedly to A005411(n), a sequence related to certain Feynman diagrams.
Some IML transform pairs, n >= 1: A000110(n) and 1/A000142(n-1); A137341(n) and A205543(n); A001044(n) and A003319(n+1); A005442(n) and A000204(n); A005443(n) and A001350(n); A007559(n) and A000244(n-1); A186685(n+1) and A131040(n-1); A061711(n) and A141151(n); A000246(n) and A000035(n); A001861(n) and A141044(n-1)/A001710(n-1); A002866(n) and A000225(n); A000262(n) and A000027(n).

Examples

			Some a(n) formulas, see A127671:
a(0) = undefined
a(1) = (1/0!) * (1*x(1))
a(2) = (1/1!) * (1*x(2) - x(1)^2)
a(3) = (1/2!) * (1*x(3) - 3*x(2)*x(1) + 2*x(1)^3)
a(4) = (1/3!) * (1*x(4) - 4*x(3)*x(1) - 3*x(2)^2 + 12*x(2)*x(1)^2 - 6*x(1)^4)
a(5) = (1/4!) * (1* x(5) - 5*x(4)*x(1) - 10*x(3)*x(2) + 20*x(3)*x(1)^2 + 30*x(2)^2*x(1) -60*x(2)*x(1)^3 + 24*x(1)^5)

References

  • Richard P. Feynman, QED, The strange theory of light and matter, 1985.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, 1995, pp. 18-23.

Links

Crossrefs

Cf. A274760, A274805, A127671, A001187, A001818, A000079, A005411.

Programs

  • Maple
    nmax:=17: b := proc(n): (doublefactorial(2*n-1))^2 end: c:= proc(n) option remember; b(n) - add(k*binomial(n, k)*b(n-k)*c(k), k=1..n-1)/n end: a := proc(n): c(n)/(n-1)! end: seq(a(n), n=1..nmax); # End first IML program.
nmax:=17: b := proc(n): (doublefactorial(2*n-1))^2 end: t1 := log(1+add(b(n)*x^n/n!, n=1..nmax+1)): t2 := series(t1, x, nmax+1): a := proc(n): n*coeff(t2, x, n) end: seq(a(n), n=1..nmax); # End second IML program.
nmax:=17: b := proc(n): (doublefactorial(2*n-1))^2 end: f := series(exp(add(t(n)*x^n/n, n=1..nmax)), x, nmax+1): d := proc(n): n!*coeff(f, x, n) end: a(1):=b(1): t(1):= b(1): for n from 2 to nmax+1 do t(n) := solve(d(n)-b(n), t(n)): a(n):=t(n): od: seq(a(n), n=1..nmax); # End third IML program.
  • Mathematica
    nMax = 22; b[n_] := ((2*n-1)!!)^2; c[n_] := c[n] = b[n] - Sum[k*Binomial[n, k]*b[n-k]*c[k], {k, 1, n-1}]/n; a[n_] := c[n]/(n-1)!; Table[a[n], {n, 1, nMax}] (* Jean-François Alcover, Feb 27 2017, translated from Maple *)

Formula

a(n) = c(n)/(n-1)! with c(n) = b(n) - Sum_{k=1..n-1}(k*binomial(n, k)*b(n-k)*c(k)), n >= 1 and a(0) = undefined, with b(n) = A001818(n) = ((2*n-1)!!)^2.
a(n) = A000079(n-1) * A005411(n), n >= 1.

A322513 Expansion of e.g.f. log(1 + Sum_{k>=1} d(k) * x^k / k!), where d(k) = number of divisors of k (A000005).

Original entry on oeis.org

0, 1, 1, -2, 1, 11, -48, -6, 1241, -6431, -15320, 452970, -2317212, -17584137, 372119776, -1552313624, -31732274313, 565880016193, -1217992446564, -90197542736656, 1400682677566587, 1990004001731140, -384348195167184028, 5109122826021406702
Offset: 0

Views

Author

Ilya Gutkovskiy, Oct 03 2019

Keywords

Comments

Logarithmic transform of A000005.

Links

Crossrefs

Cf. A000005, A028342, A274805, A294363, A295739.

Programs

  • Maple
    a:= proc(n) option remember; `if`(n=0, 0, (b-> b(n)-add(a(j)
     *binomial(n, j)*j*b(n-j), j=1..n-1)/n)(numtheory[tau]))
    end:
seq(a(n), n=0..25);  # Alois P. Heinz, Oct 06 2019
  • Mathematica
    nmax = 23; CoefficientList[Series[Log[1 + Sum[DivisorSigma[0, k] x^k/k!, {k, 1, nmax}]], {x, 0, nmax}], x] Range[0, nmax]!
a[n_] := a[n] = DivisorSigma[0, n] - Sum[Binomial[n, k] DivisorSigma[0, n - k] k a[k], {k, 1, n - 1}]/n; a[0] = 0; Table[a[n], {n, 0, 23}]
Showing 1-3 of 3 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.