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.

Previous Showing 21-25 of 25 results.

A026536 Irregular triangular array T read by rows: T(i,0 ) = T(i,2i) = 1 for i >= 0; T(i,1) = T(i,2i-1) = floor(i/2) for i >= 1; for even n >= 2, T(i,j) = T(i-1,j-2) + T(i-1,j-1) + T(i-1,j) for j = 2..2i-2, for odd n >= 3, T(i,j) = T(i-1,j-2) + T(i-1,j) for j = 2..2i-2.

Original entry on oeis.org

1, 1, 0, 1, 1, 1, 2, 1, 1, 1, 1, 3, 2, 3, 1, 1, 1, 2, 5, 6, 8, 6, 5, 2, 1, 1, 2, 6, 8, 13, 12, 13, 8, 6, 2, 1, 1, 3, 9, 16, 27, 33, 38, 33, 27, 16, 9, 3, 1, 1, 3, 10, 19, 36, 49, 65, 66, 65, 49, 36, 19, 10, 3, 1, 1, 4, 14, 32, 65, 104, 150, 180, 196, 180
Offset: 0

Views

Author

Keywords

Comments

T(n, k) is the number of strings s(0)..s(n) such that s(0) = 0, s(n) = n-k, |s(i) - s(i-1)| <= 1 if i is even, |s(i) - s(i-1)| = 1 if i is odd.

Examples

			First 5 rows:
  1
  1  0  1
  1  1  2  1  1
  1  1  3  2  3  1  1
  1  2  5  6  8  6  5  2  1
		

Crossrefs

Programs

  • Mathematica
    z = 12; t[n_, 0] := 1; t[n_, k_] := 1 /; k == 2 n; t[n_, 1] := Floor[n/2];
    t[n_, k_] := Floor[n/2] /; k == 2 n - 1; t[n_, k_] := t[n, k] =
    If[EvenQ[n], t[n - 1, k - 2] + t[n - 1, k - 1] + t[n - 1, k], t[n - 1, k -
    2] + t[n - 1, k]]; u = Table[t[n, k], {n, 0, z}, {k, 0, 2 n}];
    TableForm[u]   (* A026536 array *)
    v = Flatten[u] (* A026536 sequence *)
  • SageMath
    @cached_function
    def T(n, k):
        if k < 0 or n < 0: return 0
        elif k == 0 or k == 2*n: return 1
        elif k == 1 or k == 2*n-1: return n//2
        elif n % 2 == 1: return T(n-1, k-2) + T(n-1, k)
        return T(n-1, k-2) + T(n-1, k-1) + T(n-1, k) # Peter Luschny, Oct 13 2019

Extensions

Updated by Clark Kimberling, Aug 28 2014
Offset changed to 0 by Peter Luschny, Oct 10 2019

A026584 Irregular triangular array T read by rows: T(i,0) = T(i,2i) = 1 for i >= 0; T(i,1) = T(i,2i-1) = floor(i/2) for i >= 1; and for i >= 2 and j = 2..2i-2, T(i,j) = T(i-1,j-2) + T(i-1,j-1) + T(i-1,j) if i+j is odd, and T(i,j) = T(i-1,j-2) + T(i-1,j) if i+j is even.

Original entry on oeis.org

1, 1, 0, 1, 1, 1, 2, 1, 1, 1, 1, 4, 2, 4, 1, 1, 1, 2, 5, 7, 8, 7, 5, 2, 1, 1, 2, 8, 9, 20, 14, 20, 9, 8, 2, 1, 1, 3, 9, 19, 28, 43, 40, 43, 28, 19, 9, 3, 1, 1, 3, 13, 22, 56, 62, 111, 86, 111, 62, 56, 22, 13, 3, 1, 1, 4, 14, 38, 69, 140, 167, 259, 222, 259, 167, 140, 69, 38, 14, 4, 1
Offset: 1

Views

Author

Keywords

Comments

Row sums are in A026597. - Philippe Deléham, Oct 16 2006
T(n, k) = number of integer strings s(0)..s(n) such that s(0) = 0, s(n) = n-k, |s(i)-s(i-1)| <= 1 if s(i-1) odd, |s(i)-s(i-1)| = 1 if s(i-1) is even, for i = 1..n.

