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.

A003239 Number of rooted planar trees with n non-root nodes: circularly cycling the subtrees at the root gives equivalent trees.

Original entry on oeis.org

1, 1, 2, 4, 10, 26, 80, 246, 810, 2704, 9252, 32066, 112720, 400024, 1432860, 5170604, 18784170, 68635478, 252088496, 930138522, 3446167860, 12815663844, 47820447028, 178987624514, 671825133648, 2528212128776, 9536895064400, 36054433810102, 136583761444364, 518401146543812
Offset: 0

Views

Author

Keywords

Comments

Also number of necklaces with 2*n beads, n white and n black (to get the correspondence, start at root, walk around outside of tree, use white if move away from the root, black if towards root).
Also number of terms in polynomial expression for permanent of generic circulant matrix of order n.
a(n) is the number of equivalence classes of n-compositions of n under cyclic rotation. (Given a necklace, split it into runs of white followed by a black bead and record the lengths of the white runs. This gives an n-composition of n.) a(n) is the number of n-multisets in Z mod n whose sum is 0. - David Callan, Nov 05 2003
a(n) is the number of cyclic equivalence classes of triangulations of a once-punctured n-gon. - Esther Banaian, May 06 2025

Examples

			As _David Callan_ said, a(n) is the number of n-multisets in Z mod n whose sum is 0. So for n = 4 the a(4)=10 multisets are (0, 0, 0, 0), (1, 1, 1, 1), (0, 1, 1, 2), (0, 0, 2, 2), (2, 2, 2, 2), (0, 0, 1, 3), (1, 2, 2, 3), (1, 1, 3, 3), (0, 2, 3, 3) and (3, 3, 3, 3). - _Boas Bakker_, Apr 21 2025
		

References

  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 305 (see R(x)).
  • F. Harary and E. M. Palmer, Graphical Enumeration, Academic Press, NY, 1973; page 80, Problem 3.13.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Problem 7.112(b).

Crossrefs

Column k=2 of A208183.
Column k=1 of A261494.

Programs

  • Maple
    with(numtheory): A003239 := proc(n) local t1,t2,d; t2 := divisors(n); t1 := 0; for d in t2 do t1 := t1+phi(n/d)*binomial(2*d,d)/(2*n); od; t1; end;
    spec := [ C, {B=Union(Z,Prod(B,B)), C=Cycle(B)}, unlabeled ]; [seq(combstruct[count](spec, size=n), n=0..40)];
  • Mathematica
    a[n_] := Sum[ EulerPhi[n/k]*Binomial[2k, k]/(2n), {k, Divisors[n]}]; a[0] = 1; Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Apr 11 2012 *)
  • PARI
    C(n, k)=binomial(n,k);
    a(n) = if(n<=0, n==0, sumdiv(n, d, eulerphi(n/d) * C(2*d, d)) / (2*n) );
    /* or, second formula: */
    /* a(n) = if(n<=0, n==0, sumdiv(n, d, eulerphi(n/d) * C(2*d-1,d)) / n ); */
    /* Joerg Arndt, Oct 21 2012 */
    
  • SageMath
    def A003239(n):
        if n == 0: return 1
        return sum(euler_phi(n/d)*binomial(2*d, d)/(2*n) for d in divisors(n))
    print([A003239(n) for n in (0..29)]) # Peter Luschny, Dec 10 2020

Formula

a(n) = Sum_{d|n} (phi(n/d)*binomial(2*d, d))/(2*n) for n > 0.
a(n) = (1/n)*Sum_{d|n} (phi(n/d)*binomial(2*d-1, d)) for n > 0.
a(n) = A047996(2*n, n). - Philippe Deléham, Jul 25 2006
a(n) ~ 2^(2*n-1) / (sqrt(Pi) * n^(3/2)). - Vaclav Kotesovec, Aug 22 2015

Extensions

Sequence corrected and extended by Roderick J. Fletcher (yylee(AT)mail.ncku.edu.tw), Aug 1997
Additional comments from Michael Somos

