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 31-40 of 439 results. Next

A289780 p-INVERT of the positive integers (A000027), where p(S) = 1 - S - S^2.

Original entry on oeis.org

1, 4, 14, 47, 156, 517, 1714, 5684, 18851, 62520, 207349, 687676, 2280686, 7563923, 25085844, 83197513, 275925586, 915110636, 3034975799, 10065534960, 33382471801, 110713382644, 367182309614, 1217764693607, 4038731742156, 13394504020957, 44423039068114
Offset: 0

Views

Author

Clark Kimberling, Aug 10 2017

Keywords

Comments

Suppose s = (c(0), c(1), c(2), ...) is a sequence and p(S) is a polynomial. Let S(x) = c(0)*x + c(1)*x^2 + c(2)*x^3 + ... and T(x) = (-p(0) + 1/p(S(x)))/x. The p-INVERT of s is the sequence t(s) of coefficients in the Maclaurin series for T(x).
Taking p(S) = 1 - S gives the INVERT transform of s, so that p-INVERT is a generalization of the INVERT transform (e.g., A033453).
Guide to p-INVERT sequences using p(S) = 1 - S - S^2:
t(A000012) = t(1,1,1,1,1,1,1,...) = A001906
t(A000290) = t(1,4,9,16,25,36,...) = A289779
t(A000027) = t(1,2,3,4,5,6,7,8,...) = A289780
t(A000045) = t(1,2,3,5,8,13,21,...) = A289781
t(A000032) = t(2,1,3,4,7,11,14,...) = A289782
t(A000244) = t(1,3,9,27,81,243,...) = A289783
t(A000302) = t(1,4,16,64,256,...) = A289784
t(A000351) = t(1,5,25,125,625,...) = A289785
t(A005408) = t(1,3,5,7,9,11,13,...) = A289786
t(A005843) = t(2,4,6,8,10,12,14,...) = A289787
t(A016777) = t(1,4,7,10,13,16,...) = A289789
t(A016789) = t(2,5,8,11,14,17,...) = A289790
t(A008585) = t(3,6,9,12,15,18,...) = A289795
t(A000217) = t(1,3,6,10,15,21,...) = A289797
t(A000225) = t(1,3,7,15,31,63,...) = A289798
t(A000578) = t(1,8,27,64,625,...) = A289799
t(A000984) = t(1,2,6,20,70,252,...) = A289800
t(A000292) = t(1,4,10,20,35,56,...) = A289801
t(A002620) = t(1,2,4,6,9,12,16,...) = A289802
t(A001906) = t(1,3,8,21,55,144,...) = A289803
t(A001519) = t(1,1,2,5,13,34,...) = A289804
t(A103889) = t(2,1,4,3,6,5,8,7,,...) = A289805
t(A008619) = t(1,1,2,2,3,3,4,4,...) = A289806
t(A080513) = t(1,2,2,3,3,4,4,5,...) = A289807
t(A133622) = t(1,2,1,3,1,4,1,5,...) = A289809
t(A000108) = t(1,1,2,5,14,42,...) = A081696
t(A081696) = t(1,1,3,9,29,97,...) = A289810
t(A027656) = t(1,0,2,0,3,0,4,0,5...) = A289843
t(A175676) = t(1,0,0,2,0,0,3,0,...) = A289844
t(A079977) = t(1,0,1,0,2,0,3,...) = A289845
t(A059841) = t(1,0,1,0,1,0,1,...) = A289846
t(A000040) = t(2,3,5,7,11,13,...) = A289847
t(A008578) = t(1,2,3,5,7,11,13,...) = A289828
t(A000142) = t(1!, 2!, 3!, 4!, ...) = A289924
t(A000201) = t(1,3,4,6,8,9,11,...) = A289925
t(A001950) = t(2,5,7,10,13,15,...) = A289926
t(A014217) = t(1,2,4,6,11,17,29,...) = A289927
t(A000045*) = t(0,1,1,2,3,5,...) = A289975 (* indicates prepended 0's)
t(A000045*) = t(0,0,1,1,2,3,5,...) = A289976
t(A000045*) = t(0,0,0,1,1,2,3,5,...) = A289977
t(A290990*) = t(0,1,2,3,4,5,...) = A290990
t(A290990*) = t(0,0,1,2,3,4,5,...) = A290991
t(A290990*) = t(0,0,01,2,3,4,5,...) = A290992

Examples

			Example 1:  s = (1,2,3,4,5,6,...) = A000027 and p(S) = 1 - S.
S(x) = x + 2x^2 + 3x^3 + 4x^4 + ...
p(S(x)) = 1 - (x + 2x^2 + 3x^3 + 4x^4 + ... )
- p(0) + 1/p(S(x)) = -1 + 1 + x + 3x^2 + 8x^3 + 21x^4 + ...
T(x) = 1 + 3x + 8x^2 + 21x^3 + ...
t(s) = (1,3,8,21,...) = A001906.
***
Example 2:  s = (1,2,3,4,5,6,...) = A000027 and p(S) = 1 - S - S^2.
S(x) =  x + 2x^2 + 3x^3 + 4x^4 + ...
p(S(x)) = 1 - ( x + 2x^2 + 3x^3 + 4x^4 + ...) - ( x + 2x^2 + 3x^3 + 4x^4 + ...)^2
- p(0) + 1/p(S(x)) = -1 + 1 + x + 4x^2 + 14x^3 + 47x^4 + ...
T(x) = 1 + 4x + 14x^2 + 47x^3 + ...
t(s) = (1,4,14,47,...) = A289780.
		

Crossrefs

Cf. A000027.

Programs

  • GAP
    P:=[1,4,14,47];; for n in [5..10^2] do P[n]:=5*P[n-1]-7*P[n-2]+5*P[n-3]-P[n-4]; od; P; # Muniru A Asiru, Sep 03 2017
  • Mathematica
    z = 60; s = x/(1 - x)^2; p = 1 - s - s^2;
    Drop[CoefficientList[Series[s, {x, 0, z}], x], 1] (* A000027 *)
    Drop[CoefficientList[Series[1/p, {x, 0, z}], x], 1] (* A289780 *)
  • PARI
    x='x+O('x^99); Vec((1-x+x^2)/(1-5*x+7*x^2-5*x^3+x^4)) \\ Altug Alkan, Aug 13 2017
    

Formula

G.f.: (1 - x + x^2)/(1 - 5 x + 7 x^2 - 5 x^3 + x^4).
a(n) = 5*a(n-1) - 7*a(n-2) + 5*a(n-3) - a(n-4).

A004254 a(n) = 5*a(n-1) - a(n-2) for n > 1, a(0) = 0, a(1) = 1.

Original entry on oeis.org

0, 1, 5, 24, 115, 551, 2640, 12649, 60605, 290376, 1391275, 6665999, 31938720, 153027601, 733199285, 3512968824, 16831644835, 80645255351, 386394631920, 1851327904249, 8870244889325, 42499896542376, 203629237822555, 975646292570399, 4674602225029440, 22397364832576801
Offset: 0

Views

Author

Keywords

Comments

