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 11-20 of 25 results. Next

A222134 Decimal expansion of sqrt(5 + sqrt(5 + sqrt(5 + sqrt(5 + ... )))).

Original entry on oeis.org

2, 7, 9, 1, 2, 8, 7, 8, 4, 7, 4, 7, 7, 9, 2, 0, 0, 0, 3, 2, 9, 4, 0, 2, 3, 5, 9, 6, 8, 6, 4, 0, 0, 4, 2, 4, 4, 4, 9, 2, 2, 2, 8, 2, 8, 8, 3, 8, 3, 9, 8, 5, 9, 5, 1, 3, 0, 3, 6, 2, 1, 0, 6, 1, 9, 5, 3, 4, 3, 4, 2, 1, 2, 7, 7, 3, 8, 8, 5, 4, 4, 3, 3, 0, 2, 1, 8, 0, 7, 7, 9, 7, 4, 6, 7, 2, 2, 5, 1, 6, 3
Offset: 1

Views

Author

Jaroslav Krizek, Feb 08 2013

Keywords

Comments

c^n = A015440(n) + A015440(n-1) * A222135, where c = (1 + sqrt(21))/2 and A222135 = (-1 + sqrt(21))/2. - Gary W. Adamson, Nov 27 2023

Examples

			2.791287847477920003294023596864...
		

Crossrefs

Programs

Formula

Equals (sqrt(21) + 1)/2 = A090458 - 1 = A107905 - 2 = A222135 + 1.
sqrt(5 + sqrt(5 + sqrt(5 + sqrt(5 + ... )))) - 1 = sqrt(5 - sqrt(5 - sqrt(5 - sqrt(5 - ... )))) = A222135.
Minimal polynomial: x^2 - x - 5. - Stefano Spezia, Jul 02 2025

A180250 a(n) = 5*a(n-1) + 10*a(n-2), with a(1)=0 and a(2)=1.

Original entry on oeis.org

0, 1, 5, 35, 225, 1475, 9625, 62875, 410625, 2681875, 17515625, 114396875, 747140625, 4879671875, 31869765625, 208145546875, 1359425390625, 8878582421875, 57987166015625, 378721654296875, 2473479931640625, 16154616201171875, 105507880322265625
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    [n le 2 select n-1 else 5*Self(n-1) +10*Self(n-2): n in [1..30]]; // G. C. Greubel, Jan 16 2018
    
  • Mathematica
    Join[{a=0,b=1},Table[c=5*b+10*a;a=b;b=c,{n,100}]]
    LinearRecurrence[{5,10}, {0,1}, 30] (* G. C. Greubel, Jan 16 2018 *)
  • PARI
    a(n)=([0,1;10,5]^(n-1))[1,2] \\ Charles R Greathouse IV, Oct 03 2016
    
  • PARI
    my(x='x+O('x^30)); concat([0], Vec(x^2/(1-5*x-10*x^2))) \\ G. C. Greubel, Jan 16 2018
    
  • SageMath
    A180250= BinaryRecurrenceSequence(5,10,0,1)
    [A180250(n-1) for n in range(1,41)] # G. C. Greubel, Jul 21 2023

Formula

a(n) = ((5+sqrt(65))^(n-1) - (5-sqrt(65))^(n-1))/(2^(n-1)*sqrt(65)). - Rolf Pleisch, May 14 2011
G.f.: x^2/(1-5*x-10*x^2).
a(n) = (i*sqrt(10))^(n-1) * ChebyshevU(n-1, -i*sqrt(5/8)). - G. C. Greubel, Jul 21 2023

A015551 Expansion of x/(1 - 6*x - 5*x^2).

Original entry on oeis.org

0, 1, 6, 41, 276, 1861, 12546, 84581, 570216, 3844201, 25916286, 174718721, 1177893756, 7940956141, 53535205626, 360916014461, 2433172114896, 16403612761681, 110587537144566, 745543286675801, 5026197405777636
Offset: 0

Views

Author

Keywords

Comments

