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-2 of 2 results.

A050385 Reversion of Moebius function A008683.

Original entry on oeis.org

1, 1, 3, 10, 39, 160, 691, 3081, 14095, 65757, 311695, 1496833, 7266979, 35608419, 175875537, 874698246, 4376646808, 22016578909, 111282845162, 564886771380, 2878498888625, 14719219809915, 75505990358779, 388451973679785
Offset: 1

Views

Author

Christian G. Bower, Nov 15 1999

Keywords

Comments

From David W. Wilson, May 17 2017: (Start)
Appears to equal the number of ways to partition Z into n residue classes. For example, a(4) = 10 since we can partition Z into 4 residue classes in 10 ways:
Z = ∪ {0 (mod 2), 1 (mod 4), 3 (mod 8), 7 (mod 8)}
Z = ∪ {0 (mod 2), 3 (mod 4), 1 (mod 8), 5 (mod 8)}
Z = ∪ {0 (mod 2), 1 (mod 6), 3 (mod 6), 5 (mod 6)}
Z = ∪ {1 (mod 2), 0 (mod 4), 2 (mod 8), 6 (mod 8)}
Z = ∪ {1 (mod 2), 2 (mod 4), 0 (mod 8), 4 (mod 8)}
Z = ∪ {1 (mod 2), 0 (mod 6), 2 (mod 6), 4 (mod 6)}
Z = ∪ {0 (mod 3), 1 (mod 3), 2 (mod 6), 5 (mod 6)}
Z = ∪ {0 (mod 3), 2 (mod 3), 1 (mod 6), 4 (mod 6)}
Z = ∪ {1 (mod 3), 2 (mod 3), 0 (mod 6), 3 (mod 6)}
Z = ∪ {0 (mod 4), 1 (mod 4), 2 (mod 4), 3 (mod 4)}
(End)
Unfortunately this conjecture is not correct; it fails at a(13). - Jeffrey Shallit, Nov 19 2017
The correct statement is that a(n) is the number of "natural exact covering systems" of cardinality n. These are covering systems (like the ones in David W. Wilson's comment) that are obtained by starting with the size-1 system x == 0 (mod 1) and successively choosing a congruence and "splitting" it into r >= 2 new congruences. - Jeffrey Shallit, Dec 07 2017

Crossrefs

Programs

  • Mathematica
    InverseSeries[Sum[MoebiusMu[n] x^n, {n, 0, 25}] + O[x]^25] // CoefficientList[#, x]& // Rest (* Jean-François Alcover, Sep 29 2018 *)

Formula

G.f. A(x) satisfies: A(x) = x - Sum_{k>=2} mu(k) * A(x)^k. - Ilya Gutkovskiy, Apr 22 2020

A300673 Expansion of e.g.f. exp(Sum_{k>=1} mu(k)*x^k/k!), where mu() is the Moebius function (A008683).

Original entry on oeis.org

1, 1, 0, -3, -6, 5, 61, 126, -308, -2772, -5669, 25630, 224730, 486551, -3068155, -29264219, -72173176, 513535711, 5625869262, 16687752839, -113740116822, -1496118902963, -5508392724427, 31534346503605, 523333047780288, 2414704077547660, -10254467367668159
Offset: 0

Views

Author

Ilya Gutkovskiy, Mar 11 2018

Keywords

Comments

Exponential transform of A008683.

Examples

			E.g.f.: A(x) = 1 + x/1! - 3*x^3/3! - 6*x^4/4! + 5*x^5/5! + 61*x^6/6! + 126*x^7/7! - 308*x^8/8! - 2772*x^9/9! - 5669*x^10/10! + ...
		

Crossrefs

Programs

  • Mathematica
    nmax = 26; CoefficientList[Series[Exp[Sum[MoebiusMu[k] x^k/k!, {k, 1, nmax}]], {x, 0, nmax}], x] Range[0, nmax]!
    a[n_] := a[n] = Sum[MoebiusMu[k] Binomial[n - 1, k - 1] a[n - k], {k, 1, n}]; a[0] = 1; Table[a[n], {n, 0, 26}]
  • PARI
    a(n) = if(n==0, 1, sum(k=1, n, moebius(k)*binomial(n-1, k-1)*a(n-k))); \\ Seiichi Manyama, Feb 27 2022

Formula

E.g.f.: exp(Sum_{k>=1} A008683(k)*x^k/k!).
a(0) = 1; a(n) = Sum_{k=1..n} mu(k) * binomial(n-1,k-1) * a(n-k). - Seiichi Manyama, Feb 27 2022
Showing 1-2 of 2 results.