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 13 results. Next

A064036 Number of walks of length n on cubic lattice, starting at origin, staying in first (nonnegative) octant.

Original entry on oeis.org

1, 3, 12, 51, 234, 1110, 5460, 27405, 140490, 729918, 3845016, 20447658, 109801692, 593806356, 3234529584, 17715445605, 97567971930, 539701180590, 2998595422680, 16719506691030, 93559970043540, 525093580540620, 2955822168597480, 16680150247605390, 94365481922990460
Offset: 0

Views

Author

Henry Bottomley, Aug 23 2001

Keywords

Examples

			a(2)=12 since a(1) is obviously 3 and from each of these three positions there are four possible steps which remain in the first octant.

Links

Crossrefs

Cf. A064037. The two- and one-dimensional equivalents are A005566 and A001405. With no restriction on the walks, the number is 6^n, i.e. A000400.

Programs

  • Maple
    S:= series((BesselI(0,2*x)+BesselI(1,2*x))^3, x, 101):
seq(simplify(coeff(S,x,n))*n!, n=0..100); # Robert Israel, Oct 10 2016

Formula

a(n) = sum_j[C(n, j)B(j)B(j+1)B(n-j)] where B(k)=C(k, [k/2])=A001405(k)
E.g.f.: (BesselI(0, 2*x)+BesselI(1, 2*x))^3. - Vladeta Jovovic, Apr 28 2003
From Vaclav Kotesovec, Jun 10 2019: (Start)
Recurrence: (n+1)*(n+2)*(n+3)*a(n) = 4*(5*n^2+10*n+3)*a(n-1) + 4*(n-1)*(10*n^2+10*n-9)*a(n-2) - 144*(n-2)*(n-1)*a(n-3) - 144*(n-3)*(n-2)*(n-1)*a(n-4).
a(n) ~ 6^(n + 3/2) / (Pi*n)^(3/2). (End)

A064045 Square array read by antidiagonals of number of length 2k walks on an n-dimensional hypercubic lattice starting and finishing at the origin and staying in the nonnegative part.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 2, 2, 1, 0, 5, 10, 3, 1, 0, 14, 70, 24, 4, 1, 0, 42, 588, 285, 44, 5, 1, 0, 132, 5544, 4242, 740, 70, 6, 1, 0, 429, 56628, 73206, 16016, 1525, 102, 7, 1, 0, 1430, 613470, 1403028, 410928, 43470, 2730, 140, 8, 1, 0, 4862, 6952660, 29082339, 11925672, 1491210, 96684, 4445, 184, 9, 1
Offset: 0

Views

Author

Henry Bottomley, Aug 23 2001

Keywords

Examples

			Rows start:
1, 0,  0,   0,    0,     0,       0, ...
1, 1,  2,   5,   14,    42,     132, ...
1, 2, 10,  70,  588,  5544,   56628, ...
1, 3, 24, 285, 4242, 73206, 1403028, ...

Links

Crossrefs

Rows include A000007, A000108, A005568, A064037. Columns include A000012, A001477, A049450, A064046. Cf. A064044.

Programs

  • Maple
    a:= proc(n, k) option remember; `if`(n=0, `if`(k=0, 1, 0),
       add(binomial(2*k, 2*j)*binomial(2*j, j)/
       (j+1)*a(n-1, k-j), j=0..k))
    end:
seq(seq(a(n, d-n), n=0..d), d=0..12);  # Alois P. Heinz, May 06 2014
  • Mathematica
    a[n_, k_] := a[n, k] = If[n == 0, If[k == 0, 1, 0], Sum[Binomial[2*k, 2*j]* Binomial[2*j, j]/(j+1)*a[n-1, k-j], {j, 0, k}]]; Table[Table[a[n, d-n], {n, 0, d}], {d, 0, 12}] // Flatten (* Jean-François Alcover, Feb 26 2015, after Alois P. Heinz *)

Formula

a(n,k) = Sum_{j=0..k} C(2k,2j) c(j) a(n-1,k-j) where c(j) = C(2j,j)/(j+1) = A000108(j) with a(0,0) = 1 and a(0,k) = 0 for k>0.

A302181 Number of 3D walks of type abb.

Original entry on oeis.org

1, 5, 62, 1065, 21714, 492366, 12004740, 308559537, 8255788970, 227976044010, 6457854821340, 186814834574550, 5500292590186380, 164387681345290500, 4976887208815547640, 152378485941172462785, 4711642301137121933850, 146964278352052950118770, 4619875954522866283392300
Offset: 0

Views

Author

N. J. A. Sloane, Apr 09 2018

Keywords

Comments

See Dershowitz (2017) for precise definition.