Let the generator matrix for the ternary Golay G_12 code be [I|B], where the elements of B are taken from the set {0,1,2}. Then a(n)=(B^n)1,2 for instance. - _Paul Barry, Feb 13 2004
Pisano period lengths: 1, 2, 4, 4, 1, 4, 42, 8, 12, 2, 10, 4, 12, 42, 4, 16, 96, 12, 360, 4, ... - R. J. Mathar, Aug 10 2012

Crossrefs

Programs

  • Magma
    I:=[0,1]; [n le 2 select I[n] else 6*Self(n-1)+5*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Nov 14 2011
    
  • Mathematica
    Join[{a=0,b=1},Table[c=6*b+5*a;a=b;b=c,{n,100}]] (* Vladimir Joseph Stephan Orlovsky, Jan 16 2011 *)
    CoefficientList[Series[x/(1-6x-5x^2),{x,0,20}],x] (* or *) LinearRecurrence[ {6,5},{0,1},30] (* Harvey P. Dale, Oct 30 2017 *)
  • PARI
    a(n)=([0,1; 5,6]^n*[0;1])[1,1] \\ Charles R Greathouse IV, Oct 03 2016
  • Sage
    [lucas_number1(n,6,-5) for n in range(0, 21)] # Zerinvary Lajos, Apr 24 2009
    

Formula

a(n) = 6*a(n-1) + 5*a(n-2).
a(n) = sqrt(14)*(3+sqrt(14))^n/28 - sqrt(14)*(3-sqrt(14))^n/28. - Paul Barry, Feb 13 2004

A185540 T(n,k)=1/6 the number of nXk 0..5 arrays with every element equal to exactly one or two of its horizontal and vertical neighbors.

Original entry on oeis.org

0, 1, 1, 1, 11, 1, 6, 110, 110, 6, 11, 1155, 2990, 1155, 11, 41, 12105, 127570, 127570, 12105, 41, 96, 127190, 4403225, 17526450, 4403225, 127190, 96, 301, 1336205, 167398110, 2241330570, 2241330570, 167398110, 1336205, 301, 781, 14039455
Offset: 1

Views

Author

R. H. Hardin Jan 30 2011

Keywords

Comments

Table starts
....0.........1............1..............6..............11..............41
....1........11..........110...........1155...........12105..........127190
....1.......110.........2990.........127570.........4403225.......167398110
....6......1155.......127570.......17526450......2241330570....295044602325
...11.....12105......4403225.....2241330570....994633804270.468362874353295
...41....127190....167398110...295044602325.468362874353295
...96...1336205...6079401960.38418831900875
..301..14039455.225574291125
..781.147510990
.2286

Examples

			Some solutions for 4X3 with a(1,1)=0
..0..1..1....0..0..1....0..0..2....0..0..1....0..0..1....0..4..4....0..0..0
..0..0..5....5..5..1....4..0..2....3..4..1....0..4..1....0..3..3....0..1..0
..1..5..5....0..2..4....4..4..3....3..4..1....4..4..0....1..0..0....0..1..0
..1..4..4....0..2..4....1..1..3....3..4..4....2..2..0....1..3..3....0..5..5
		

Crossrefs

Column 1 is A015440(n-2)

A193376 T(n,k) = number of ways to place any number of 2 X 1 tiles of k distinguishable colors into an n X 1 grid; array read by descending antidiagonals, with n, k >= 1.

Original entry on oeis.org

1, 1, 2, 1, 3, 3, 1, 4, 5, 5, 1, 5, 7, 11, 8, 1, 6, 9, 19, 21, 13, 1, 7, 11, 29, 40, 43, 21, 1, 8, 13, 41, 65, 97, 85, 34, 1, 9, 15, 55, 96, 181, 217, 171, 55, 1, 10, 17, 71, 133, 301, 441, 508, 341, 89, 1, 11, 19, 89, 176, 463, 781, 1165, 1159, 683, 144, 1, 12, 21, 109, 225, 673
Offset: 1

Views

Author

R. H. Hardin, Jul 24 2011

Keywords

Comments

Transposed variant of A083856. - R. J. Mathar, Aug 23 2011
As to the sequences by columns beginning (1, N, ...), let m = (N-1). The g.f. for the sequence (1, N, ...) is 1/(1 - x - m*x^2). Alternatively, the corresponding matrix generator is [[1,1], [m,0]]. Another equivalency is simply: The sequence beginning (1, N, ...) is the INVERT transform of (1, m, 0, 0, 0, ...). Convergents to the sequences a(n)/a(n-1) are (1 + sqrt(4*m+1))/2. - Gary W. Adamson, Feb 25 2014

Examples

			Array T(n,k) (with rows n >= 1 and column k >= 1) begins as follows:
  ..1...1....1....1.....1.....1.....1......1......1......1......1......1...
  ..2...3....4....5.....6.....7.....8......9.....10.....11.....12.....13...
  ..3...5....7....9....11....13....15.....17.....19.....21.....23.....25...
  ..5..11...19...29....41....55....71.....89....109....131....155....181...
  ..8..21...40...65....96...133...176....225....280....341....408....481...
  .13..43...97..181...301...463...673....937...1261...1651...2113...2653...
  .21..85..217..441...781..1261..1905...2737...3781...5061...6601...8425...
  .34.171..508.1165..2286..4039..6616..10233..15130..21571..29844..40261...
  .55.341.1159.2929..6191.11605.19951..32129..49159..72181.102455.141361...
  .89.683.2683.7589.17621.35839.66263.113993.185329.287891.430739.624493...
  ...
Some solutions for n = 5 and k = 3 with colors = 1, 2, 3 and empty = 0:
..0....2....3....2....0....1....0....0....2....0....0....2....3....0....0....0
..0....2....3....2....2....1....2....3....2....1....0....2....3....1....1....1
..1....0....0....0....2....0....2....3....2....1....0....1....0....1....1....1
..1....2....2....0....3....2....2....3....2....0....3....1....3....3....2....1
..0....2....2....0....3....2....2....3....0....0....3....0....3....3....2....1
		

Crossrefs

Column 1 is A000045(n+1), column 2 is A001045(n+1), column 3 is A006130, column 4 is A006131, column 5 is A015440, column 6 is A015441(n+1), column 7 is A015442(n+1), column 8 is A015443, column 9 is A015445, column 10 is A015446, column 11 is A015447, and column 12 is A053404,
Row 2 is A000027(n+1), row 3 is A004273(n+1), row 4 is A028387, row 5 is A000567(n+1), and row 6 is A106734(n+2).
Diagonal is A171180, superdiagonal 1 is A083859(n+1), and superdiagonal 2 is A083860(n+1).

Programs

  • Maple
    T:= proc(n,k) option remember; `if`(n<0, 0,
          `if`(n<2 or k=0, 1, k*T(n-2, k) +T(n-1, k)))
        end;
    seq(seq(T(n, d+1-n), n=1..d), d=1..12); # Alois P. Heinz, Jul 29 2011
  • Mathematica
    T[n_, k_] := T[n, k] = If[n < 0, 0, If[n < 2 || k == 0, 1, k*T[n-2, k]+T[n-1, k]]]; Table[Table[T[n, d+1-n], {n, 1, d}], {d, 1, 12}] // Flatten (* Jean-François Alcover, Mar 04 2014, after Alois P. Heinz *)

Formula

With z X 1 tiles of k colors on an n X 1 grid (with n >= z), either there is a tile (of any of the k colors) on the first spot, followed by any configuration on the remaining (n-z) X 1 grid, or the first spot is vacant, followed by any configuration on the remaining (n-1) X 1. Thus, T(n,k) = T(n-1,k) + k*T(n-z,k), with T(n,k) = 1 for n = 0, 1, ..., z-1.
The solution is T(n,k) = Sum_r r^(-n-1)/(1 + z*k*r^(z-1)), where the sum is over the roots r of the polynomial k*x^z + x - 1.
For z = 2, T(n,k) = ((2*k / (sqrt(1 + 4*k) - 1))^(n+1) - (-2*k/(sqrt(1 + 4*k) + 1))^(n+1)) / sqrt(1 + 4*k).
T(n,k) = Sum_{s=0..[n/2]} binomial(n-s,s) * k^s.
For z X 1 tiles, T(n,k,z) = Sum_{s = 0..[n/z]} binomial(n-(z-1)*s, s) * k^s. - R. H. Hardin, Jul 31 2011

Extensions

Formula and proof from Robert Israel in the Sequence Fans mailing list.

A189800 a(n) = 6*a(n-1) + 8*a(n-2), with a(0)=0, a(1)=1.

Original entry on oeis.org

0, 1, 6, 44, 312, 2224, 15840, 112832, 803712, 5724928, 40779264, 290475008, 2069084160, 14738305024, 104982503424, 747801460736, 5326668791808, 37942424436736, 270267896954880, 1925146777223168, 13713023838978048, 97679317251653632, 695780094221746176
Offset: 0

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    I:=[0,1]; [n le 2 select I[n] else 6*Self(n-1)+8*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Nov 14 2011
    
  • Mathematica
    LinearRecurrence[{6, 8}, {0, 1}, 50]
    CoefficientList[Series[-(x/(-1+6 x+8 x^2)),{x,0,50}],x] (* Harvey P. Dale, Jul 26 2011 *)
  • PARI
    a(n)=([0,1; 8,6]^n*[0;1])[1,1] \\ Charles R Greathouse IV, Oct 03 2016

Formula

G.f.: x/(1 - 2*x*(3+4*x)). - Harvey P. Dale, Jul 26 2011

A060959 Table by antidiagonals of generalized Fibonacci numbers: T(n,k) = T(n,k-1) + n*T(n,k-2) with T(n,0)=0 and T(n,1)=1.

Original entry on oeis.org

0, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 2, 1, 1, 0, 1, 3, 3, 1, 1, 0, 1, 5, 5, 4, 1, 1, 0, 1, 8, 11, 7, 5, 1, 1, 0, 1, 13, 21, 19, 9, 6, 1, 1, 0, 1, 21, 43, 40, 29, 11, 7, 1, 1, 0, 1, 34, 85, 97, 65, 41, 13, 8, 1, 1, 0, 1, 55, 171, 217, 181, 96, 55, 15, 9, 1, 1, 0, 1, 89, 341, 508, 441, 301, 133, 71, 17, 10, 1, 1, 0
Offset: 0

Views

Author

Henry Bottomley, May 10 2001

Keywords

Examples

			Square array begins as:
  0, 1, 1, 1,  1,  1,  1, ...
  0, 1, 1, 2,  3,  5,  8, ...
  0, 1, 1, 3,  5, 11, 21, ...
  0, 1, 1, 4,  7, 19, 40, ...
  0, 1, 1, 5,  9, 29, 65, ...
  0, 1, 1, 6, 11, 41, 96, ...
		

Crossrefs

Programs

  • GAP
    Flat(List([0..12], n-> List([0..n], k-> (((1+Sqrt(1+4*k))/2)^(n-k) - ((1-Sqrt(1+4*k))/2)^(n-k))/Sqrt(1+4*k) ))); # G. C. Greubel, Jan 15 2020
  • Magma
    [Round( (((1+Sqrt(1+4*k))/2)^(n-k) - ((1-Sqrt(1+4*k))/2)^(n-k) )/Sqrt(1+4*k) ): k in [0..n], n in [0..12]]; // G. C. Greubel, Jan 15 2020
    
  • Maple
    seq(seq( round((((1+sqrt(1+4*k))/2)^(n-k) - ((1-sqrt(1+4*k))/2)^(n-k) )/sqrt(1+4*k)), k=0..n), n=0..12); # G. C. Greubel, Jan 15 2020
  • Mathematica
    T[n_, k_]:= If[n==k==0, 0, Round[(((1+Sqrt[1+4n])/2)^k - ((1-Sqrt[1+4n])/2)^k)/Sqrt[1+4n]]]; Table[T[k, n-k], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Jan 15 2020 *)
  • PARI
    T(n,k) = ( ((1+sqrt(1+4*n))/2)^k - ((1-sqrt(1+4*n))/2)^k )/sqrt(1+4*n);
    for(n=0,12, for(k=0,n, print1( round(T(k,n-k)), ", "))) \\ G. C. Greubel, Jan 15 2020
    
  • Sage
    [[ round( (((1+sqrt(1+4*k))/2)^(n-k) - ((1-sqrt(1+4*k))/2)^(n-k) )/sqrt(1+4*k) ) for k in (0..n)] for n in (0..12)] # G. C. Greubel, Jan 15 2020
    

Formula

T(n, k) = ( ((1+sqrt(1+4*n))/2)^k - ((1-sqrt(1+4*n))/2)^k )/sqrt(1+4*n).

A365824 a(n) = a(n-1) + 5*a(n-2), for n >= 0, with a(0) = 1 and a(1) = 0.

Original entry on oeis.org

1, 0, 5, 5, 30, 55, 205, 480, 1505, 3905, 11430, 30955, 88105, 242880, 683405, 1897805, 5314830, 14803855, 41378005, 115397280, 322287305, 899273705, 2510710230, 7007078755, 19560629905, 54596023680, 152399173205, 425379291605, 1187375157630, 3314271615655
Offset: 0

Views

Author

Wolfdieter Lang, Nov 20 2023

Keywords

Comments

This sequence {a(n)} appears in the formula for powers of phi21 := (1 + sqrt(21))/2 = A222134 = 2.791287..., together with b(n) = A015440(n-1), with A015440(-1) = 0, as phi21^n = a(n) + b(n)*phi21(n), for n >= 0. But the later given formulas in terms of scaled Chebyshev polynomials, called here {S21(n)}, are valid also for negative n values, i.e., for nonnegative powers of 1/phi21 = (-1 + sqrt(21))/10 = 0.35825756949... = A367453.
Limit_{n->oo} a(n)/a(n-1) = (1 + sqrt(21))/2 = A222134 = 2.791287...

Examples

			phi21^2 = a(2) + b(2)*phi(n) = 5 + phi21 = 7.79128784..., a real algebraic integer in Q(sqrt(21)).
(1/phi21)^2 = a(-2) + b(-2)*phi21 = (1/25)*(6 - phi21) = 0.12834848..., a real algebraic number in Q(sqrt(21)).
		

Crossrefs

Cf. A010477 (sqrt(21)), A015440, A049310, A222134, A367453.

Programs

  • Mathematica
    LinearRecurrence[{1,5},{1,0},50] (* Paolo Xausa, Nov 21 2023 *)
  • PARI
    a(n) = abs([1, 3; 1, -2]^(n-2)*[5; 5])[2, 1] \\ Thomas Scheuerle, Nov 20 2023

Formula

a(n) = a(n-1) + 5*a(n-2), for n >= 0, with a(0) = 1 and a(1) = 0.
G.f.: (1 - x)/(1 - x - 5*x^2).
a(n) = S21(n+1) - S21(n), for n >= 0, where S21(n) = sqrt(-5)^(n-1)*S(n-1, 1/sqrt(-5)), with the Chebyshev polynomials {S(n, x)} (see A049310).
The above mentioned sequence {b(n)} has terms b(n) = A015440(n-1) = S21(n), for n >= 0, with the same recurrence as {a(n)} but with b(0) = 0 and b(1) = 1, and g.f. x/(1 - x - 5*x^2).
The formula for negative indices of S is: S(-1, 0) = 0 and S(-n, x) = -S(n-2, x) for n >= 2.

A109447 Binomial coefficients C(n,k) with n-k odd, read by rows.

Original entry on oeis.org

1, 2, 1, 3, 4, 4, 1, 10, 5, 6, 20, 6, 1, 21, 35, 7, 8, 56, 56, 8, 1, 36, 126, 84, 9, 10, 120, 252, 120, 10, 1, 55, 330, 462, 165, 11, 12, 220, 792, 792, 220, 12, 1, 78, 715, 1716, 1287, 286, 13, 14, 364, 2002, 3432, 2002, 364, 14, 1, 105, 1365, 5005, 6435, 3003, 455, 15
Offset: 1

Views

Author

Philippe Deléham, Aug 27 2005

Keywords

Comments

The same as A119900 without 0's. A reflected version of A034867 or A202064. - Alois P. Heinz, Feb 07 2014
From Vladimir Shevelev, Feb 07 2014: (Start)
Also table of coefficients of polynomials P_1(x)=1, P_2(x)=2, for n>=2, P_(n+1)(x) = 2*P_n(x)+(x-1)* P_(n-1)(x). The polynomials P_n(x)/2^(n-1) are connected with sequences A000045 (x=5), A001045 (x=9), A006130 (x=13), A006131 (x=17), A015440 (x=21), A015441 (x=25), A015442 (x=29), A015443 (x=33), A015445 (x=37), A015446 (x=41), A015447 (x=45), A053404 (x=49); also the polynomials P_n(x) are connected with sequences A000129, A002605, A015518, A063727, A085449, A002532, A083099, A015519, A003683, A002534, A083102, A015520. (End)

Examples

			Starred terms in Pascal's triangle (A007318), read by rows:
1;
1*, 1;
1, 2*, 1;
1*, 3, 3*, 1;
1, 4*, 6, 4*, 1;
1*, 5, 10*, 10, 5*, 1;
1, 6*, 15, 20*, 15, 6*, 1;
1*, 7, 21*, 35, 35*, 21, 7*, 1;
1, 8*, 28, 56*, 70, 56*, 28, 8*, 1;
1*, 9, 36*, 84, 126*, 126, 84*, 36, 9*, 1;
Triangle T(n,k) begins:
1;
2;
1,    3;
4,    4;
1,   10,  5;
6,   20,  6;
1,   21,  35,   7;
8,   56,  56,   8;
1,   36, 126,  84,  9;
10, 120, 252, 120, 10;
		

Crossrefs

Cf. A109446.

Programs

  • Maple
    T:= (n, k)-> binomial(n, 2*k+1-irem(n, 2)):
    seq(seq(T(n, k), k=0..ceil((n-2)/2)), n=1..20);  # Alois P. Heinz, Feb 07 2014
  • Mathematica
    Flatten[ Table[ If[ OddQ[n - k], Binomial[n, k], {}], {n, 0, 15}, {k, 0, n}]] (* Robert G. Wilson v *)

Extensions

More terms from Robert G. Wilson v, Aug 30 2005
Corrected offset by Alois P. Heinz, Feb 07 2014

A128100 Triangle read by rows: T(n,k) is the number of ways to tile a 2 X n rectangle with k pieces of 2 X 2 tiles and n-2k pieces of 1 X 2 tiles (0 <= k <= floor(n/2)).

Original entry on oeis.org

1, 1, 2, 1, 3, 2, 5, 5, 1, 8, 10, 3, 13, 20, 9, 1, 21, 38, 22, 4, 34, 71, 51, 14, 1, 55, 130, 111, 40, 5, 89, 235, 233, 105, 20, 1, 144, 420, 474, 256, 65, 6, 233, 744, 942, 594, 190, 27, 1, 377, 1308, 1836, 1324, 511, 98, 7, 610, 2285, 3522, 2860, 1295, 315, 35, 1, 987, 3970
Offset: 0

Views

Author

Emeric Deutsch, Feb 18 2007

Keywords

Comments

Row sums are the Jacobsthal numbers (A001045). Column 0 yields the Fibonacci numbers (A000045); the other columns yield convolved Fibonacci numbers (A001629, A001628, A001872, A001873, etc.). Sum_{k=0..floor(n/2)} k*T(n,k) = A073371(n-2).
Triangle T(n,k), with zeros omitted, given by (1, 1, -1, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (0, 1, -1, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Jan 24 2012
Riordan array (1/(1-x-x^2), x^2/(1-x-x^2)), with zeros omitted. - Philippe Deléham, Feb 06 2012
Diagonal sums are A000073(n+2) (tribonacci numbers). - Philippe Deléham, Feb 16 2014
Number of induced subgraphs of the Fibonacci cube Gamma(n-1) that are isomorphic to the hypercube Q_k. Example: row n=4 is 5, 5, 1; indeed, the Fibonacci cube Gamma(3) is a square with an additional pendant edge attached to one of its vertices; it has 5 vertices (i.e., Q_0's), 5 edges (i.e., Q_1's) and 1 square (i.e., Q_2). - Emeric Deutsch, Aug 12 2014
Row n gives the coefficients of the polynomial p(n,x) defined as the numerator of the rational function given by f(n,x) = 1 + (x + 1)/f(n-1,x), where f(x,0) = 1. Conjecture: for n > 2, p(n,x) is irreducible if and only if n is a (prime - 2). - Clark Kimberling, Oct 22 2014

Examples

			Triangle starts:
   1;
   1;
   2,  1;
   3,  2;
   5,  5,  1;
   8, 10,  3;
  13, 20,  9,  1;
  21, 38, 22,  4;
From _Philippe Deléham_, Jan 24 2012: (Start)
Triangle (1, 1, -1, 0, 0, ...) DELTA (0, 1, -1, 0, 0, 0, ...) begins:
   1;
   1,  0;
   2,  1,  0;
   3,  2,  0,  0;
   5,  5,  1,  0,  0;
   8, 10,  3,  0,  0,  0;
  13, 20,  9,  1,  0,  0,  0;
  21, 38, 22,  4,  0,  0,  0,  0; (End)
From _Clark Kimberling_, Oct 22 2014: (Start)
Here are the first 4 polynomials p(n,x) as in Comment and generated by Mathematica program:
  1
  2 +  x
  3 + 2x
  5 + 5x + x^2. (End)
		

Crossrefs

Programs

  • Maple
    G:=1/(1-z-(1+t)*z^2): Gser:=simplify(series(G,z=0,19)): for n from 0 to 16 do P[n]:=sort(coeff(Gser,z,n)) od: for n from 0 to 16 do seq(coeff(P[n],t,j),j=0..floor(n/2)) od; # yields sequence in triangular form
  • Mathematica
    p[x_, n_] := 1 + (x + 1)/p[x, n - 1]; p[x_, 1] = 1;
    Numerator[Table[Factor[p[x, n]], {n, 1, 20}]]  (* Clark Kimberling, Oct 22 2014 *)

Formula

G.f.: 1/(1-z-(1+t)z^2).
Sum_{k=0..n} T(n,k)*x^k = A053404(n), A015447(n), A015446(n), A015445(n), A015443(n), A015442(n), A015441(n), A015440(n), A006131(n), A006130(n), A001045(n+1), A000045(n+1), A000012(n), A010892(n), A107920(n+1), A106852(n), A106853(n), A106854(n), A145934(n), A145976(n), A145978(n), A146078(n), A146080(n), A146083(n), A146084(n) for x = 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, -1, -2, -3, -4, -5, -6, -7, -8, -9, -10, -11, -12, and -13, respectively. - Philippe Deléham, Jan 24 2012
T(n,k) = T(n-1,k) + T(n-2,k) + T(n-2,k-1). - Philippe Deléham, Jan 24 2012
G.f.: T(0)/2, where T(k) = 1 + 1/(1 - (2*k+1+ x*(1+y))*x/((2*k+2+ x*(1+y))*x + 1/T(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Nov 06 2013
T(n,k) = Sum_{i=k..floor(n/2)} binomial(n-i,i)*binomial(i,k). See Corollary 3.3 in the Klavzar et al. link. - Emeric Deutsch, Aug 12 2014
Previous Showing 11-20 of 25 results. Next