Nonnegative values of y satisfying x^2 - 21*y^2 = 4; values of x are in A003501. - Wolfdieter Lang, Nov 29 2002
a(n) is equal to the permanent of the (n-1) X (n-1) Hessenberg matrix with 5's along the main diagonal, i's along the superdiagonal and the subdiagonal (i is the imaginary unit), and 0's everywhere else. - John M. Campbell, Jun 09 2011
For n >= 1, a(n) equals the number of 01-avoiding words of length n-1 on alphabet {0,1,2,3,4}. - Milan Janjic, Jan 25 2015
From Klaus Purath, Jul 26 2024: (Start)
For any three consecutive terms (x, y, z), y^2 - xz = 1 always applies.
a(n) = (t(i+2n) - t(i))/(t(i+n+1) - t(i+n-1)) where (t) is any recurrence t(k) = 4t(k-1) + 4t(k-2) - t(k-3) or t(k) = 5t(k-1) - t(k-2) without regard to initial values.
In particular, if the recurrence (t) of the form (4,4,-1) has the same three initial values as the current sequence, a(n) = t(n) applies.
a(n) = (t(k+1)*t(k+n) - t(k)*t(k+n+1))/(y^2 - xz) where (t) is any recurrence of the current family with signature (5,-1) and (x, y, z) are any three consecutive terms of (t), for integer k >= 0. (End)

Examples

			G.f. = x + 5*x^2 + 24*x^3 + 115*x^4 + 551*x^5 + 2640*x^6 + 12649*x^7 + ...
		

References

  • F. A. Haight, On a generalization of Pythagoras' theorem, pp. 73-77 of J. C. Butcher, editor, A Spectrum of Mathematics. Auckland University Press, 1971.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Partial sums of A004253.
Cf. A000027, A001906, A001353, A003501, A030221. a(n) = sqrt((A003501(n)^2 - 4)/21).
First differences of a(n) are in A004253, partial sums in A089817.
Cf. A004253.
INVERT transformation yields A001109. - R. J. Mathar, Sep 11 2008

Programs

  • Magma
    [ n eq 1 select 0 else n eq 2 select 1 else 5*Self(n-1)-Self(n-2): n in [1..30] ]; // Vincenzo Librandi, Aug 19 2011
  • Maple
    A004254:=1/(1-5*z+z**2); # Simon Plouffe in his 1992 dissertation
  • Mathematica
    a[n_]:=(MatrixPower[{{1,3},{1,4}},n].{{1},{1}})[[2,1]]; Table[a[n],{n,0,40}] (* Vladimir Joseph Stephan Orlovsky, Feb 19 2010 *)
    a[ n_] := ChebyshevU[2 n - 1, Sqrt[7]/2] / Sqrt[7]; (* Michael Somos, Jan 22 2017 *)
  • PARI
    {a(n) = subst(4*poltchebi(n+1) - 10*poltchebi(n), x, 5/2) / 21}; /* Michael Somos, Dec 04 2002 */
    
  • PARI
    {a(n) = imag((5 + quadgen(84))^n) / 2^(n-1)}; /* Michael Somos, Dec 04 2002 */
    
  • PARI
    {a(n) = polchebyshev(n - 1, 2, 5/2)}; /* Michael Somos, Jan 22 2017 */
    
  • PARI
    {a(n) = simplify( polchebyshev( 2*n - 1, 2, quadgen(28)/2) / quadgen(28))}; /* Michael Somos, Jan 22 2017 */
    
  • Sage
    [lucas_number1(n,5,1) for n in range(27)] # Zerinvary Lajos, Jun 25 2008
    

Formula

G.f.: x/(1-5*x+x^2).
a(n) = (((5+sqrt(21))/2)^n-((5-sqrt(21))/2)^n)/sqrt(21). - Barry E. Williams, Aug 29 2000
a(n) = S(2*n-1, sqrt(7))/sqrt(7) = S(n-1, 5); S(n, x)=U(n, x/2), Chebyshev polynomials of 2nd kind, A049310.
A003501(n) = sqrt(21*a(n)^2 + 4).
a(n) = Sum_{k=0..n-1} binomial(n+k, 2*k+1)*2^k. - Paul Barry, Nov 30 2004
[A004253(n), a(n)] = [1,3; 1,4]^n * [1,0]. - Gary W. Adamson, Mar 19 2008
a(n+1) = Sum_{k=0..n} Gegenbauer_C(n-k,k+1,2). - Paul Barry, Apr 21 2009
a(n+1) = Sum_{k=0..n} A101950(n,k)*4^k. - Philippe Deléham, Feb 10 2012
From Peter Bala, Dec 23 2012: (Start)
Product {n >= 1} (1 + 1/a(n)) = (1/3)*(3 + sqrt(21)).
Product {n >= 2} (1 - 1/a(n)) = (1/10)*(3 + sqrt(21)). (End)
From Michael Somos, Jan 22 2017: (Start)
A054493(2*n - 1) = 7 * a(n)^2 for all n in Z.
a(n) = -a(-n) for all n in Z.
0 = -1 + a(n)*(+a(n) - 5*a(n+1)) + a(n+1)*(+a(n+1)) for all n in Z. (End)
Limit_{n->oo} a(n+1)/a(n) = (5 + sqrt(21))/2 = A107905. - Wolfdieter Lang, Nov 15 2023
From Klaus Purath, Jul 26 2024: (Start)
a(n) = 4(a(n-1) + a(n-2)) - a(n-3).
a(n) = 6(a(n-1) - a(n-2)) + a(n-3).
In general, for all sequences of the form U(n) = P*U(n-1) - U(n-2) the following applies:
U(n) = (P-1)*U(n-1) + (P-1)*U(n-2) - U(n-3).
U(n) = (P+1)*U(n-1) - (P+1)*U(n-2) + U(n-3). (End)
a(n) = (5*a(n-1)+sqrt(21*a(n-1)^2+4))/2 for n>0. - Alexandru Petrescu, Apr 15 2025
From Peter Bala, May 22 2025: (Start)
Product_{n >= 1} ((a(2*n) + 1)/(a(2*n) - 1))^2 = 7/3.
Product_{n >= 1} ((a(2*n+1) + 1)/(a(2*n+1) - 1))^2 = 25/21.
The o.g.f. A(x) satisfies A(x) + A(-x) + 10*A(x)*A(-x) = 0. The o.g.f. for A097778 equals -1/x * A(sqrt(x))*A(-sqrt(x)). (End)
E.g.f.: 2*exp(5*x/2)*sinh(sqrt(21)*x/2)/sqrt(21). - Stefano Spezia, Jul 02 2025

A027306 a(n) = 2^(n-1) + ((1 + (-1)^n)/4)*binomial(n, n/2).

Original entry on oeis.org

1, 1, 3, 4, 11, 16, 42, 64, 163, 256, 638, 1024, 2510, 4096, 9908, 16384, 39203, 65536, 155382, 262144, 616666, 1048576, 2449868, 4194304, 9740686, 16777216, 38754732, 67108864, 154276028, 268435456, 614429672, 1073741824, 2448023843
Offset: 0

Views

Author

Keywords

Comments