Links

  • Nachum Dershowitz, Touchard's Drunkard, Journal of Integer Sequences, Vol. 20 (2017), #17.1.5.

Crossrefs

Cf. A000108, A000984, A002212, A002896, A005572, A026375, A064037, A081671, A138547, A145847, A145867, A150500, A202814.
Cf. A005558, A126869.

Programs

  • Maple
    C := n-> binomial(2*n, n)/(n+1): # Catalan numbers
A302181 := n-> add(binomial(2*n, k)*C(iquo(k+1, 2))*C(iquo(k, 2))*(2*iquo(k, 2)+1)*add((-1)^(k+j)*binomial(2*n-k, iquo(j,2)), j=0..2*n-k), k=0..2*n): seq(A302181(n), n = 0 .. 18); # Mélika Tebni, Nov 06 2024

Formula

a(n) = Sum_{k=0..2*n} binomial(2*n, k) * A005558(k) * A126869(2*n-k). - Mélika Tebni, Nov 06 2024

Extensions

a(8)-a(18) from Nachum Dershowitz, Aug 03 2020

A302180 Number of 3D walks of type aad.

Original entry on oeis.org

1, 1, 3, 7, 23, 71, 251, 883, 3305, 12505, 48895, 193755, 783355, 3205931, 13302329, 55764413, 236174933, 1008773269, 4343533967, 18834033443, 82201462251, 360883031291, 1592993944723, 7066748314147, 31493800133173, 140953938878821, 633354801073571, 2856369029213263
Offset: 0

Views

Author

N. J. A. Sloane, Apr 09 2018

Keywords

Comments

See Dershowitz (2017) for precise definition.
Number of 3D walks of length n in the first octant using steps (1, 1, 0), (1, -1, 0), (1, 0, 1), (1, 0, -1) and (1, 0, 0) that start at the origin and end at (n, 0, 0). The analogous problem in 2D is given by the Motzkin numbers A001006. - Farzan Byramji, Mar 06 2021
Inverse binomial transform of A145867 (Number of 3D walks of type aae). - Mélika Tebni, Nov 05 2024

Links

  • Nachum Dershowitz, Touchard's Drunkard, Journal of Integer Sequences, Vol. 20 (2017), #17.1.5.

Crossrefs

Cf. A000108, A000984, A001006, A002212, A002896, A005572, A026375, A064037, A081671, A138547, A145847, A145867, A150500, A202814.

Programs

  • Maple
    M := n-> add(binomial(n, 2*k)*binomial(2*k, k)/(k+1), k = 0 .. iquo(n,2)): # Motzkin numbers
A302180 := n-> add((-1)^(n-k)*binomial(n, k)*add(binomial(k, j)*M(j)*M(k-j), j=0..k), k=0..n):  seq(A302180(n), n = 0 .. 26); # Mélika Tebni, Nov 05 2024

Extensions

a(14)-a(26) from Farzan Byramji, Mar 06 2021

A302182 Number of 3D walks of type abc.

Original entry on oeis.org

1, 1, 5, 12, 62, 200, 1065, 3990, 21714, 89082, 492366, 2147376, 12004740, 54718092, 308559537, 1454116950, 8255788970, 39935276810, 227976044010, 1126178350440, 6457854821340, 32456552441040, 186814834574550, 952569927106980, 5500292590186380, 28391993275117500
Offset: 0

Views

Author

N. J. A. Sloane, Apr 09 2018

Keywords

Comments

See Dershowitz (2017) for precise definition.

Links

  • Nachum Dershowitz, Touchard's Drunkard, Journal of Integer Sequences, Vol. 20 (2017), #17.1.5.

Crossrefs

Cf. A000108, A000984, A002212, A002896, A005572, A026375, A064037, A081671, A138547, A145847, A145867, A150500, A202814.
Cf. A018224, A126120, A135394, A302181.

Programs

  • Python
    from math import comb as binomial
