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.

A022553 Number of binary Lyndon words containing n letters of each type; periodic binary sequences of period 2n with n zeros and n ones in each period.

Original entry on oeis.org

1, 1, 1, 3, 8, 25, 75, 245, 800, 2700, 9225, 32065, 112632, 400023, 1432613, 5170575, 18783360, 68635477, 252085716, 930138521, 3446158600, 12815663595, 47820414961, 178987624513, 671825020128, 2528212128750, 9536894664375, 36054433807398, 136583760011496
Offset: 0

Views

Author

Keywords

Comments

Also number of asymmetric rooted plane trees with n+1 nodes. - Christian G. Bower
Conjecturally, number of irreducible alternating Euler sums of depth n and weight 3n.
a(n+1) is inverse Euler transform of A000108. Inverse Witt transform of A006177.
Dimension of the degree n part of the primitive Lie algebra of the Hopf algebra CQSym (Catalan Quasi-Symmetric functions). - Jean-Yves Thibon (jyt(AT)univ-mlv.fr), Oct 22 2006
For n>0, 2*a(n) is divisible by n (cf. A268619), 12*a(n) is divisible by n^2 (cf. A268592). - Max Alekseyev, Feb 09 2016

Examples

			a(3)=3 counts 6-periodic 000111, 001011 and 001101. a(4)=8 counts 00001111, 00010111, 00011011, 00011101, 00100111, 00101011, 00101101, and 00110101. - _R. J. Mathar_, Oct 20 2021
		

References

  • F. Bergeron, G. Labelle and P. Leroux, Combinatorial Species and Tree-Like Structures, Cambridge, 1998, p. 336 (4.4.64)

Crossrefs

Cf. A003239, A005354, A000740, A007727, A086655, A289978 (multiset trans.), A001037 (binary Lyndon wds.), A074655 (3 letters), A074656 (4 letters).
A diagonal of the square array described in A051168.

Programs

  • Maple
    with(numtheory):
    a:= n-> `if`(n=0, 1,
            add(mobius(n/d)*binomial(2*d, d), d=divisors(n))/(2*n)):
    seq(a(n), n=0..30);  # Alois P. Heinz, Jan 21 2011
  • Mathematica
    a[n_] := Sum[MoebiusMu[n/d]*Binomial[2d, d], {d, Divisors[n]}]/(2n); a[0] = 1; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Feb 02 2015 *)
  • PARI
    a(n)=if(n<1,n==0,sumdiv(n,d,moebius(n/d)*binomial(2*d,d))/2/n)
    
  • Python
    from sympy import mobius, binomial, divisors
    def a(n):
        return 1 if n == 0 else sum(mobius(n//d)*binomial(2*d, d) for d in divisors(n))//(2*n)
    print([a(n) for n in range(31)]) # Indranil Ghosh, Aug 05 2017
    
  • Sage
    def a(n):
        return 1 if n ==0 else sum(moebius(n//d)*binomial(2*d, d) for d in divisors(n))//(2*n)
    # F. Chapoton, Apr 23 2020

Formula

a(n) = A060165(n)/2 = A007727(n)/(2*n) = A045630(n)/n.
Product_n (1-x^n)^a(n) = 2/(1+sqrt(1-4*x)); a(n) = 1/(2*n) * Sum_{d|n} mu(n/d)*C(2*d,d). Also Moebius transform of A003239. - Christian G. Bower
a(n) ~ 2^(2*n-1) / (sqrt(Pi) * n^(3/2)). - Vaclav Kotesovec, Sep 11 2014
G.f.: 1 + Sum_{k>=1} mu(k)*log((1 - sqrt(1 - 4*x^k))/(2*x^k))/k. - Ilya Gutkovskiy, May 18 2019

A292287 Number of multisets of exactly n nonempty balanced binary Lyndon words with a total of 4n letters (2n zeros and 2n ones).

Original entry on oeis.org

1, 1, 4, 12, 43, 142, 508, 1781, 6414, 23124, 84296, 308613, 1137129, 4207456, 15636927, 58322808, 218272766, 819319778, 3083913810, 11636761924, 44010780075, 166802192488, 633420816341, 2409731688860, 9182826866499, 35048239457878, 133965833871427
Offset: 0

Views

Author

Alois P. Heinz, Sep 20 2017

Keywords

Crossrefs

Programs

  • Maple
    with(numtheory):
    g:= proc(n) option remember; `if`(n=0, 1, add(
          mobius(n/d)*binomial(2*d, d), d=divisors(n))/(2*n))
        end:
    a:= proc(n) option remember; `if`(n=0, 1, add(add(
          d*g(d+1), d=divisors(j))*a(n-j), j=1..n)/n)
        end:
    seq(a(n), n=0..30);
  • Mathematica
    g[n_] := g[n] = If[n == 0, 1, Sum[MoebiusMu[n/d] Binomial[2d, d], {d, Divisors[n]}]/(2n)];
    a[n_] := a[n] = If[n == 0, 1, Sum[Sum[d*g[d+1], {d, Divisors[j]}]*a[n-j], {j, 1, n}]/n];
    Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Nov 23 2023, after Alois P. Heinz *)

Formula

G.f.: Product_{j>=1} 1/(1-x^j)^A022553(j+1).
a(n) = A289978(2n,n).
Showing 1-2 of 2 results.