Inverse binomial transform of A027914. Hankel transform (see A001906 for definition) is {1, 2, 3, 4, ..., n, ...}. - Philippe Deléham, Jul 21 2005
Number of walks of length n on a line that starts at the origin and ends at or above 0. - Benjamin Phillabaum, Mar 05 2011
Number of binary integers (i.e., with a leading 1 bit) of length n+1 which have a majority of 1-bits. E.g., for n+1=4: (1011, 1101, 1110, 1111) a(3)=4. - Toby Gottfried, Dec 11 2011
Number of distinct symmetric staircase walks connecting opposite corners of a square grid of side n > 1. - Christian Barrientos, Nov 25 2018
From Gus Wiseman, Aug 20 2021: (Start)
Also the number of integer compositions of n + 1 with alternating sum > 0, where the alternating sum of a sequence (y_1,...,y_k) is Sum_i (-1)^(i-1) y_i. These compositions are ranked by A345917. For example, the a(0) = 1 through a(4) = 11 compositions are:
(1) (2) (3) (4) (5)
(21) (31) (32)
(111) (112) (41)
(211) (113)
(122)
(212)
(221)
(311)
(1121)
(2111)
(11111)
The following relate to these compositions:
- The unordered version is A027193.
- The complement is counted by A058622.
- The reverse unordered version is A086543.
- The version for alternating sum >= 0 is A116406.
- The version for alternating sum < 0 is A294175.
- Ranked by A345917. (End)
The Gauss congruences a(n*p^k) == a(n^p^(k-1)) (mod p^k) hold for prime p and positive integers n and k. - Peter Bala, Jan 07 2022

Examples

			From _Gus Wiseman_, Aug 20 2021: (Start)