def row(n: int) -> list[int]:
    return sum(binomial(n, k)*binomial(k, k//2)//(k//2+1)*((k+1) %2)*binomial(n-k, (n-k)//2)**2 for k in range(n+1))
for n in range(26): print(row(n)) # Mélika Tebni, Nov 27 2024

Formula

From Mélika Tebni, Nov 27 2024: (Start)
a(n) = Sum_{k=0..n} binomial(n, k)*A126120(k)*A018224(n-k).
a(2*n+1) = A135394(n) / (2*n+2).
a(2*n) = A302181(n). (End)

Extensions

a(13)-a(25) from Mélika Tebni, Nov 27 2024

A302184 Number of 3D walks of type abe.

Original entry on oeis.org

1, 2, 7, 26, 108, 472, 2159, 10194, 49396, 244328, 1229308, 6273896, 32410096, 169181664, 891181607, 4731912082, 25302648644, 136150941064, 736747902236, 4007011320808, 21893702201648, 120125750018656, 661630546993116, 3656966382542984, 20278320788680912, 112782556853239712
Offset: 0

Views

Author

N. J. A. Sloane, Apr 09 2018

Keywords

Comments

See Dershowitz (2017) for precise definition.

Links

  • Nachum Dershowitz, Touchard's Drunkard, Journal of Integer Sequences, Vol. 20 (2017), #17.1.5.

Crossrefs

Cf. A000108, A000984, A002212, A002896, A005572, A026375, A064037, A081671, A126120, A138547, A145847, A145867, A150500, A202814.

Programs

  • Maple
    a := n -> 2*add(binomial(n, k)*binomial(k, k/2)*binomial(2*(n-k), n-k)/(k+2), k = 0..n, 2): seq(a(n), n = 0..25);  # Peter Luschny, Nov 30 2024
  • Python
    from math import comb as binomial
def a(n: int):
    return sum(binomial(n, k)*binomial(k, k//2)//(k//2+1)*((k+1) %2)*binomial(2*(n-k), n-k) for k in range(n+1))
print([a(n) for n in range(26)]) # Mélika Tebni, Nov 30 2024

Formula

a(n) = Sum_{k=0..n} binomial(n, k)*A126120(k)*A000984(n-k). - Mélika Tebni, Nov 30 2024

Extensions

a(12)-a(25) from Mélika Tebni, Nov 30 2024

A302178 The number of 3D walks of semilength n in a quadrant returning to the origin.

Original entry on oeis.org

1, 4, 40, 570, 9898, 195216, 4209084, 96941130, 2349133930, 59272544760, 1545550116240, 41416083787260, 1135679731004700, 31760915181412800, 903492759037272480, 26086451983000501410, 763124703525758894490, 22585374873810849150600, 675419388009799152812400
Offset: 0

Views

Author

N. J. A. Sloane, Apr 09 2018

Keywords

Links

  • Nachum Dershowitz, Touchard's Drunkard, Journal of Integer Sequences, Vol. 20 (2017), #17.1.5. The sequence is type aab in Table 3.

Crossrefs

Cf. A000108, A000984, A002212, A002896, A005572, A026375, A064037, A081671, A138547, A145847, A145867, A150500, A202814.

Formula

a(n) = Sum_{i=0..n,j=0..n-i} A000108(i) * A000108(j) * A000984_(n-i-j) * (2n)!/((2i)!*(2j)!*(2n-2i-2j)!). - Nachum Dershowitz, Aug 13 2020

Extensions

a(8)-a(18) from Nachum Dershowitz, Aug 03 2020
Name edited by Nachum Dershowitz, Aug 13 2020

A302179 The number of 3D walks of length n in an octant returning to axis of origin.

Original entry on oeis.org

1, 1, 4, 9, 40, 120, 570, 1995, 9898, 38178, 195216, 805266, 4209084, 18239364, 96941130, 436235085, 2349133930, 10891439130, 59272544760, 281544587610, 1545550116240, 7489973640240, 41416083787260, 204122127237210, 1135679731004700, 5678398655023500, 31760915181412800, 160789633105902300
Offset: 0

Views

Author

N. J. A. Sloane, Apr 09 2018

Keywords

Links

Crossrefs

Cf. A000108, A000984, A002212, A002896, A005572, A026375, A064037, A081671, A138547, A145847, A145867, A150500, A202814.

Programs

  • PARI
    C(n) = binomial(2*n, n)/(n+1); \\ A000108
f(n) = binomial(n, floor(n/2)); \\ A001405
a(n) = sum(i=0, n, if (!(i%2), sum(j=0, n-i, if (!(j%2), C(i/2)*C(j/2)*f(n-i-j)*n!/(i! * j! * (n-i-j)!))))); \\ Michel Marcus, Aug 07 2020

Formula

a(n) = Sum_{i=0..n, j=0..n-i, i,j even} A126120(i) * A126120(j) * A001405(n-i-j) * n!/(i! * j! * (n-i-j)!). - Nachum Dershowitz, Aug 06 2020
E.g.f.: (BesselI(1, 2*x)/x)^2*(BesselI(0, 2*x) + BesselI(1, 2*x)). - Mélika Tebni, Jan 06 2025

Extensions

a(13)-a(27) from Nachum Dershowitz, Aug 04 2020

A302183 Number of 3D n-step walks of type abd.

Original entry on oeis.org

1, 1, 4, 10, 39, 131, 521, 1989, 8149, 33205, 139870, 592120, 2552155, 11079303, 48639722, 214997228, 957817013, 4292316197, 19349957108, 87663905954, 399038606291, 1823961268751, 8369603968599, 38540835938335, 178056111047329, 825079806039121, 3833960405339446
Offset: 0

Views

Author

N. J. A. Sloane, Apr 09 2018

Keywords

Comments

See Dershowitz (2017) for precise definition.

Links

  • Nachum Dershowitz, Touchard's Drunkard, Journal of Integer Sequences, Vol. 20 (2017), #17.1.5.

Crossrefs

Cf. A000108, A000984, A002212, A002896, A005572, A026375, A064037, A081671, A138547, A145847, A145867, A150500, A202814.
Cf. A001006, A126869, A302184.

Programs

  • Python
    from math import comb as binomial
def M(n): return sum(binomial(n, 2*k)*binomial(2*k, k)//(k+1) for k in range(n//2+1)) # Motzkin numbers
def a(n):
    return sum(binomial(n, k)*binomial(k, k//2)*((k+1) %2)*M(n-k) for k in range(n+1))
print([a(n) for n in range(27)]) # Mélika Tebni, Dec 03 2024

Formula

From Mélika Tebni, Dec 03 2024: (Start)
a(n) = Sum_{k=0..n} binomial(n, k)*A126869(k)*A001006(n-k).
Inverse binomial transform of A302184. (End)

Extensions

a(13)-a(26) from Mélika Tebni, Dec 03 2024

A302185 Number of 3D n-step walks of type acc.

Original entry on oeis.org

1, 2, 7, 24, 98, 400, 1785, 7980, 37674, 178164, 874146, 4294752, 21667932, 109436184, 563910633, 2908233900, 15235550330, 79870553620, 424021948350, 2252356700880, 12088746573540, 64913104882080, 351594254659830, 1905139854213960, 10399223643879420, 56783986550235000
Offset: 0

Views

Author

N. J. A. Sloane, Apr 09 2018

Keywords

Comments

See Dershowitz (2017) for precise definition.

Links

  • Nachum Dershowitz, Touchard's Drunkard, Journal of Integer Sequences, Vol. 20 (2017), #17.1.5.

Crossrefs

Cf. A000108, A000984, A002212, A002896, A005572, A026375, A064037, A081671, A138547, A145847, A145867, A150500, A202814.
Cf. A000918, A001405, A005558, A005566, A135394, A302182.

Programs

  • Maple
    b:= n-> binomial(n, floor(n/2))*binomial(n+1, floor((n+1)/2)):
C:= n-> binomial(2*n, n)/(n+1):
a:= n-> add(binomial(n, 2*k)*C(k)*b(n-2*k), k=0..n/2):
seq(a(n), n=0..25);  # Alois P. Heinz, Dec 06 2024
# second Maple program:
a:= proc(n) option remember; `if`(n<4, [1, 2, 7, 24][n+1],
      (8*(14*n^4+85*n^3+190*n^2+188*n+63)*a(n-1)+4*(n-1)*
      (80*n^4+418*n^3+676*n^2+269*n-108)*a(n-2)-96*(n-1)*(n-2)*
      (10*n^2+31*n+27)*a(n-3)-144*(n-1)*(n-2)*(n-3)*(8*n^2+33*n+36)*
       a(n-4))/((n+4)*(n+3)*(n+2)*(8*n^2+17*n+11)))
    end:
seq(a(n), n=0..25);  # Alois P. Heinz, Dec 06 2024
  • Mathematica
    b[n_] := Binomial[n, Floor[n/2]]*Binomial[n+1, Floor[(n+1)/2]];
c[n_] := Binomial[2*n, n]/(n+1);
a[n_] := Sum[Binomial[n, 2*k]*c[k]*b[n - 2*k], {k, 0, n/2}];
Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Jan 28 2025, after Alois P. Heinz *)
  • Python
    from math import comb as binomial
def C(n): return (binomial(2*n, n)//(n+1)) # Catalan numbers
def a(n):
    return sum(binomial(n, k)*C((k+1)//2)*C(k//2)*(2*(k//2)+1)*binomial(n-k, (n-k)//2) for k in range(n+1))
print([a(n) for n in range(26)]) # Mélika Tebni, Dec 06 2024

Formula

From Mélika Tebni, Dec 06 2024: (Start)
E.g.f.: (BesselI(0, 2*x) + BesselI(1, 2*x))^2*BesselI(1, 2*x) / x.
a(n) = Sum_{k=0..n} binomial(n, k)*A005558(k)*A001405(n-k).
a(2*n+1) = 2*A302182(2*n+1) = A135394(n) / (n+1).
For n > 0, a(A000918(n)) is odd. (End)

Extensions

a(13)-a(25) from Mélika Tebni, Dec 06 2024
Showing 1-10 of 13 results. Next

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.