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.

A007561 Number of asymmetric rooted connected graphs where every block is a complete graph.

Original entry on oeis.org

0, 1, 1, 1, 3, 6, 16, 43, 120, 339, 985, 2892, 8606, 25850, 78347, 239161, 734922, 2271085, 7054235, 22010418, 68958139, 216842102, 684164551, 2165240365, 6871792256, 21865189969, 69737972975, 222915760126, 714001019626, 2291298553660, 7366035776888
Offset: 0

Views

Author

Keywords

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Column k=2 of A316101.

Programs

  • Maple
    g:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          add(binomial(a(i), j)*g(n-i*j, i-1), j=0..n/i)))
        end:
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          add(binomial(g(i, i), j)*b(n-i*j, i-1), j=0..n/i)))
        end:
    a:= n-> `if`(n<1, 0, b(n-1, n-1)):
    seq(a(n), n=0..40); # Alois P. Heinz, May 19 2013
  • Mathematica
    g[n_, i_] := g[n, i] = If[n == 0, 1, If[i<1, 0, Sum[Binomial[a[i], j]*g[n-i*j, i-1], {j, 0, n/i}]]]; b[n_, i_] := b[n, i] = If[n == 0, 1, If[i<1, 0, Sum[Binomial[g[i, i], j]*b[n-i*j, i-1], {j, 0, n/i}]]]; a[n_] := If[n<1, 0, b[n-1, n-1]]; Table[a[n] // FullSimplify, {n, 0, 40}] (* Jean-François Alcover, Feb 11 2014, after Alois P. Heinz *)

Formula

Shifts left when weigh-transform applied twice.
a(n) ~ c * d^n / n^(3/2), where d = 3.382016466020272807429818743..., c = 0.161800727760188847021075748... . - Vaclav Kotesovec, Jul 26 2014

Extensions

Additional comments from Christian G. Bower