The a(0) = 1 through a(4) = 11 binary numbers with a majority of 1-bits (Gottfried's comment) are:
  1   11   101   1011   10011
           110   1101   10101
           111   1110   10110
                 1111   10111
                        11001
                        11010
                        11011
                        11100
                        11101
                        11110
                        11111
The version allowing an initial zero is A058622.
(End)
		

References

  • A. P. Prudnikov, Yu. A. Brychkov and O.I. Marichev, "Integrals and Series", Volume 1: "Elementary Functions", Chapter 4: "Finite Sums", New York, Gordon and Breach Science Publishers, 1986-1992, Eq. (4.2.1.6)

Crossrefs

a(n) = Sum{(k+1)T(n, m-k)}, 0<=k<=[ (n+1)/2 ], T given by A008315.
Column k=2 of A226873. - Alois P. Heinz, Jun 21 2013
The even bisection is A000302.
The odd bisection appears to be A032443.

Programs

  • GAP
    List([0..35],n->Sum([0..Int(n/2)],k->Binomial(n,k))); # Muniru A Asiru, Nov 27 2018
  • Haskell
    a027306 n = a008949 n (n `div` 2)  -- Reinhard Zumkeller, Nov 14 2014
    
  • Magma
    [2^(n-1)+(1+(-1)^n)/4*Binomial(n, n div 2): n in [0..40]]; // Vincenzo Librandi, Jun 19 2016
    
  • Maple
    a:= proc(n) add(binomial(n, j), j=0..n/2) end:
    seq(a(n), n=0..32); # Zerinvary Lajos, Mar 29 2009
  • Mathematica
    Table[Sum[Binomial[n, k], {k, 0, Floor[n/2]}], {n, 1, 35}]
    (* Second program: *)
    a[0] = a[1] = 1; a[2] = 3; a[n_] := a[n] = (2(n-1)(2a[n-2] + a[n-1]) - 8(n-2) a[n-3])/n; Array[a, 33, 0] (* Jean-François Alcover, Sep 04 2016 *)
  • PARI
    a(n)=if(n<0,0,(2^n+if(n%2,0,binomial(n, n/2)))/2)
    

Formula

a(n) = Sum_{k=0..floor(n/2)} binomial(n,k).
Odd terms are 2^(n-1). Also a(2n) - 2^(2n-1) is given by A001700. a(n) = 2^n + (n mod 2)*binomial(n, (n-1)/2).
E.g.f.: (exp(2x) + I_0(2x))/2.
O.g.f.: 2*x/(1-2*x)/(1+2*x-((1+2*x)*(1-2*x))^(1/2)). - Vladeta Jovovic, Apr 27 2003
a(n) = A008949(n, floor(n/2)); a(n) + a(n-1) = A248574(n), n > 0. - Reinhard Zumkeller, Nov 14 2014
From Peter Bala, Jul 21 2015: (Start)
a(n) = [x^n]( 2*x - 1/(1 - x) )^n.
O.g.f.: (1/2)*( 1/sqrt(1 - 4*x^2) + 1/(1 - 2*x) ).
Inverse binomial transform is (-1)^n*A246437(n).
exp( Sum_{n >= 1} a(n)*x^n/n ) = 1 + x + 2*x^2 + 3*x^3 + 6*x^4 + 10*x^5 + ... is the o.g.f. for A001405. (End)
a(n) = Sum_{k=1..floor((n+1)/2)} binomial(n-1,(2n+1-(-1)^n)/4 -k). - Anthony Browne, Jun 18 2016
D-finite with recurrence: n*a(n) + 2*(-n+1)*a(n-1) + 4*(-n+1)*a(n-2) + 8*(n-2)*a(n-3) = 0. - R. J. Mathar, Aug 09 2017

Extensions

Better description from Robert G. Wilson v, Aug 30 2000 and from Yong Kong (ykong(AT)curagen.com), Dec 28 2000

A027941 a(n) = Fibonacci(2*n + 1) - 1.

Original entry on oeis.org

0, 1, 4, 12, 33, 88, 232, 609, 1596, 4180, 10945, 28656, 75024, 196417, 514228, 1346268, 3524577, 9227464, 24157816, 63245985, 165580140, 433494436, 1134903169, 2971215072, 7778742048, 20365011073, 53316291172, 139583862444, 365435296161, 956722026040
Offset: 0

Views

Author

Keywords

Comments

Also T(2n+1,n+1), T given by A027935. Also first row of Inverse Stolarsky array.
Third diagonal of array defined by T(i, 1)=T(1, j)=1, T(i, j)=Max(T(i-1, j)+T(i-1, j-1); T(i-1, j-1)+T(i, j-1)). - Benoit Cloitre, Aug 05 2003
Number of Schroeder paths of length 2(n+1) having exactly one up step starting at an even height (a Schroeder path is a lattice path starting from (0,0), ending at a point on the x-axis, consisting only of steps U=(1,1) (up steps), D=(1,-1) (down steps) and H=(2,0) (level steps) and never going below the x-axis). Schroeder paths are counted by the large Schroeder numbers (A006318). Example: a(1)=4 because among the six Schroeder paths of length 4 only the paths (U)HD, (U)UDD, H(U)D, (U)DH have exactly one U step that starts at an even height (shown between parentheses). - Emeric Deutsch, Dec 19 2004
Also: smallest number not writeable as the sum of fewer than n positive Fibonacci numbers. E.g., a(5)=88 because it is the smallest number that needs at least 5 Fibonacci numbers: 88 = 55 + 21 + 8 + 3 + 1. - Johan Claes, Apr 19 2005 [corrected for offset and clarification by Mike Speciner, Sep 19 2023] In general, a(n) is the sum of n positive Fibonacci numbers as a(n) = Sum_{i=1..n} A000045(2*i). See A001076 when negative Fibonacci numbers can be included in the sum. - Mike Speciner, Sep 24 2023
Except for first term, numbers a(n) that set a new record in the number of Fibonacci numbers needed to sum up to n. Position of records in sequence A007895. - Ralf Stephan, May 15 2005
Successive extremal petal bends beta(n) = a(n-2). See the Ring Lemma of Rodin and Sullivan in K. Stephenson, Introduction to Circle Packing (Cambridge U. P., 2005), pp. 73-74 and 318-321. - David W. Cantrell (DWCantrell(AT)sigmaxi.net)
a(n+1)= AAB^(n)(1), n>=1, with compositions of Wythoff's complementary A(n):=A000201(n) and B(n)=A001950(n) sequences. See the W. Lang link under A135817 for the Wythoff representation of numbers (with A as 1 and B as 0 and the argument 1 omitted). E.g., 4=`110`, 12=`1100`, 33=`11000`, 88=`110000`, ..., in Wythoff code. AA(1)=1=a(1) but for uniqueness reason 1=A(1) in Wythoff code. - N. J. A. Sloane, Jun 29 2008
Start with n. Each n generates a sublist {n-1,n-1,n-2,..,1}. Each element of each sublist also generates a sublist. Add numbers in all terms. For example, 3->{2,2,1} and both 2->{1,1}, so a(3) = 3 + 2 + 2 + 1 + 1 + 1 + 1 + 1 = 12. - Jon Perry, Sep 01 2012
For n>0: smallest number such that the inner product of Zeckendorf binary representation and its reverse equals n: A216176(a(n)) = n, see also A189920. - Reinhard Zumkeller, Mar 10 2013
Also, numbers m such that 5*m*(m+2)+1 is a square. - Bruno Berselli, May 19 2014
Also, number of nonempty submultisets of multisets of weight n that span an initial interval of integers (see 2nd example). - Gus Wiseman, Feb 10 2015
From Robert K. Moniot, Oct 04 2020: (Start)
Including a(-1):=0, consecutive terms (a(n-1),a(n))=(u,v) or (v,u) give all points on the hyperbola u^2-u+v^2-v-4*u*v=0 with both coordinates nonnegative integers. Note that this follows from identifying (1,u+1,v+1) with the Markov triple (1,Fibonacci(2n-1),Fibonacci(2n+1)). See A001519 (comments by Robert G. Wilson, Oct 05 2005, and Wolfdieter Lang, Jan 30 2015).
Let T(n) denote the n-th triangular number. If i, j are any two successive elements of the above sequence then (T(i-1) + T(j-1))/T(i+j-1) = 3/5. (End)

Examples

			a(5) = 88 = 2*33 + 12 + 4 + 1 + 5. a(6) = 232 = 2*88 + 33 + 12 + 4 + 1 + 6. - _Jon Perry_, Sep 01 2012
a(4) = 33 counts all nonempty submultisets of the last row: [1][2][3][4], [11][12][13][14][22][23][24][33][34], [111][112][113][122][123][124][133][134][222][223][233][234], [1111][1112][1122][1123][1222][1223][1233][1234]. - _Gus Wiseman_, Feb 10 2015
		

References

  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 12.

Crossrefs

Related to partial sums of Fibonacci(k*n) over n: A000071, A099919, A058038, A138134, A053606; this sequence is the case k=2.
Cf. A212336 for more sequences with g.f. of the type 1/(1 - k*x + k*x^2 - x^3).
Cf. A000225 (sublist connection).
Cf. A258993 (row sums, n > 0), A000967.

Programs

Formula

a(n) = Sum_{i=1..n} binomial(n+i, n-i). - Benoit Cloitre, Oct 15 2002
G.f.: Sum_{k>=1} x^k/(1-x)^(2*k+1). - Benoit Cloitre, Apr 21 2003
a(n) = Sum_{k=1..n} F(2*k), i.e., partial sums of A001906. - Benoit Cloitre, Oct 27 2003
a(n) = Sum_{k=0..n-1} U(k, 3/2) = Sum_{k=0..n-1} S(k, 3), with S(k, 3) = A001906(k+1). - Paul Barry, Nov 14 2003
G.f.: x/((1-x)*(1-3*x+x^2)) = x/(1-4*x+4*x^2-x^3).
a(n) = 4*a(n-1) - 4*a(n-2) + a(n-3) with n>=2, a(-1)=0, a(0)=0, a(1)=1.
a(n) = 3*a(n-1) - a(n-2) + 1 with n>=1, a(-1)=0, a(0)=0.
a(n) = Sum_{k=1..n} F(k)*L(k), where L(k) = Lucas(k) = A000032(k) = F(k-1) + F(k+1). - Alexander Adamchuk, May 18 2007
a(n) = 2*a(n-1) + (Sum_{k=1..n-2} a(k)) + n. - Jon Perry, Sep 01 2012
Sum {n >= 1} 1/a(n) = 3 - phi, where phi = 1/2*(1 + sqrt(5)) is the golden ratio. The ratio of adjacent terms r(n) := a(n)/a(n-1) satisfies the recurrence r(n+1) = (4*r(n) - 1)/(r(n) + 1) for n >= 2. - Peter Bala, Dec 05 2013
a(n) = S(n, 3) - S(n-1, 3) - 1, n >= 0, with Chebyshev's S-polynomials (see A049310), where S(-1, x) = 0. - Wolfdieter Lang, Aug 28 2014
a(n) = -1 + (2^(-1-n)*((3-sqrt(5))^n*(-1+sqrt(5)) + (1+sqrt(5))*(3+sqrt(5))^n)) / sqrt(5). - Colin Barker, Jun 03 2016
E.g.f.: (sqrt(5)*sinh(sqrt(5)*x/2) + 5*cosh(sqrt(5)*x/2))*exp(3*x/2)/5 - exp(x). - Ilya Gutkovskiy, Jun 03 2016
a(n) = Sum_{k=0..n} binomial(n+1,k+1)*Fibonacci(k). - Vladimir Kruchinin, Oct 14 2016
a(n) = Sum_{k=0..n-1} Sum_{i=0..n-1} C(k+i+1,k-i). - Wesley Ivan Hurt, Sep 21 2017
a(n)*a(n-2) = a(n-1)*(a(n-1) - 1) for n>1. - Robert K. Moniot, Aug 23 2020
a(n) = Sum_{k=1..n} C(2*n-k,k). - Wesley Ivan Hurt, Dec 22 2020
a(n) = Sum_{k = 1..2*n+2} (-1)^k*Fibonacci(k). - Peter Bala, Nov 14 2021
a(n) = (2*cosh((1 + 2*n)*arccsch(2)))/sqrt(5) - 1. - Peter Luschny, Nov 21 2021
a(n) = F(n + (n mod 2)) * L(n+1 - (n mod 2)), where L(n) = A000032(n) and F(n) = A000045(n) (Euler and Sadek, 2001). - Amiram Eldar, Jan 13 2022

Extensions

More terms from James Sellers, Sep 08 2000
Paul Barry's Nov 14 2003 formula, recurrences and g.f. corrected for offset 0 and index link for Chebyshev polynomials added by Wolfdieter Lang, Aug 28 2014

A053763 a(n) = 2^(n^2 - n).

Original entry on oeis.org

1, 1, 4, 64, 4096, 1048576, 1073741824, 4398046511104, 72057594037927936, 4722366482869645213696, 1237940039285380274899124224, 1298074214633706907132624082305024, 5444517870735015415413993718908291383296, 91343852333181432387730302044767688728495783936
Offset: 0

Views

Author

Stephen G Penrice, Mar 29 2000

Keywords

Comments

Nilpotent n X n matrices over GF(2). Also number of simple digraphs (without self-loops) on n labeled nodes (see also A002416).
For n >= 1 a(n) is the size of the Sylow 2-subgroup of the Chevalley group A_n(4) (sequence A053291). - Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 30 2001
(-1)^ceiling(n/2) * resultant of the Chebyshev polynomial of first kind of degree n and Chebyshev polynomial of first kind of degree (n+1) (cf. A039991). - Benoit Cloitre, Jan 26 2003
The number of reflexive binary relations on an n-element set. - Justin Witt (justinmwitt(AT)gmail.com), Jul 12 2005
From Rick L. Shepherd, Dec 24 2008: (Start)
Number of gift exchange scenarios where, for each person k of n people,
i) k gives gifts to g(k) of the others, where 0 <= g(k) <= n-1,
ii) k gives no more than one gift to any specific person,
iii) k gives no single gift to two or more people and
iv) there is no other person j such that j and k jointly give a single gift.
(In other words -- but less precisely -- each person k either gives no gifts or gives exactly one gift per person to 1 <= g(k) <= n-1 others.) (End)
In general, sequences of the form m^((n^2 - n)/2) enumerate the graphs with n labeled nodes with m types of edge. a(n) therefore is the number of labeled graphs with n nodes with 4 types of edge. To clarify the comment from Benoit Cloitre, dated Jan 26 2003, in this context: simple digraphs (without self-loops) have four types of edge. These types of edges are as follows: the absent edge, the directed edge from A -> B, the directed edge from B -> A and the bidirectional edge, A <-> B. - Mark Stander, Apr 11 2019