Examples

			First 5 rows:
  1
  1  0  1
  1  1  2  1  1
  1  1  4  2  4  1  1
  1  2  5  7  8  7  5  2  1
		

Crossrefs

Programs

  • Mathematica
    z = 12; t[n_, 0] := 1; t[n_, k_] := 1 /; k == 2 n; t[n_, 1] := Floor[n/2]; t[n_, k_] := Floor[n/2] /; k == 2 n - 1; t[n_, k_] := t[n, k] = If[EvenQ[n + k], t[n - 1, k - 2] + t[n - 1, k], t[n - 1, k - 2] + t[n - 1, k - 1] + t[n - 1, k]]; u = Table[t[n, k], {n, 0, z}, {k, 0, 2 n}];
    TableForm[u]   (* A026584 array *)
    v = Flatten[u] (* A026584 sequence *)
  • Sage
    @CachedFunction
    def T(n,k):
        if (k==0 or k==2*n): return 1
        elif (k==1 or k==2*n-1): return (n//2)
        else: return T(n-1, k-2) + T(n-1, k) if ((n+k)%2==0) else T(n-1, k-2) + T(n-1, k-1) + T(n-1, k)
    flatten([[T(n,k) for k in (0..2*n)] for n in (0..12)]) # G. C. Greubel, Dec 11 2021

Formula

T(n, k) = T(n-1, k-2) + T(n-1, k) if ( (n+k) mod 2 ) = 0, otherwise T(n-1, k-2) + T(n-1, k-1) + T(n-1, k), where T(n, 0) = T(n, 2*n) = 1, T(n, 1) = T(n, 2*n-1) = floor(n/2).

Extensions

Updated by Clark Kimberling, Aug 29 2014

A026568 Irregular triangular array T read by rows: T(i,0) = T(i,2i) = 1 for i >= 0; T(i,1) = T(i,2i-1) = [ (i+1)/2 ] for i >= 1; and for i >= 2 and 2 <=j <= i - 2, T(i,j) = T(i-1,j-2) + T(i-1,j-1) + T(i-1,j) if i + j is even, T(i,j) = T(i-1,j-2) + T(i-1,j) if i + j is odd.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 2, 4, 5, 4, 2, 1, 1, 2, 7, 7, 13, 7, 7, 2, 1, 1, 3, 8, 16, 20, 27, 20, 16, 8, 3, 1, 1, 3, 12, 19, 44, 43, 67, 43, 44, 19, 12, 3, 1, 1, 4, 13, 34, 56, 106, 111, 153, 111, 106, 56, 34, 13, 4, 1, 1, 4, 18, 38, 103, 140, 273
Offset: 1

Views

Author

Keywords

Comments

T(n, k) = number of strings s(0)..s(n) such that s(0) = 0, s(n) = n - k, |s(i)-s(i-1)| <= 1 if s(i-1) is even, |s(i)-s(i-1)| = 1 if s(i-1) is odd, for 1 <= i <= n.

Examples

			First 5 rows:
  1
  1  1  1
  1  1  3  1  1
  1  2  4  5  4  2  1
  1  2  7  7 13  7  7  2  1
		

Crossrefs

Cf. T(n,n) is A026569.

Programs

  • Mathematica
    z = 12; t[n_, 0] := 1; t[n_, 1] := Floor[(n + 1)/2]; t[n_, k_] := t[n, k] = Which[k == 2 n, 1, k == 2 n - 1, Floor[(n + 1)/2], EvenQ[n + k], t[n - 1, k - 2] + t[n - 1, k - 1] + t[n - 1, k], OddQ[n + k], t[n - 1, k - 2] + t[n - 1, k]]; u = Table[t[n, k], {n, 0, z}, {k, 0, 2 n}];
    TableForm[u] (* A026568 array *)
    Flatten[u]   (* A026568 sequence *)
  • PARI
    T(k,n)=if(n<0||n>2*k,0,if(n==0||n==2*k,1,if(k>0&&(n==1||n==2*k-1),(k+1)\2,T(k-1,n-2)+T(k-1,n)+if((k+n)%2==0,T(k-1,n-1))))) \\ Ralf Stephan

Extensions

Updated by Clark Kimberling, Aug 28 2014

A026532 Ratios of successive terms are 3, 2, 3, 2, 3, 2, 3, 2, ...

Original entry on oeis.org

1, 3, 6, 18, 36, 108, 216, 648, 1296, 3888, 7776, 23328, 46656, 139968, 279936, 839808, 1679616, 5038848, 10077696, 30233088, 60466176, 181398528, 362797056, 1088391168, 2176782336, 6530347008, 13060694016, 39182082048, 78364164096, 235092492288, 470184984576
Offset: 1

Views

Author

Keywords

Comments

Preface the series with a 1: (1, 1, 3, 6, 18, 36, ...); then the next term in the series = (1, 1, 3, 6, ...) dot (1, 2, 1, 2, ...). Example: 36 = (1, 1, 3, 6, 18) dot (1, 2, 1, 2, 1) = (1 + 2 + 3 + 12 + 18). - Gary W. Adamson, Apr 18 2009
Partial products of A176059. - Reinhard Zumkeller, Apr 04 2012

Crossrefs

Cf. A038730, A038792, and A134511 for incomplete Fibonacci sequences, and A324242 for incomplete Lucas sequences.

Programs

  • Haskell
    a026532 n = a026532_list !! (n-1)
    a026532_list = scanl (*) 1 $ a176059_list
    -- Reinhard Zumkeller, Apr 04 2012
    
  • Magma
    [(1/4)*(3-(-1)^n)*6^Floor(n/2) : n in [1..30]]; // Vincenzo Librandi, Jun 08 2011
    
  • Mathematica
    FoldList[(2 + Boole[EvenQ@ #2]) #1 &, Range@ 28] (* or *)
    CoefficientList[Series[x*(1+3x)/(1-6x^2), {x,0,31}], x] (* Michael De Vlieger, Aug 02 2017 *)
    LinearRecurrence[{0,6},{1,3},30] (* Harvey P. Dale, Jul 11 2018 *)
  • PARI
    a(n)=if(n%2,3,1)*6^(n\2) \\ Charles R Greathouse IV, Jul 02 2013
    
  • Python
    def a(n): return (3 if n%2 else 1)*6**(n//2)
    print([a(n) for n in range(31)]) # Indranil Ghosh, Aug 02 2017
    
  • Sage
    [(1/2)*6^((n-2)/2)*(3*(1+(-1)^n) + sqrt(6)*(1-(-1)^n)) for n in (1..30)] # G. C. Greubel, Dec 21 2021

Formula

a(n) = T(n, 0) + T(n, 1) + ... + T(n, 2n-2), T given by A026519.
From Benoit Cloitre, Nov 14 2003: (Start)
a(n) = (1/2)*(5+(-1)^n)*a(n-1) for n>1, a(1) = 1.
a(n) = (1/4)*(3-(-1)^n)*6^floor(n/2). (End)
From Ralf Stephan, Feb 03 2004: (Start)
G.f.: x*(1+3*x)/(1-6*x^2).
a(n+2) = 6*a(n). (End)
a(n+3) = a(n+2)*a(n+1)/a(n). - Reinhard Zumkeller, Mar 04 2011
a(n) = (1/2)*6^((n-2)/2)*(3*(1+(-1)^n) + sqrt(6)*(1-(-1)^n)). - G. C. Greubel, Dec 21 2021
Sum_{n>=1} 1/a(n) = 8/5. - Amiram Eldar, Feb 13 2023

Extensions

New definition from Ralf Stephan, Dec 01 2004
Offset changed from 0 to 1 by Vincenzo Librandi, Jun 08 2011

A026268 Triangle, T(n, k): T(n,k) = 1 for n < 3, T(3,1) = T(3,2) = T(3,3) = 2, T(n,0) = 1, T(n,1) = n-1, T(n,n) = T(n-1,n-2) + T(n-1,n-1), otherwise T(n,k) = T(n-1,k-2) + T(n-1,k-1) + T(n-1,k), read by rows.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 1, 3, 5, 6, 4, 1, 4, 9, 14, 15, 10, 1, 5, 14, 27, 38, 39, 25, 1, 6, 20, 46, 79, 104, 102, 64, 1, 7, 27, 72, 145, 229, 285, 270, 166, 1, 8, 35, 106, 244, 446, 659, 784, 721, 436, 1, 9, 44, 149, 385, 796, 1349, 1889, 2164, 1941, 1157, 1, 10, 54, 202, 578, 1330, 2530, 4034, 5402, 5994, 5262, 3098
Offset: 0

Views

Author

Keywords

Comments

a(n) = number of strings s(0)..s(n) such that s(n) = n-k, where s(0) = 0, s(1) = 1, |s(i)-s(i-1)| <= 1 for i >= 2; |s(2)-s(1)| = 1, and |s(3)-s(2)| = 1 if s(2) = 1.

Examples

			Triangle begins as:
  1;
  1, 1;
  1, 1,  1;
  1, 2,  2,   2;
  1, 3,  5,   6,   4;
  1, 4,  9,  14,  15,  10;
  1, 5, 14,  27,  38,  39,   25;
  1, 6, 20,  46,  79, 104,  102,   64;
  1, 7, 27,  72, 145, 229,  285,  270,  166;
  1, 8, 35, 106, 244, 446,  659,  784,  721,  436;
  1, 9, 44, 149, 385, 796, 1349, 1889, 2164, 1941, 1157;
		

Crossrefs

Programs

  • Magma
    f:= func< n | n eq 2 select 1 else (n^2 -n -2)/2 >;
    function T(n,k) // T = A026268
      if k eq 0 or n lt 3 then return 1;
      elif k eq 1 then return n-1;
      elif k eq 2 then return f(n);
      elif k eq n then return T(n-1, n-2) + T(n-1, n-1);
      else return T(n-1, k-2) + T(n-1, k-1) + T(n-1, k);
      end if; return T;
    end function;
    [T(n,k): k in [0..n], n in [0..14]]; // G. C. Greubel, Sep 24 2022
    
  • Mathematica
    T[n_, k_]:= T[n, k]= If[n<3 || k==0, 1, If[k==1, n-1, If[k==2, (n^2-n-2)/2 + Boole[n==2], If[k==n, T[n-1, n-2] +T[n-1, n-1], T[n-1, k-2] + T[n-1, k-1] + T[n -1, k] ]]]];
    Table[T[n, k], {n,0,14}, {k,0,n}]//Flatten (* corrected by G. C. Greubel, Sep 24 2022 *)
  • SageMath
    def T(n,k): # T = A026268
        if n<3 or k==0: return 1
        elif k==1: return n-1
        elif k==2: return (n^2 -n -2)//2 + int(n==2)
        elif k==n: return T(n-1, n-2) + T(n-1, n-1)
        else: return T(n-1, k-2) + T(n-1, k-1) + T(n-1, k)
    flatten([[T(n,k) for k in range(n+1)] for n in range(14)]) # G. C. Greubel, Sep 24 2022

Formula

From G. C. Greubel, Sep 24 2022: (Start)
T(n, 1) = A000027(n-1), n >= 1.
T(n, 2) = A212342(n-1), n >= 2.
T(n, n-1) = A026270(n), n >= 2.
T(n, n-2) = A026288(n), n >= 2.
T(n, n-3) = A026289(n), n >= 3.
T(n, n-4) = A026290(n), n >= 4.
T(n, n) = A026269(n), n >= 2.
T(n, floor(n/2)) = A026297(n), n >= 0.
T(2*n, n) = A026292(n).
T(2*n, n-1) = A026295(n), n >= 1.
T(2*n, n+1) = A026296(n), n >= 1.
T(2*n-1, n-1) = A026291(n), n >= 2.
T(3*n, n) = A026293(n), n >= 0.
T(4*n, n) = A026294(n), n >= 0.
Sum_{k=0..n} T(n, k) = A026299(n-1), n >= 3.(End)

Extensions

Updated by Clark Kimberling, Aug 29 2014
Indices of b-file corrected by Sidney Cadot, Jan 06 2023.
Previous Showing 21-25 of 25 results.