A261494 Number A(n,k) of necklaces with n white beads and k*n black beads; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 4, 1, 1, 1, 4, 10, 10, 1, 1, 1, 5, 19, 43, 26, 1, 1, 1, 6, 31, 116, 201, 80, 1, 1, 1, 7, 46, 245, 776, 1038, 246, 1, 1, 1, 8, 64, 446, 2126, 5620, 5538, 810, 1, 1, 1, 9, 85, 735, 4751, 19811, 42288, 30667, 2704, 1
Offset: 0

Views

Author

Alois P. Heinz, Aug 21 2015

Keywords

Comments

For k>=1 is column k asymptotic to (k+1)^((k+1)*n-1/2) / (sqrt(2*Pi) * k^(k*n+1/2) * n^(3/2)). - Vaclav Kotesovec, Aug 22 2015

Examples

			A(2,2) = 3: 000011, 000101, 001001.
A(3,2) = 10: 000000111, 000001011, 000010011, 000100011, 001000011, 010000011, 000010101, 000100101, 001000101, 001001001.
Square array A(n,k) begins:
  1,  1,    1,    1,     1,     1,      1, ...
  1,  1,    1,    1,     1,     1,      1, ...
  1,  2,    3,    4,     5,     6,      7, ...
  1,  4,   10,   19,    31,    46,     64, ...
  1, 10,   43,  116,   245,   446,    735, ...
  1, 26,  201,  776,  2126,  4751,   9276, ...
  1, 80, 1038, 5620, 19811, 54132, 124936, ...
		

Crossrefs

Main diagonal gives A261495.
Lower diagonal gives A261496.
Cf. A000010.

Programs

  • Maple
    with(numtheory):
    A:= (n, k)-> `if`(n=0, 1, add(binomial((k+1)*n/d, n/d)
                        *phi(d), d=divisors(n))/((k+1)*n)):
    seq(seq(A(n, d-n), n=0..d), d=0..14);
  • Mathematica
    A[n_, k_] := If[n==0, 1, DivisorSum[n, Binomial[(k+1)*n/#, n/#]*EulerPhi[#] /((k+1)*n)&]]; Table[A[n, d-n], {d, 0, 14}, {n, 0, d}] // Flatten (* Jean-François Alcover, Feb 19 2017, translated from Maple *)
  • PARI
    a(n,k) = if(n<1, 1, sumdiv(n, d, binomial((k + 1)*n/d, n/d) * eulerphi(d)) / ((k + 1)*n));
    for(d=0, 14, for(n=0, d, print1(a(n, d - n),", ");); print();) \\ Indranil Ghosh, Mar 25 2017

Formula

A(n,k) = 1/((k+1)*n) * Sum_{d|n} C((k+1)*n/d,n/d) * A000010(d) for n>0, A(0,k) = 1.
A(n,k) = 1/((k+1)*n)*Sum_{i=1..n} C((k+1)*gcd(n,i),gcd(n,i)) = 1/((k+1)*n)*Sum_{i=1..n} C((k+1)*n/gcd(n,i),n/gcd(n,i))*phi(gcd(n,i))/phi(n/gcd(n,i)) for n >= 1, where phi = A000010. - Richard L. Ollerton, May 19 2021

A346577 a(n) = (1/(3*n)) * Sum_{d|n} mu(n/d) * binomial(3*d,d).

Original entry on oeis.org

1, 2, 9, 40, 200, 1026, 5537, 30624, 173583, 1001400, 5864749, 34768296, 208267319, 1258574114, 7663720500, 46976003712, 289628805622, 1794932293950, 11175157356521, 69864074596000, 438403736543598, 2760351027094298, 17433869214973753, 110420300844952992
Offset: 1

Views

Author

Ilya Gutkovskiy, Jul 24 2021

Keywords

Comments

Inverse Euler transform of A001764.
Moebius transform of A082936.

Crossrefs

Programs

  • Mathematica
    Table[(1/(3 n)) Sum[MoebiusMu[n/d] Binomial[3 d, d], {d, Divisors[n]}], {n, 24}]
  • PARI
    a(n) = sumdiv(n, d, moebius(n/d)*binomial(3*d,d))/(3*n); \\ Michel Marcus, Jul 24 2021

Formula

a(n) = A060170(n)/3. - Hugo Pfoertner, Jul 24 2021
Showing 1-3 of 3 results.