Examples

			a(2)=4 because there are four 2 x 2 nilpotent matrices over GF(2):{{0,0},{0,0}},{{0,1},{0,0}},{{0,0},{1,0}},{{1,1,},{1,1}} where 1+1=0. - _Geoffrey Critzer_, Oct 05 2012
		

References

  • J. L. Gross and J. Yellen, eds., Handbook of Graph Theory, CRC Press, 2004; p. 521.
  • F. Harary and E. M. Palmer, Graphical Enumeration, Academic Press, NY, 1973, p. 5, Eq. (1.1.5).

Crossrefs

Programs

Formula

Sequence given by the Hankel transform (see A001906 for definition) of A059231 = {1, 1, 5, 29, 185, 1257, 8925, 65445, 491825, ...}; example: det([1, 1, 5, 29; 1, 5, 29, 185; 5, 29, 185, 1257; 29, 185, 1257, 8925]) = 4^6 = 4096. - Philippe Deléham, Aug 20 2005
a(n) = 4^binomial(n, n-2). - Zerinvary Lajos, Jun 16 2007
a(n) = Sum_{i=0..n^2-n} binomial(n^2-n, i). - Rick L. Shepherd, Dec 24 2008
G.f. A(x) satisfies: A(x) = 1 + x * A(4*x). - Ilya Gutkovskiy, Jun 04 2020
Sum_{n>=1} 1/a(n) = A319016. - Amiram Eldar, Oct 27 2020
Sum_{n>=0} a(n)*u^n/A002884(n) = Product_{r>=1} 1/(1-u/q^r). - Geoffrey Critzer, Oct 28 2021

A133494 Diagonal of the array of iterated differences of A047848.

Original entry on oeis.org

1, 1, 3, 9, 27, 81, 243, 729, 2187, 6561, 19683, 59049, 177147, 531441, 1594323, 4782969, 14348907, 43046721, 129140163, 387420489, 1162261467, 3486784401, 10460353203, 31381059609, 94143178827, 282429536481, 847288609443, 2541865828329, 7625597484987, 22876792454961, 68630377364883
Offset: 0

Views

Author

Paul Barry, Paul Curtz, Dec 23 2007

Keywords

Comments

a(n) is the number of ways to choose a composition C, and then choose a composition of each part of C. - Geoffrey Critzer, Mar 19 2012
a(n) is the top left entry of the n-th power of the 3 X 3 matrix [1, 1, 1; 1, 1, 1; 1, 1, 1]. - R. J. Mathar, Feb 03 2014
a(n) is the reptend length of 1/3^(n+1) in decimal. - Jianing Song, Nov 14 2018
Also the number of pairs of integer compositions, the first summing to n and the second with sum equal to the length of the first. If an integer composition is regarded as an arrow from sum to length, these are composable pairs, and the obvious composition operation founds a category of integer compositions. For example, we have (2,1,1,4) . (1,2,1) . (1,2) = (2,6), where dots represent the composition operation. The version without empty compositions is A000244. Composable triples are counted by 1 followed by A000302. The unordered version is A022811. - Gus Wiseman, Jul 14 2022

Examples

			From _Gus Wiseman_, Jul 15 2020: (Start)
The a(0) = 1 through a(3) = 9 ways to choose a composition of each part of a composition:
  ()  (1)  (2)      (3)
           (1,1)    (1,2)
           (1),(1)  (2,1)
                    (1,1,1)
                    (1),(2)
                    (2),(1)
                    (1),(1,1)
                    (1,1),(1)
                    (1),(1),(1)
(End)
		

Crossrefs

The strict version is A336139.
Splittings of partitions are A323583.
Multiset partitions of partitions are A001970.
Partitions of each part of a partition are A063834.
Compositions of each part of a partition are A075900.
Strict partitions of each part of a strict partition are A279785.
Compositions of each part of a strict partition are A304961.
Strict compositions of each part of a composition are A307068.
Compositions of each part of a strict composition are A336127.

Programs

Formula

