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-10 of 30 results. Next

A238357 Number of genus-7 rooted maps with n edges.

Original entry on oeis.org

14230536445125, 3128879373858000, 360626952084151500, 29001816720933903504, 1828003659229082834100, 96187365300257285300064, 4395215998078319892167640, 179153431308203084149883760, 6641365771586560905099092466, 227189907562197156785567456832, 7252879937219595844346639732688
Offset: 14

Views

Author

Joerg Arndt, Feb 26 2014

Keywords

Crossrefs

Column g=7 of A269919.
Cf. A239921 (unrooted sensed), A348800 (unrooted unsensed).
Rooted maps with n edges of genus g for 0 <= g <= 10: A000168, A006300, A006301, A104742, A215402, A238355, A238356, this sequence, A238358, A238359, A238360.

Programs

  • Mathematica
    T[0, 0] = 1; T[n_, g_] /; g < 0 || g > n/2 = 0; T[n_, g_] := T[n, g] = ((4 n - 2)/3 T[n - 1, g] + (2 n - 3) (2 n - 2) (2 n - 1)/12 T[n - 2, g - 1] + 1/2 Sum[(2 k - 1) (2 (n - k) - 1) T[k - 1, i] T[n - k - 1, g - i], {k, 1, n - 1}, {i, 0, g}])/((n + 1)/6);
    a[n_] := T[n, 7];
    Table[a[n], {n, 14, 30}] (* Jean-François Alcover, Jul 20 2018 *)
  • PARI
    \\ see A238396
    
  • PARI
    system("wget http://oeis.org/A238357/a238357.txt");
    A005159_ser(N) = my(x='x+O('x^(N+1))); (1 - sqrt(1-12*x))/(6*x);
    A238357_ser(N) = subst(read("a238357.txt"), 'y, A005159_ser(N+14));
    Vec(A238357_ser(11)) \\ Gheorghe Coserea, Jun 03 2017

A238358 Number of genus-8 rooted maps with n edges.

Original entry on oeis.org

11288163762500625, 2927974178219879250, 394372363395179602125, 36751560969705187643982, 2663973075006196131775590, 160098273686603663417293308, 8303278159618015743881266599, 381958851175370643701603049354, 15896435050196091382215375181044, 607566907750822335161584110201960
Offset: 16

Views

Author

Joerg Arndt, Feb 26 2014

Keywords

Crossrefs

Column g=8 of A269919.
Cf. A239922 (unrooted sensed), A348801 (unrooted unsensed).
Rooted maps with n edges of genus g for 0 <= g <= 10: A000168, A006300, A006301, A104742, A215402, A238355, A238356, A238357, this sequence, A238359, A238360.

Programs

  • Mathematica
    T[0, 0] = 1; T[n_, g_] /; g < 0 || g > n/2 = 0; T[n_, g_] := T[n, g] = ((4 n - 2)/3 T[n - 1, g] + (2 n - 3) (2 n - 2) (2 n - 1)/12 T[n - 2, g - 1] + 1/2 Sum[(2 k - 1) (2 (n - k) - 1) T[k - 1, i] T[n - k - 1, g - i], {k, 1, n - 1}, {i, 0, g}])/((n + 1)/6);
    a[n_] := T[n, 8];
    Table[a[n], {n, 16, 30}] (* Jean-François Alcover, Jul 20 2018 *)
  • PARI
    \\ see A238396

A053764 a(n) = 3^(n^2 - n).

Original entry on oeis.org

1, 1, 9, 729, 531441, 3486784401, 205891132094649, 109418989131512359209, 523347633027360537213511521, 22528399544939174411840147874772641, 8727963568087712425891397479476727340041449, 30432527221704537086371993251530170531786747066637049, 955004950796825236893190701774414011919935138974343129836853841, 269721605590607563262106870407286853611938890184108047911269431464974473521
Offset: 0

Views

Author

Stephen G Penrice, Mar 29 2000

Keywords

Comments

Number of nilpotent n X n matrices X over GF(3), that is, the number of n X n matrices X over GF(3) satisfying X^k = 0 for some k >= 1.
More generally, Fine and Herstein prove that the probability that an n X n matrix over GF(p^m) is nilpotent is 1/p^(mn) and the probability that an n X n matrix over Z/mZ is nilpotent is 1/k^n, where k is the product of the distinct prime factors of m.
Is this the same sequence (apart from the initial term) as A053854? - Philippe Deléham, Dec 09 2007
[1,9,729,531441,3486784401,...] is the Hankel transform of A005159. - Philippe Deléham, Dec 10 2007

Crossrefs

Programs

Formula

Sequence given by the Hankel transform (see A001906 for definition) of A082181 = {1, 1, 10, 109, 1270, 15562, 198100, ...}; example: det([1, 1, 10, 109; 1, 10, 109, 1270; 10, 109, 1270, 15562; 109, 1270, 15562, 198100]) = 9^6 = 531441. - Philippe Deléham, Aug 20 2005

Extensions

More terms from James Sellers, Apr 08 2000

A238359 Number of genus-9 rooted maps with n edges.

Original entry on oeis.org

11665426077721040625, 3498878057690404966500, 540996834819906946713375, 57494374008560749302297480, 4724172886681078698955547790, 320061005837218582787265273000, 18618409220753939214291224549409, 956146512935178711199035220384800, 44232688287025023758415781081779828, 1871678026675570344184400604255444240
Offset: 18

Views

Author

Joerg Arndt, Feb 26 2014

Keywords

Crossrefs

Column g=9 of A269919.
Rooted maps with n edges of genus g for 0 <= g <= 10: A000168, A006300, A006301, A104742, A215402, A238355, A238356, A238357, A238358, this sequence, A238360.

Programs

  • Mathematica
    T[0, 0] = 1; T[n_, g_] /; g < 0 || g > n/2 = 0; T[n_, g_] := T[n, g] = ((4 n - 2)/3 T[n - 1, g] + (2 n - 3) (2 n - 2) (2 n - 1)/12 T[n - 2, g - 1] + 1/2 Sum[(2 k - 1) (2 (n - k) - 1) T[k - 1, i] T[n - k - 1, g - i], {k, 1, n - 1}, {i, 0, g}])/((n + 1)/6);
    a[n_] := T[n, 9];
    Table[a[n], {n, 18, 30}] (* Jean-François Alcover, Jul 20 2018 *)
  • PARI
    \\ see A238396

A238360 Number of genus-10 rooted maps with n edges.

Original entry on oeis.org

15230046989184655753125, 5199629454143883380553750, 909887917857275652461097750, 108861830345440643086946970900, 10021124647635764856828690342402, 757187906770815991999545249667404, 48918614114003431712044170834572688, 2779227352989564224315657269511192976, 141720718575991991799057452443438430580
Offset: 20

Views

Author

Joerg Arndt, Feb 26 2014

Keywords

Crossrefs

Column g=10 of A269919.
Rooted maps with n edges of genus g for 0 <= g <= 10: A000168, A006300, A006301, A104742, A215402, A238355, A238356, A238357, A238358, A238359, this sequence.

Programs

  • Mathematica
    T[0, 0] = 1; T[n_, g_] /; g < 0 || g > n/2 = 0; T[n_, g_] := T[n, g] = ((4 n - 2)/3 T[n - 1, g] + (2 n - 3) (2 n - 2) (2 n - 1)/12 T[n - 2, g - 1] + 1/2 Sum[(2 k - 1) (2 (n - k) - 1) T[k - 1, i] T[n - k - 1, g - i], {k, 1, n - 1}, {i, 0, g}])/((n + 1)/6);
    a[n_] := T[n, 10];
    Table[a[n], {n, 20, 30}] (* Jean-François Alcover, Jul 20 2018 *)
  • PARI
    \\ see A238396

A085880 Triangle T(n,k) read by rows: multiply row n of Pascal's triangle (A007318) by the n-th Catalan number (A000108).

Original entry on oeis.org

1, 1, 1, 2, 4, 2, 5, 15, 15, 5, 14, 56, 84, 56, 14, 42, 210, 420, 420, 210, 42, 132, 792, 1980, 2640, 1980, 792, 132, 429, 3003, 9009, 15015, 15015, 9009, 3003, 429, 1430, 11440, 40040, 80080, 100100, 80080, 40040, 11440, 1430, 4862, 43758, 175032, 408408, 612612, 612612, 408408, 175032, 43758, 4862
Offset: 0

Views

Author

N. J. A. Sloane, Aug 17 2003

Keywords

Comments

Coefficients of terms in the series reversion of (1-k*x-(k+1)*x^2)/(1+x). - Paul Barry, May 21 2005
Equals A131427 * A007318 as infinite lower triangular matrices. [Philippe Deléham, Sep 15 2008]
Sum_{k=0..n} T(n,k)*x^k = A168491(n), A000007(n), A000108(n), A151374(n), A005159(n), A151403(n), A156058(n), A156128(n), A156266(n), A156270(n), A156273(n), A156275(n) for x = -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively. - Philippe Deléham, Nov 15 2013
Diagonal sums are A052709(n+1). - Philippe Deléham, Nov 15 2013

Examples

			Triangle starts:
[ 1]     1;
[ 2]     1,     1;
[ 3]     2,     4,      2;
[ 4]     5,    15,     15,      5;
[ 5]    14,    56,     84,     56,     14;
[ 6]    42,   210,    420,    420,    210,     42;
[ 7]   132,   792,   1980,   2640,   1980,    792,    132;
[ 8]   429,  3003,   9009,  15015,  15015,   9009,   3003,    429;
[ 9]  1430, 11440,  40040,  80080, 100100,  80080,  40040,  11440,  1430;
[10]  4862, 43758, 175032, 408408, 612612, 612612, 408408, 175032, 43758, 4862;
...
		

Programs

  • GAP
    Flat(List([0..10], n-> List([0..n], k-> Binomial(n,k)*Binomial(2*n,n)/( n+1) ))); # G. C. Greubel, Feb 07 2020
  • Magma
    [Binomial(n,k)*Catalan(n): k in [0..n], n in [0..10]]; // G. C. Greubel, Feb 07 2020
    
  • Maple
    seq(seq(binomial(n, k)*binomial(2*n, n)/(n+1), k = 0..n), n = 0..10); # G. C. Greubel, Feb 07 2020
  • Mathematica
    Table[Binomial[n, k]*CatalanNumber[n], {n,0,10}, {k,0,n}]//Flatten (* G. C. Greubel, Feb 07 2020 *)
  • PARI
    tabl(nn) = {for (n=0, nn, c =  binomial(2*n,n)/(n+1); for (k=0, n, print1(c*binomial(n, k), ", ");); print(););} \\ Michel Marcus, Apr 09 2015
    
  • Sage
    [[binomial(n,k)*catalan_number(n) for k in (0..n)] for n in (0..10)] # G. C. Greubel, Feb 07 2020
    

Formula

Triangle given by [1, 1, 1, 1, 1, 1, ...] DELTA [1, 1, 1, 1, 1, 1, ...] where DELTA is Deléham's operator defined in A084938.
Sum_{k>=0} T(n, k) = A151374(n) (row sums). - Philippe Deléham, Aug 11 2005
G.f.: (1-sqrt(1-4*(x+y)))/(2*(x+y)). - Vladimir Kruchinin, Apr 09 2015

A107264 Expansion of (1 - 3*x - sqrt((1-3*x)^2 - 4*3*x^2))/(2*3*x^2).

Original entry on oeis.org

1, 3, 12, 54, 261, 1323, 6939, 37341, 205011, 1143801, 6466230, 36960300, 213243435, 1240219269, 7263473148, 42799541886, 253556163243, 1509356586897, 9023497273548, 54154973176074, 326154592965879, 1970575690572297
Offset: 0

Views

Author

Paul Barry, May 15 2005

Keywords

Comments

Series reversion of x/(1+3x+3x^2). Transform of 3^n under the matrix A107131. A row of A107267.
Counts colored Motzkin paths, where H(1,0) and U(1,1) each have 3 colors and D(1,-1) one color. - Paul Barry, May 18 2005
Number of Motzkin paths of length n in which both the "up" and the "level" steps come in three colors. - Paul Barry, May 18 2005
Third binomial transform of 1,0,3,0,18,0,... or 3^n*C(n) (A005159) with interpolated zeros. - Paul Barry, May 24 2005
As a continued fraction, the g.f. is 1/(1-3*x-3*x^2/(1-3*x-3*x^2/(1-3*x-3*x^2/(1-3*x-3*x^2/(.... - Paul Barry, Dec 02 2008

Programs

  • Mathematica
    CoefficientList[Series[(1-3*x-Sqrt[1-6*x-3*x^2])/(6*x^2), {x, 0, 20}], x] (* Vaclav Kotesovec, Oct 17 2012 *)

Formula

G.f.: (1 - 3x - sqrt(1-6x-3x^2))/(6x^2);
a(n) = Sum_{k=0..n} (1/(k+1))*C(k+1, n-k+1)*C(n, k)3^k.
a(n) = Sum_{k=0..floor(n/2)} C(n, 2k)*C(k)*3^(n-k). - Paul Barry, May 18 2005
E.g.f.: exp(3x)*Bessel_I(1, sqrt(3)*2*x)/(sqrt(3)*x). - Paul Barry, May 24 2005
a(n) = (1/Pi)*Integral_{x=3-2*sqrt(3)..3+2*sqrt(3)} x^n*sqrt(-x^2 + 6*x + 3)/6. - Paul Barry, Sep 16 2006
a(n) = A156016(n+1)/3. - Philippe Deléham, Feb 04 2009
D-finite with recurrence: (n+2)*a(n) = 3*(2*n+1)*a(n-1) + 3*(n-1)*a(n-2). - Vaclav Kotesovec, Oct 17 2012
a(n) ~ (5+3*sqrt(3))*(3+2*sqrt(3))^n/(2*sqrt(2*Pi)*n^(3/2)). - Vaclav Kotesovec, Oct 17 2012
G.f.: Let F(x) be the g.f. of A348189 with offset 1, then F(x) = x + 2*x^2*F(x)^2*A(x*F(x)). - Alexander Burstein, Feb 14 2022

A156058 a(n) = 5^n * Catalan(n).

Original entry on oeis.org

1, 5, 50, 625, 8750, 131250, 2062500, 33515625, 558593750, 9496093750, 164023437500, 2870410156250, 50784179687500, 906860351562500, 16323486328125000, 295863189697265625, 5395152282714843750, 98911125183105468750, 1822047042846679687500
Offset: 0

Views

Author

Philippe Deléham, Feb 03 2009

Keywords

Comments

From Joerg Arndt, Oct 22 2012: (Start)
Number of strings of length 2*n of five different types of balanced parentheses.
The number of strings of length 2*n of t different types of balanced parentheses is given by t^n * A000108(n): there are n opening parentheses in the strings, giving t^n choices for the type (the closing parentheses are chosen to match). (End)
Number of Dyck paths of length 2n in which the step U=(1,1) come in 5 colors. [José Luis Ramírez Ramírez, Jan 31 2013]

Crossrefs

Programs

  • Magma
    [5^n*Catalan(n): n in [0..20]]; // Vincenzo Librandi, Jul 19 2011
  • Maple
    A156058_list := proc(n) local j, a, w; a := array(0..n); a[0] := 1;
    for w from 1 to n do a[w] := 5*(a[w-1]+add(a[j]*a[w-j-1],j=1..w-1)) od; convert(a,list)end: A156058_list(16); # Peter Luschny, May 19 2011
    A156058 := proc(n)
        5^n*A000108(n) ;
    end proc: # R. J. Mathar, Oct 06 2012
  • Mathematica
    Table[5^n CatalanNumber[n],{n,0,20}]  (* Harvey P. Dale, Mar 13 2011 *)

Formula

a(n) = 5^n*A000108(n).
From Gary W. Adamson, Jul 18 2011: (Start)
a(n) is the upper left term in M^n, M = the infinite square production matrix:
5, 5, 0, 0, 0, 0,...
5, 5, 5, 0, 0, 0,...
5, 5, 5, 5, 0, 0,...
5, 5, 5, 5, 5, 0,...
... (End)
E.g.f.: KummerM(1/2, 2, 20*x). - Peter Luschny, Aug 26 2012
D-finite with recurrence (n+1)*a(n) -10*(2*n-1)*a(n-1)=0. - R. J. Mathar, Oct 06 2012
G.f.: c(5*x) with c(x) the o.g.f. of A000108 (Catalan). - Philippe Deléham, Nov 15 2013
a(n) = Sum_{k=0..n} A085880(n,k)*4^k. - Philippe Deléham, Nov 15 2013
G.f.: 1/(1 - 5*x/(1 - 5*x/(1 - 5*x/(1 - ...)))), a continued fraction. - Ilya Gutkovskiy, Apr 19 2017
Sum_{n>=0} 1/a(n) = 410/361 + 600*arctan(1/sqrt(19)) / (361*sqrt(19)). - Vaclav Kotesovec, Nov 23 2021
Sum_{n>=0} (-1)^n/a(n) = 130/147 - 200*arctanh(1/sqrt(21)) / (147*sqrt(21)). - Amiram Eldar, Jan 25 2022

A151162 Number of walks within N^3 (the first octant of Z^3) starting at (0,0,0) and consisting of n steps taken from {(-1, 0, 0), (1, 0, 0), (1, 0, 1), (1, 1, 0)}.

Original entry on oeis.org

1, 3, 12, 45, 180, 702, 2808, 11097, 44388, 176418, 705672, 2812482, 11249928, 44903484, 179613936, 717517521, 2870070084, 11470898106, 45883592424, 183438670950, 733754683800, 2934026948196, 11736107792784, 46934017407594, 187736069630376, 750833732416212, 3003334929664848
Offset: 0

Views

Author

Manuel Kauers, Nov 18 2008

Keywords

Comments

Hankel transform is 3^C(n+1,2). - Philippe Deléham, Feb 01 2009
Inverse binomial transform of A151253. - Philippe Deléham, Feb 03 2009

Crossrefs

Programs

  • Magma
    [n le 3 select Factorial(n+1)/2 else (4*n*Self(n-1) + 12*(n-3)*Self(n-2) - 48*(n-3)*Self(n-3))/n: n in [1..41]]; // G. C. Greubel, Nov 09 2022
    
  • Mathematica
    aux[i_, j_, k_, n_]:= Which[Min[i, j, k, n]<0 || Max[i, j, k]>n, 0, n==0, KroneckerDelta[i, j, k, n], True, aux[i, j, k, n]= aux[-1+i, -1+j, k, -1+n] + aux[-1+i, j, -1+k, -1+n] + aux[-1+i, j, k, -1+n] + aux[1+i, j, k, -1+n]]; Table[Sum[aux[i,j,k,n], {i,0,n}, {j,0,n}, {k,0,n}], {n,0,30}]
    a[n_]:= a[n]= If[n<3, (n+2)!/2, (4*(n+1)*a[n-1] +12*(n-2)*a[n-2] -48*(n-2)*a[n- 3])/(n+1)]; Table[a[n], {n,0,40}] (* G. C. Greubel, Nov 09 2022 *)
  • SageMath
    def a(n): # a = A151162
        if (n==0): return 1
        elif (n%2==1): return 4*a(n-1) - 3^((n-1)/2)*catalan_number((n-1)/2)
        else: return 4*a(n-1)
    [a(n) for n in (0..40)] # G. C. Greubel, Nov 09 2022

Formula

a(n) = Sum_{k=0..n} A120730(n,k)*3^k. - Philippe Deléham, Feb 01 2009
From Philippe Deléham, Feb 02 2009: (Start)
a(2*n+2) = 4*a(2*n+1), a(2*n+1) = 4*a(2*n) - 3^n*A000108(n)
a(2*n+1) = 4*a(2*n) - A005159(n).
G.f.: (sqrt(1 - 12*x^2) + 6*x - 1)/(6*x*(1 - 4*x)). (End)
a(n) ~ 2^(2*n+1)/3. - Vaclav Kotesovec, Oct 31 2017
a(n) = (4*(n+1)*a(n-1) + 12*(n-2)*a(n-2) - 48*(n-2)*a(n-3))/(n+1). - G. C. Greubel, Nov 09 2022

A290605 Square array A(n,k), n>=0, k>=0, read by antidiagonals, where column k is the expansion of 2/(1 + sqrt(1 - 4*k*x)).

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 2, 2, 0, 1, 3, 8, 5, 0, 1, 4, 18, 40, 14, 0, 1, 5, 32, 135, 224, 42, 0, 1, 6, 50, 320, 1134, 1344, 132, 0, 1, 7, 72, 625, 3584, 10206, 8448, 429, 0, 1, 8, 98, 1080, 8750, 43008, 96228, 54912, 1430, 0, 1, 9, 128, 1715, 18144, 131250, 540672, 938223, 366080, 4862, 0
Offset: 0

Views

Author

Ilya Gutkovskiy, Aug 07 2017

Keywords

Comments

Number of 2n-length strings of balanced parentheses of at most k different types. Also number of binary trees with n inner nodes of at most k different dimensions. - Alois P. Heinz, Oct 28 2019

Examples

			G.f. of column k: A(x) = 1 + k*x + 2*k^2*x^2 + 5*k^3*x^3 + 14*k^4*x^4 + 42*k^5*x^5 + 132*k^6*x^6 + ...
Square array begins:
  1,   1,     1,      1,      1,       1,  ...
  0,   1,     2,      3,      4,       5,  ...
  0,   2,     8,     18,     32,      50,  ...
  0,   5,    40,    135,    320,     625,  ...
  0,  14,   224,   1134,   3584,    8750,  ...
  0,  42,  1344,  10206,  43008,  131250,  ...
		

Crossrefs

Rows n=0-2 give: A000012, A001477, A001105.
Main diagonal gives A291699.

Programs

  • Maple
    ctln:= proc(n) option remember; binomial(2*n, n)/(n+1) end:
    A:= proc(n, k) option remember; k^n*ctln(n) end:
    seq(seq(A(n, d-n), n=0..d), d=0..10);  # Alois P. Heinz, Oct 28 2019
  • Mathematica
    Table[Function[k, SeriesCoefficient[2/(1 + Sqrt[1 - 4 k x]), {x, 0, n}]][j - n], {j, 0, 10}, {n, 0, j}] // Flatten
    Table[Function[k, SeriesCoefficient[1/(1 + ContinuedFractionK[-k x, 1, {i, 1, n}]), {x, 0, n}]][j - n], {j, 0, 10}, {n, 0, j}] // Flatten

Formula

A(n,k) = k^n*(2*n)!/(n!*(n + 1)!).
A(n,k) = k^n*A000108(n).
G.f. of column k: 2/(1 + sqrt(1 - 4*k*x)).
G.f. of column k: 1/(1 - k*x/(1 - k*x/(1 - k*x/(1 - k*x/(1 - k*x/(1 - ...)))))), a continued fraction.
E.g.f. of column k: (BesselI(0,2*k*x) - BesselI(1,2*k*x))*exp(2*k*x).
If g.f. = 2/(1 + sqrt(1 - 4*k*x)), then a(n) ~ k^n*4^n/(sqrt(Pi)*n^(3/2)).
A(n,k) = Sum_{i=0..k} binomial(k,i) * A256061(n,k-i). - Alois P. Heinz, Oct 28 2019
For fixed k >= 1, Sum_{n>=0} 1/A(n,k) = 2*k*(8*k + 1) / (4*k - 1)^2 + 24 * k^2 * arcsin(1/(2*sqrt(k))) / (4*k - 1)^(5/2). - Vaclav Kotesovec, Nov 23 2021
For fixed k >= 1, Sum_{n>=0} (-1)^n / A(n,k) = 2*k*(8*k - 1) / (4*k + 1)^2 - 24 * k^2 * log((1 + sqrt(4*k + 1))/(2*sqrt(k))) / (4*k + 1)^(5/2). - Vaclav Kotesovec, Nov 24 2021
Showing 1-10 of 30 results. Next