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.

A026907 Triangular array T read by rows (9-diamondization of Pascal's triangle). Step 1: t(n,k) = sum of 9 entries in diamond-shaped subarray of Pascal's triangle having vertices C(n,k), C(n+4,k+2), C(n+2,k), C(n+2,k+2). Step 2: T(n,k) = t(n,k) - t(0,0) + 1.

Original entry on oeis.org

1, 13, 13, 28, 44, 28, 46, 90, 90, 46, 67, 154, 198, 154, 67, 91, 239, 370, 370, 239, 91, 118, 348, 627, 758, 627, 348, 118, 148, 484, 993, 1403, 1403, 993, 484, 148, 181, 650, 1495, 2414, 2824, 2414, 1495, 650, 181, 217, 849, 2163, 3927, 5256, 5256, 3927, 2163, 849, 217
Offset: 0

Views

Author

Keywords

Examples

			Triangle starts:
   1;
  13,  13;
  28,  44,  28;
  46,  90,  90,  46;
  67, 154, 198, 154,  67;
  91, 239, 370, 370, 239,  91;
  ...
		

Crossrefs

Sums: A026915 (row), A026916, A026917, A026918 (diagonal).

Programs

  • Magma
    A026907:= func< n,k | Binomial(n,k) + 3*Binomial(n+4,k+2) - 18 >;
    [A026907(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Aug 21 2025
    
  • Mathematica
    t[n_, k_]:=Binomial[n + 4, k + 2 ] + Binomial[n + 3, k + 1] + Binomial[n + 3, k + 2] + Binomial[n + 2, k] + Binomial[n + 2, k + 1] + Binomial[n + 2, k + 2] + Binomial[n + 1, k] + Binomial[n + 1, k + 1] + Binomial[n, k] ; T[n_, k_]:=t[n,k] - t[0, 0] + 1; Flatten[Table[T[n, k], {n, 0, 9},{k, 0, n}]] (* Indranil Ghosh, Mar 13 2017 *)
  • PARI
    alias(C, binomial);
    t(n,k) = C(n+4,k+2) + C(n+3,k+1) + C(n+3,k+2) + C(n+2,k) + C(n+2,k+1) + C(n+2,k+2) + C(n+1,k) + C(n+1,k+1) + C(n,k);
    T(n,k) = t(n,k)-t(0,0)+1;
    tabl(nn) = for (n=0, nn, for (k=0, n, print1(T(n,k), ", ")); print());
    \\ Michel Marcus, Mar 13 2017
    
  • SageMath
    def A026907(n,k): return binomial(n,k) +3*binomial(n+4,k+2) -18
    print(flatten([[A026907(n,k) for k in range(n+1)] for n in range(13)])) # G. C. Greubel, Aug 21 2025

Formula

From G. C. Greubel, Aug 21 2025: (Start)
T(n, k) = binomial(n,k) + 3*binomial(n+4, k+2) - 18.
Sum_{k=0..n} (-1)^k*T(n, k) = 6*(1+(-1)^n)* floor((n+1)/2) + [n=0]. (End)

A026909 a(n) = (1/2)*A026907(2*n, n).

Original entry on oeis.org

22, 99, 379, 1412, 5265, 19758, 74637, 283560, 1082449, 4148603, 15953607, 61526969, 237876571, 921678876, 3577968081, 13913243136, 54184698801, 211307360871, 825059443551, 3225071709981, 12619275028611, 49423455006501, 193732625020419, 760001601263697, 2983614441691035, 11720995167614703, 46074369462135607
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    A026909:= func< n | ((n+1)*Catalan(n) +3*(n+3)*Catalan(n+2))/2 -9 >;
    [A026909(n): n in [1..40]]; // G. C. Greubel, Aug 22 2025
    
  • Mathematica
    With[{b=Binomial}, Table[(b[2*n,n] +3*b[2*n+4,n+2] -18)/2, {n,40}]] (* G. C. Greubel, Aug 22 2025 *)
  • SageMath
    def A026909(n): return (binomial(2*n,n) +3*binomial(2*n+4,n+2))//2 -9
    print([A026909(n) for n in range(1,41)]) # G. C. Greubel, Aug 22 2025

Formula

From G. C. Greubel, Aug 22 2025: (Start)
a(n) = (binomial(2*n, n) + 3*binomial(2*n+4, n+2) - 18)/2.
a(n) = ((n+2)*A000108(n) + 3*(n+3)*A000108(n+2))/2 - 9.
G.f.: 2*(3 + 25*x - x^2 - (3 - 13*x + x^2)*sqrt(1-4*x))/((1-x)*sqrt(1-4*x)*(1 + sqrt(1 - 4*x))^2).
E.g.f.: (1/(2*x))*exp(2*x)*( 25*x*BesselI(0, 2*x) - 6*(1-4*x)*BesselI(1, 2*x) ) - 9*exp(x) - 1/2. (End)

Extensions

More terms added by G. C. Greubel, Aug 22 2025
Showing 1-2 of 2 results.