Binomial transform of A078008. - Paul Curtz, Aug 04 2008
From R. J. Mathar, Nov 11 2008: (Start)
G.f.: (1 - 2*x)/(1 - 3*x).
a(n) = A000244(n-1), n > 0. (End)
From Philippe Deléham, Nov 13 2008: (Start)
a(n) = Sum_{k=0..n} A112467(n,k)*2^k.
a(n) = Sum_{k=0..n} A071919(n,k)*2^k. (End)
Let A(x) be the g.f. Then B(x) = x*A(x) satisfies B(x/(1-x)) = x/(1 - 2*B(x)). - Vladimir Kruchinin, Dec 05 2011
G.f.: 1/(1 - (Sum_{k>=1} (x/(1 - x))^k)). - Joerg Arndt, Sep 30 2012
For n > 0, a(n) = 2*(Sum_{k=0..n-1} a(k)) - 1 = 3^(n-1). - J. Conrad, Oct 29 2015
G.f.: 1 + x/(1 + x)*(1 + 4*x/(1 + 4*x)*(1 + 7*x/(1 + 7*x)*(1 + 10*x/(1 + 10*x)*(1 + .... - Peter Bala, May 27 2017
Invert transform of A011782(n) = 2^(n-1). Second invert transform of A000012. - Gus Wiseman, Jul 19 2020
a(n) = ceiling(3^(n-1)). - Alois P. Heinz, Jul 26 2020
From Elmo R. Oliveira, Mar 31 2025: (Start)
E.g.f.: (2 + exp(3*x))/3.
a(n) = 3*a(n-1) for n > 1. (End)

Extensions

Definition clarified by R. J. Mathar, Nov 11 2008

A004187 a(n) = 7*a(n-1) - a(n-2) with a(0) = 0, a(1) = 1.

Original entry on oeis.org

0, 1, 7, 48, 329, 2255, 15456, 105937, 726103, 4976784, 34111385, 233802911, 1602508992, 10983760033, 75283811239, 516002918640, 3536736619241, 24241153416047, 166151337293088, 1138818207635569, 7805576116155895, 53500214605455696, 366695926122033977
Offset: 0

Views

Author

Keywords

Comments

Define the sequence T(a_0,a_1) by a_{n+2} is the greatest integer such that a_{n+2}/a_{n+1}= 0 . A004187 (with initial 0 omitted) is T(1,7).
This is a divisibility sequence.
For n>=2, a(n) equals the permanent of the (n-1) X (n-1) tridiagonal matrix with 7's along the main diagonal, and i's along the superdiagonal and the subdiagonal (i is the imaginary unit). - John M. Campbell, Jul 08 2011
a(n) and b(n) := A056854(n) are the proper and improper nonnegative solutions of the Pell equation b(n)^2 - 5*(3*a(n))^2 = +4. see the cross-reference to A056854 below. - Wolfdieter Lang, Jun 26 2013
For n>=1, a(n) equals the number of 01-avoiding words of length n-1 on alphabet {0,1,2,3,4,5,6}. - Milan Janjic, Jan 25 2015
The digital root is A253298, which shares its digital root with A253368. - Peter M. Chema, Jul 04 2016
Lim_{n->oo} a(n+1)/a(n) = 2 + 3*phi = 1+ A090550 = 6.854101... - Wolfdieter Lang, Nov 16 2023

Examples

			a(2) = 7*a(1) - a(0) = 7*7 - 1 = 48. - _Michael B. Porter_, Jul 04 2016
		

Crossrefs

Cf. A000027, A001906, A001353, A004254, A001109, A049685, A033888. a(n)=sqrt((A056854(n)^2 - 4)/45).
Second column of array A028412.

Programs

  • Magma
    [Fibonacci(4*n)/3 : n in [0..30]]; // Vincenzo Librandi, Jun 07 2011
    
  • Magma
    /* By definition: */ [n le 2 select n-1 else 7*Self(n-1)-Self(n-2): n in [1..23]]; // Bruno Berselli, Dec 24 2012
  • Maple
    seq(combinat:-fibonacci(4*n)/3, n = 0 .. 30); # Robert Israel, Jan 26 2015
  • Mathematica
    LinearRecurrence[{7,-1},{0,1},30] (* Harvey P. Dale, Jul 13 2011 *)
    CoefficientList[Series[x/(1 - 7*x + x^2), {x, 0, 50}], x] (* Vincenzo Librandi, Dec 23 2012 *)
  • Maxima
    a[0]:0$ a[1]:1$ a[n]:=7*a[n-1] - a[n-2]$ A004187(n):=a[n]$ makelist(A004187(n),n,0,30); /* Martin Ettl, Nov 11 2012 */
    
  • MuPAD
    numlib::fibonacci(4*n)/3 $ n = 0..25; // Zerinvary Lajos, May 09 2008
    
  • PARI
    a(n)=fibonacci(4*n)/3 \\ Charles R Greathouse IV, Mar 09 2012
    
  • PARI
    concat(0, Vec(x/(1-7*x+x^2) + O(x^99))) \\ Altug Alkan, Jul 03 2016
    
  • Sage
    [lucas_number1(n,7,1) for n in range(27)] # Zerinvary Lajos, Jun 25 2008
    
  • Sage
    [fibonacci(4*n)/3 for n in range(0, 21)] # Zerinvary Lajos, May 15 2009
    

Formula

G.f.: x/(1-7*x+x^2).
a(n) = F(4*n)/3 = A033888(n)/3, where F=A000045 (the Fibonacci sequence).
a(n) = S(2*n-1, sqrt(9))/sqrt(9) = S(n-1, 7); S(n, x) := U(n, x/2), Chebyshev polynomials of the 2nd kind, A049310.
a(n) = Sum_{i = 0..n-1} C(2*n-1-i, i)*5^(n-i-1). - Mario Catalani (mario.catalani(AT)unito.it), Jul 23 2004
[A049685(n-1), a(n)] = [1,5; 1,6]^n * [1,0]. - Gary W. Adamson, Mar 21 2008
a(n) = A167816(4*n). - Reinhard Zumkeller, Nov 13 2009
a(n) = (((7+sqrt(45))/2)^n-((7-sqrt(45))/2)^n)/sqrt(45). - Noureddine Chair, Aug 31 2011
a(n+1) = Sum_{k = 0..n} A101950(n,k)*6^k. - Philippe Deléham, Feb 10 2012
a(n) = (A081072(n)/3)-1. - Martin Ettl, Nov 11 2012
From Peter Bala, Dec 23 2012: (Start)
Product {n >= 1} (1 + 1/a(n)) = (1/5)*(5 + 3*sqrt(5)).
Product {n >= 2} (1 - 1/a(n)) = (1/14)*(5 + 3*sqrt(5)). (End)
From Peter Bala, Apr 02 2015: (Start)
Sum_{n >= 1} a(n)*x^(2*n) = -A(x)*A(-x), where A(x) = Sum_{n >= 1} Fibonacci(2*n)* x^n.
1 + 5*Sum_{n >= 1} a(n)*x^(2*n) = F(x)*F(-x) = G(x)*G(-x), where F(x) = 1 + A(x) and G(x) = 1 + 5*A(x).
1 + Sum_{n >= 1} a(n)*x^(2*n) = H(x)*H(-x) = I(x)*I(-x), where H(x) = 1 + Sum_{n >= 1} Fibonacci(2*n + 3)*x^n and I(x) = 1 + x + x*Sum_{n >= 1} Fibonacci(2*n - 1)*x^n. (End)
E.g.f.: 2*exp(7*x/2)*sinh(3*sqrt(5)*x/2)/(3*sqrt(5)). - Ilya Gutkovskiy, Jul 03 2016
a(n) = Sum_{k = 0..n-1} (-1)^(n+k+1)*9^k*binomial(n+k, 2*k+1). - Peter Bala, Jul 17 2023
a(n) = Sum_{k = 0..floor(n/2)} (-1)^k*7^(n-2*k)*binomial(n-k, k). - Greg Dresden, Aug 03 2024
From Peter Bala, Jul 22 2025: (Start)
The following products telescope:
Product {n >= 2} (1 + (-1)^n/a(n)) = (3/14)*(3 + sqrt(5)).
Product {n >= 1} (1 - (-1)^n/a(n)) = (1/3)*(3 + sqrt(5)).
Product_{n >= 1} (a(2*n) + 1)/(a(2*n) - 1) = (3/5)*sqrt(5). (End)

Extensions

Entry improved by comments from Michael Somos and Wolfdieter Lang, Aug 02 2000

A290890 p-INVERT of the positive integers, where p(S) = 1 - S^2.

Original entry on oeis.org

0, 1, 4, 11, 28, 72, 188, 493, 1292, 3383, 8856, 23184, 60696, 158905, 416020, 1089155, 2851444, 7465176, 19544084, 51167077, 133957148, 350704367, 918155952, 2403763488, 6293134512, 16475640049, 43133785636, 112925716859, 295643364940, 774004377960
Offset: 0

Views

Author

Clark Kimberling, Aug 15 2017

Keywords

Comments

Suppose s = (c(0), c(1), c(2), ...) is a sequence and p(S) is a polynomial. Let S(x) = c(0)*x + c(1)*x^2 + c(2)*x^3 + ... and T(x) = (-p(0) + 1/p(S(x)))/x. The p-INVERT of s is the sequence t(s) of coefficients in the Maclaurin series for T(x). Taking p(S) = 1 - S gives the "INVERT" transform of s, so that p-INVERT is a generalization of the "INVERT" transform (e.g., A033453).
Note that in A290890, s = (1,2,3,4,...); i.e., A000027(n+1) for n>=0, whereas in A290990, s = (0,1,2,3,4,...); i.e., A000027(n) for n>=0.
Guide to p-INVERT sequences using s = (1,2,3,4,5,...) = A000027:
p(S) t(1,2,3,4,5,...)
1 - S A001906
1 - S^2 A290890; see A113067 for signed version
1 - S^3 A290891
1 - S^4 A290892
1 - S^5 A290893
1 - S^6 A290894
1 - S^7 A290895
1 - S^8 A290896
1 - S - S^2 A289780
1 - S - S^3 A290897
1 - S - S^4 A290898
1 - S^2 - S^4 A290899
1 - S^2 - S^3 A290900
1 - S^3 - S^4 A290901
1 - 2S A052530; (1/2)*A052530 = A001353
1 - 3S A290902; (1/3)*A290902 = A004254
1 - 4S A003319; (1/4)*A003319 = A001109
1 - 5S A290903; (1/5)*A290903 = A004187
1 - 2*S^2 A290904; (1/2)*A290904 = A290905
1 - 3*S^2 A290906; (1/3)*A290906 = A290907
1 - 4*S^2 A290908; (1/4)*A290908 = A099486
1 - 5*S^2 A290909; (1/5)*A290909 = A290910
1 - 6*S^2 A290911; (1/6)*A290911 = A290912
1 - 7*S^2 A290913; (1/7)*A290913 = A290914
1 - 8*S^2 A290915; (1/8)*A290915 = A290916
(1 - S)^2 A290917
(1 - S)^3 A290918
(1 - S)^4 A290919
(1 - S)^5 A290920
(1 - S)^6 A290921
1 - S - 2*S^2 A290922
1 - 2*S - 2*S^2 A290923; (1/2)*A290923 = A290924
1 - 3*S - 2*S^2 A290925
(1 - S^2)^2 A290926
(1 - S^2)^3 A290927
(1 - S^3)^2 A290928
(1 - S)(1 - S^2) A290929
(1 - S^2)(1 - S^4) A290930
1 - 3 S + S^2 A291025
1 - 4 S + S^2 A291026
1 - 5 S + S^2 A291027
1 - 6 S + S^2 A291028
1 - S - S^2 - S^3 A291029
1 - S - S^2 - S^3 - S^4 A201030
1 - 3 S + 2 S^3 A291031
1 - S - S^2 - S^3 + S^4 A291032
1 - 6 S A291033
1 - 7 S A291034
1 - 8 S A291181
1 - 3 S + 2 S^3 A291031
1 - 3 S + 2 S^2 A291182
1 - 4 S + 2 S^3 A291183
1 - 4 S + 3 S^3 A291184

Examples

			(See the examples at A289780.)
		

Crossrefs

Cf. A000027, A113067, A289780, A113067 (signed version of same sequence).

Programs

  • Mathematica
    z = 60; s = x/(1 - x)^2; p = 1 - s^2;
    Drop[CoefficientList[Series[s, {x, 0, z}], x], 1] (* A000027 *)
    Drop[CoefficientList[Series[1/p, {x, 0, z}], x], 1] (* A290890 *)

Formula

G.f.: x/(1 - 4 x + 5 x^2 - 4 x^3 + x^4).
a(n) = 4*a(n-1) - 5*a(n-2) + 4*a(n-3) - a(n-4).

A104457 Decimal expansion of 1 + phi = phi^2 = (3 + sqrt(5))/2.

Original entry on oeis.org

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

Views

Author

Eric W. Weisstein, Mar 08 2005

Keywords

Comments

Only first term differs from the decimal expansion of phi.
Zelo extends work of D. Roy by showing that the square of the golden ratio is the optimal exponent of approximation by algebraic numbers of degree 4 with bounded denominator and trace. - Jonathan Vos Post, Mar 02 2009 (Cf. last sentence in the Zelo reference. - Joerg Arndt, Jan 04 2014)
Hawkes asks: "What two numbers are those whose product, difference of their squares, and the ratio or quotient of their cubes, are all equal to each other?". - Charles R Greathouse IV, Dec 11 2012
This is the case n=10 in (Gamma(1/n)/Gamma(3/n))*(Gamma((n-1)/n)/Gamma((n-3)/n)) = 1+2*cos(2*Pi/n). - Bruno Berselli, Dec 14 2012
An algebraic integer of degree 2, with minimal polynomial x^2 - 3x + 1. - Charles R Greathouse IV, Nov 12 2014 [The other root is 2 - phi = A132338 - Wolfdieter Lang, Aug 29 2022]
To eight digits: 5*(((Pi+1)/e)-1) = 2.61803395481182... - Dan Graham, Nov 21 2017
The ratio diagonal/side of the second smallest diagonal in a regular 10-gon. - Mohammed Yaseen, Nov 04 2020
phi^2/10 is the moment of inertia of a solid regular icosahedron with a unit mass and a unit edge length (see A341906). - Amiram Eldar, Jun 08 2021

Examples

			2.6180339887498948482045868343656381177203091798...
		

References

  • Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, 2003, Section 2.17.1, p. 153.
  • Alfred S. Posamentier, Math Charmers, Tantalizing Tidbits for the Mind, Prometheus Books, NY, 2003, pages 138-139.
  • Damien Roy. Diophantine Approximation in Small Degree. Centre de Recherches Mathématiques. CRM Proceedings and Lecture Notes. Volume 36 (2004), 269-285.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 45.

Crossrefs

2 + 2*cos(2*Pi/n): A116425 (n = 7), A332438 (n = 9), A296184 (n = 10), A019973 (n = 12).

Programs

Formula

Equals 2 + A094214 = 1 + A001622. - R. J. Mathar, May 19 2008
Satisfies these three equations: x-sqrt(x)-1 = 0; x-1/sqrt(x)-2 = 0; x^2-3*x+1 = 0. - Richard R. Forberg, Oct 11 2014
Equals the nested radical sqrt(phi^2+sqrt(phi^4+sqrt(phi^8+...))). For a proof, see A094885. - Stanislav Sykora, May 24 2016
From Christian Katzmann, Mar 19 2018: (Start)
Equals Sum_{n>=0} (5*(2*n)!+8*n!^2)/(2*n!^2*3^(2*n+1)).
Equals 3/2 + Sum_{n>=0} 5*(2*n)!/(2*n!^2*3^(2*n+1)). (End)
Equals 1/A132338 = 2*A239798 = 5*A229780. - Mohammed Yaseen, Nov 04 2020
Equals Product_{k>=1} 1 + 1/(phi + phi^k), where phi is the golden ratio (A001622) (Ohtsuka, 2018). - Amiram Eldar, Dec 02 2021
c^n = phi * A001906(n) + A001519(n), where c = phi^2. - Gary W. Adamson, Sep 08 2023
Equals lim_{n->oo} S(n, 3)/S(n-1, 3) with the S-Chebyshev polynomials (see A049310), S(3, n) = A000045(2*(n+1)) = A001906(n+1). - Wolfdieter Lang, Nov 15 2023
From Peter Bala, May 08 2024: (Start)
Constant c = 2 + 2*cos(2*Pi/5).
The linear fractional transformation z -> c - c/z has order 5, that is, z = c - c/(c - c/(c - c/(c - c/(c - c/z)))). (End)
Equals Product_{k>=1} (1 + 1/A032908(k)). - Amiram Eldar, Nov 28 2024

A004189 a(n) = 10*a(n-1) - a(n-2); a(0) = 0, a(1) = 1.

Original entry on oeis.org

0, 1, 10, 99, 980, 9701, 96030, 950599, 9409960, 93149001, 922080050, 9127651499, 90354434940, 894416697901, 8853812544070, 87643708742799, 867583274883920, 8588189040096401, 85014307126080090, 841554882220704499, 8330534515080964900, 82463790268588944501, 816307368170808480110
Offset: 0

Views

Author

Keywords

Comments

Indices of square numbers which are also generalized pentagonal numbers.
If t(n) denotes the n-th triangular number, t(A105038(n))=a(n)*a(n+1). - Robert Phillips (bobanne(AT)bellsouth.net), May 25 2008
The n-th term is a(n) = ((5+sqrt(24))^n - (5-sqrt(24))^n)/(2*sqrt(24)). - Sture Sjöstedt, May 31 2009
For n >= 2, a(n) equals the permanent of the (n-1) X (n-1) tridiagonal matrix with 10's along the main diagonal, and i's along the superdiagonal and the subdiagonal (i is the imaginary unit). - John M. Campbell, Jul 08 2011
a(n) and b(n) (A001079) are the nonnegative proper solutions of the Pell equation b(n)^2 - 6*(2*a(n))^2 = +1. See the cross reference to A001079 below. - Wolfdieter Lang, Jun 26 2013
For n >= 1, a(n) equals the number of 01-avoiding words of length n-1 on alphabet {0,1,...,9}. - Milan Janjic, Jan 25 2015
For n > 1, this also gives the number of (n-1)-decimal-digit numbers which avoid a particular two-digit number with distinct digits. For example, there are a(5) = 9701 4-digit numbers which do not include "39" as a substring; see Wikipedia link. - Charles R Greathouse IV, Jan 14 2016
All possible solutions for y in Pell equation x^2 - 24*y^2 = 1. The values for x are given in A001079. - Herbert Kociemba, Jun 05 2022
Dickson on page 384 gives the Diophantine equation "(20) 24x^2 + 1 = y^2" and later states "... three consecutive sets (x_i, y_i) of solutions of (20) or 2x^2 + 1 = 3y^2 satisfy x_{n+1} = 10x_n - x_{n-1}, y_{n+1} = 10y_n - y_{n-1} with (x_1, y_1) = (0, 1) or (1, 1), (x_2, y_2) = (1, 5) or (11, 9), respectively." The first set of values (x_n, y_n) = (A001079(n-1), a(n-1)). - Michael Somos, Jun 19 2023

Examples

			a(2)=10 and (3(-8)^2-(-8))/2=10^2, a(3)=99 and (3(81)^2-(81))/2=99^2. - _Michael Somos_, Sep 05 2006
G.f. = x + 10*x^2 + 99*x^3 + 980*x^4 + 9701*x^5 + 96030*x^6 + ...
		

References

  • L. E. Dickson, History of the Theory of Numbers, Vol. II, Diophantine Analysis. AMS Chelsea Publishing, Providence, Rhode Island, 1999, p. 384.

Crossrefs

Chebyshev sequence U(n, m): A000027 (m=1), A001353 (m=2), A001109 (m=3), A001090 (m=4), this sequence (m=5), A004191 (m=6), A007655 (m=7), A077412 (m=8), A049660 (m=9), A075843 (m=10), A077421 (m=11), A077423 (m=12), A097309 (m=13), A097311 (m=14), A097313 (m=15), A029548 (m=16), A029547 (m=17), A144128 (m=18), A078987 (m=19), A097316 (m=33).
Cf. A323182.

Programs

  • GAP
    m:=5;; a:=[0,1];; for n in [3..20] do a[n]:=2*m*a[n-1]-a[n-2]; od; a; # G. C. Greubel, Dec 23 2019
  • Magma
    [ n eq 1 select 0 else n eq 2 select 1 else 10*Self(n-1)-Self(n-2): n in [1..20] ]; // Vincenzo Librandi, Aug 19 2011
    
  • Maple
    A004189 := proc(n)
        option remember;
        if n <= 1 then
            n ;
        else
            10*procname(n-1)-procname(n-2) ;
        end if;
    end proc:
    seq(A004189(n),n=0..20) ; # R. J. Mathar, Apr 30 2017
    seq( simplify(ChebyshevU(n-1, 5)), n=0..20); # G. C. Greubel, Dec 23 2019
  • Mathematica
    Table[GegenbauerC[n-1,1,5], {n,0,30}] (* Vladimir Joseph Stephan Orlovsky, Sep 11 2008; modified by G. C. Greubel, Jun 06 2019 *)
    LinearRecurrence[{10, -1}, {0, 1}, 20] (* Jean-François Alcover, Nov 15 2017 *)
    ChebyshevU[Range[21] -2, 5] (* G. C. Greubel, Dec 23 2019 *)
  • PARI
    {a(n) = subst(poltchebi(n+1) - 5*poltchebi(n), 'x, 5) / 24}; /* Michael Somos, Sep 05 2006 */
    
  • PARI
    a(n)=([9,1;8,1]^(n-1)*[1;1])[1,1] \\ Charles R Greathouse IV, Jan 14 2016
    
  • PARI
    vector(21, n, n--; polchebyshev(n-1, 2, 5) ) \\ G. C. Greubel, Dec 23 2019
    
  • Sage
    [lucas_number1(n,10,1) for n in range(22)] # Zerinvary Lajos, Jun 25 2008
    
  • Sage
    [chebyshev_U(n-1,5) for n in (0..20)] # G. C. Greubel, Dec 23 2019
    

Formula

a(n) = S(2*n-1, sqrt(12))/sqrt(12) = S(n-1, 10); S(n, x) := U(n, x/2), Chebyshev polynomials of 2nd kind, A049310. S(-1, x) := 0.
A001079(n) = sqrt(24*(a(n)^2)+1), that is a(n) = sqrt((A001079(n)^2-1)/24).
From Barry E. Williams, Aug 18 2000: (Start)
a(n) = ( (5+2*sqrt(6))^n - (5-2*sqrt(6))^n )/(4*sqrt(6)).
G.f.: x/(1-10*x+x^2). (End)
a(-n) = -a(n). - Michael Somos, Sep 05 2006
From Mohamed Bouhamida, May 26 2007: (Start)
a(n) = 9*(a(n-1) + a(n-2)) - a(n-3).
a(n) = 11*(a(n-1) - a(n-2)) + a(n-3).
a(n) = 10*a(n-1) - a(n-2). (End)
a(n+1) = Sum_{k=0..n} A101950(n,k)*9^k. - Philippe Deléham, Feb 10 2012
From Peter Bala, Dec 23 2012: (Start)
Product {n >= 1} (1 + 1/a(n)) = 1/2*(2 + sqrt(6)).
Product {n >= 2} (1 - 1/a(n)) = 1/5*(2 + sqrt(6)). (End)
a(n) = (A054320(n-1) + A072256(n))/2. - Richard R. Forberg, Nov 21 2013
a(2*n - 1) = A046173(n).
E.g.f.: exp(5*x)*sinh(2*sqrt(6)*x)/(2*sqrt(6)). - Stefano Spezia, Dec 12 2022
a(n) = Sum_{k = 0..n-1} binomial(n+k, 2*k+1)*8^k = Sum_{k = 0..n-1} (-1)^(n+k+1)* binomial(n+k, 2*k+1)*12^k. - Peter Bala, Jul 18 2023
Previous Showing 31-40 of 439 results. Next