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

A006190 a(n) = 3*a(n-1) + a(n-2), with a(0)=0, a(1)=1.

Original entry on oeis.org

0, 1, 3, 10, 33, 109, 360, 1189, 3927, 12970, 42837, 141481, 467280, 1543321, 5097243, 16835050, 55602393, 183642229, 606529080, 2003229469, 6616217487, 21851881930, 72171863277, 238367471761, 787274278560, 2600190307441, 8587845200883, 28363725910090
Offset: 0

Views

Author

Keywords

Comments

Denominators of continued fraction convergents to (3+sqrt(13))/2. - Benoit Cloitre, Jun 14 2003
a(n) and A006497(n) occur in pairs: (a,b): (1,3), (3,11), (10,36), (33,119), (109,393), ... such that b^2 - 13a^2 = 4(-1)^n. - Gary W. Adamson, Jun 15 2003
Form the 4-node graph with matrix A = [1,1,1,1; 1,1,1,0; 1,1,0,1; 1,0,1,1]. Then this sequence counts the walks of length n from the vertex with degree 5 to one (any) of the other vertices. - Paul Barry, Oct 02 2004
a(n+1) is the binomial transform of A006138. - Paul Barry, May 21 2006
a(n+1) is the diagonal sum of the exponential Riordan array (exp(3x),x). - Paul Barry, Jun 03 2006
Number of paths in the right half-plane from (0,0) to the line x=n-1, consisting of steps U=(1,1), D=(1,-1), h=(1,0) and H=(2,0). Example: a(3)=10 because we have hh, H, UD, DU, hU, Uh, UU, hD, Dh and DD. - Emeric Deutsch, Sep 03 2007
Equals INVERT transform of A000129. Example: a(5) = 109 = (29, 12, 5, 2, 1) dot (1, 1, 3, 10, 33) = (29 + 12 + 15 + 20 + 33). - Gary W. Adamson, Aug 06 2010
For n >= 2, a(n) equals the permanent of the (n-1) X (n-1) tridiagonal matrix with 3's along the main diagonal, 1's along the superdiagonal and subdiagonal, and 0's everywhere else. - John M. Campbell, Jul 08 2011
These numbers could also be called "bronze Fibonacci numbers". Indeed, for n >= 1, F(n+1) = ceiling(phi*F(n)), if n is even and F(n+1) = floor(phi*F(n)), if n is odd, where phi is the golden ratio; analogously, for Pell numbers (A000129), or "silver Fibonacci numbers", P(n+1) = ceiling(delta*a(n)), if n is even and P(n+1) = floor(delta*a(n)), if n is odd, where delta = delta_S = 1 + sqrt(2) is the silver ratio. Here, for n >= 1, we have a(n+1) = ceiling(c*a(n)), if n is even and a(n+1) = floor(c*a(n)), if n is odd, where c = (3 + sqrt(13))/2 is the bronze ratio (cf. comment in A098316). - Vladimir Shevelev, Feb 23 2013
Let p(n,x) denote the Fibonacci polynomial, defined by p(1,x) = 1, p(2,x) = x, p(n,x) = x*p(n-1,x) + p(n-2,x). Let q(n,x) be the numerator polynomial of the rational function p(n, x + 1 + 1/x). Then q(n,1) = a(n). - Clark Kimberling, Nov 04 2013
The (1,1)-entry of the matrix A^n where A = [0,1,0; 1,2,1; 1,1,2]. - David Neil McGrath, Jul 18 2014
a(n+1) counts closed walks on K2, containing three loops on the other vertex. Equivalently the (1,1)-entry of A^(n+1) where the adjacency matrix of digraph is A = (0,1; 1,3). - David Neil McGrath, Oct 29 2014
For n >= 1, a(n) equals the number of words of length n-1 on alphabet {0,1,2,3} avoiding runs of zeros of odd lengths. - Milan Janjic, Jan 28 2015
With offset 1 is the INVERTi transform of A001076. - Gary W. Adamson, Jul 24 2015
Apart from the initial 0, this is the p-INVERT transform of (1,0,1,0,1,0,...) for p(S) = 1 - 3 S. See A291219. - Clark Kimberling, Sep 02 2017
From Rogério Serôdio, Mar 30 2018: (Start)
This is a divisibility sequence (i.e., if n|m then a(n)|a(m)).
gcd(a(n),a(n+k)) = a(gcd(n, k)) for all positive integers n and k. (End)
Numbers of straight-chain fatty acids involving oxo and/or hydroxy groups, if cis-/trans isomerism and stereoisomerism are neglected. - Stefan Schuster, Apr 04 2018
Number of 3-compositions of n restricted to odd parts (and allowed zeros); see Hopkins & Ouvry reference. - Brian Hopkins, Aug 17 2020
From Michael A. Allen, Jan 25 2023: (Start)
Also called the 3-metallonacci sequence; the g.f. 1/(1-k*x-x^2) gives the k-metallonacci sequence.
a(n+1) is the number of tilings of an n-board (a board with dimensions n X 1) using unit squares and dominoes (with dimensions 2 X 1) if there are 3 kinds of squares available. (End)
a(n) is the number of compositions of n when there are P(k) sorts of parts k, with k,n >= 1, P(k) = A000129(k) is the k-th Pell number (see example below). - Enrique Navarrete, Dec 15 2023

Examples

			From _Enrique Navarrete_, Dec 15 2023: (Start)
From the comment on compositions with Pell number of parts, A000129(k), there are A000129(1)=1 type of 1, A000129(2)=2 types of 2, A000129(3)=5 types of 3, A000129(4)=12 types of 4, A000129(5)=29 types of 5 and A000129(6)=70 types of 6.
The following table gives the number of compositions of n=6:
Composition, number of such compositions, number of compositions of this type:
  6,              1,     70;
  5+1,            2,     58;
  4+2,            2,     48;
  3+3,            1,     25;
  4+1+1,          3,     36;
  3+2+1,          6,     60;
  2+2+2,          1,      8;
  3+1+1+1,        4,     20;
  2+2+1+1,        6,     24;
  2+1+1+1+1,      5,     10;
  1+1+1+1+1+1,    1,      1;
for a total of a(6)=360 compositions of n=6. (End).
		

References

  • H. L. Abbott and D. Hanson, A lattice path problem, Ars Combin., 6 (1978), 163-178.
  • A. Brousseau, Fibonacci and Related Number Theoretic Tables. Fibonacci Association, San Jose, CA, 1972, p. 128.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • L.-N. Machaut, Query 3436, L'Intermédiaire des Mathématiciens, 16 (1909), 62-63. - N. J. A. Sloane, Mar 08 2022

Crossrefs

Row sums of Pascal's rhombus (A059317). Also row sums of triangle A054456(n, m).
Sequences with g.f. 1/(1-k*x-x^2) or x/(1-k*x-x^2): A000045 (k=1), A000129 (k=2), this sequence (k=3), A001076 (k=4), A052918 (k=5), A005668 (k=6), A054413 (k=7), A041025 (k=8), A099371 (k=9), A041041 (k=10), A049666 (k=11), A041061 (k=12), A140455 (k=13), A041085 (k=14), A154597 (k=15), A041113 (k=16), A178765 (k=17), A041145 (k=18), A243399 (k=19), A041181 (k=20).
Cf. A006497, A052906, A175182 (Pisano periods), A201001 (prime subsequence), A092936 (squares).
Cf. A243399.

Programs

  • GAP
    a:=[0,1];; for n in [3..30] do a[n]:=3*a[n-1]+a[n-2]; od; a; # Muniru A Asiru, Mar 31 2018
  • Haskell
    a006190 n = a006190_list !! n
    a006190_list = 0 : 1 : zipWith (+) (map (* 3) $ tail a006190_list) a006190_list
    -- Reinhard Zumkeller, Feb 19 2011
    
  • Magma
    [ n eq 1 select 0 else n eq 2 select 1 else 3*Self(n-1)+Self(n-2): n in [1..30] ]; // Vincenzo Librandi, Aug 19 2011
    
  • Maple
    a[0]:=0: a[1]:=1: for n from 2 to 35 do a[n]:= 3*a[n-1]+a[n-2] end do: seq(a[n],n=0..30); # Emeric Deutsch, Sep 03 2007
    A006190:=-1/(-1+3*z+z**2); # Simon Plouffe in his 1992 dissertation, without the leading 0
    seq(combinat[fibonacci](n,3),n=0..30); # R. J. Mathar, Dec 07 2011
  • Mathematica
    a[n_] := (MatrixPower[{{1, 3}, {1, 2}}, n].{{1}, {1}})[[2, 1]]; Table[ a[n], {n, -1, 24}] (* Robert G. Wilson v, Jan 13 2005 *)
    LinearRecurrence[{3,1},{0,1},30] (* or *) CoefficientList[Series[x/ (1-3x-x^2), {x,0,30}], x] (* Harvey P. Dale, Apr 20 2011 *)
    Table[If[n==0, a1=1; a0=0, a2=a1; a1=a0; a0=3*a1+a2], {n, 0, 30}] (* Jean-François Alcover, Apr 30 2013 *)
    Table[Fibonacci[n, 3], {n, 0, 30}] (* Vladimir Reshetnikov, May 08 2016 *)
  • PARI
    a(n)=if(n<1,0,contfracpnqn(vector(n,i,2+(i>1)))[2,1])
    
  • PARI
    a(n)=([1,3;1,2]^n)[2,1] \\ Charles R Greathouse IV, Mar 06 2014
    
  • PARI
    concat([0],Vec(x/(1-3*x-x^2)+O(x^30))) \\ Joerg Arndt, Apr 30 2013
    
  • Sage
    [lucas_number1(n,3,-1) for n in range(0, 30)] # Zerinvary Lajos, Apr 22 2009
    

Formula

G.f.: x/(1 - 3*x - x^2).
From Benoit Cloitre, Jun 14 2003: (Start)
a(3*n) = 2*A041019(5*n-1), a(3*n+1) = A041019(5*n), a(3*n+2) = A041019(5*n+3).
a(2*n) = 3*A004190(n-1); a(3*n) = 10*A041613(n-1) for n >= 1. (End)
From Gary W. Adamson, Jun 15 2003: (Start)
a(n-1) + a(n+1) = A006497(n).
A006497(n)^2 - 13*a(n)^2 = 4(-1)^n. (End)
a(n) = U(n-1, (3/2)i)(-i)^(n-1), i^2 = -1. - Paul Barry, Nov 19 2003
a(n) = Sum_{k=0..n-1} binomial(n-k-1,k)*3^(n-2*k-1). - Paul Barry, Oct 02 2004
a(n) = F(n, 3), the n-th Fibonacci polynomial evaluated at x=3.
Let M = {{0, 1}, {1, 3}}, v[1] = {0, 1}, v[n] = M.v[n - 1]; then a(n) = Abs[v[n][[1]]]. - Roger L. Bagula, May 29 2005 [Or a(n) = [M^(n+1)]{1,1}. - _L. Edson Jeffery, Aug 27 2013]
From Paul Barry, May 21 2006: (Start)
a(n+1) = Sum_{k=0..n} Sum_{j=0..n-k} C(k,j)*C(n-j,k)*2^(k-j).
a(n) = Sum_{k=0..n} Sum_{j=0..n-k} C(k,j)*C(n-j,k)*2^(n-j-k).
a(n+1) = Sum_{k=0..floor(n/2)} C(n-k,k)*3^(n-2*k).
a(n) = Sum_{k=0..n} C(k,n-k)*3^(2*k-n). (End)
E.g.f.: exp(3*x/2)*sinh(sqrt(13)*x/2)/(sqrt(13)/2). - Paul Barry, Jun 03 2006
a(n) = (ap^n - am^n)/(ap - am), with ap = (3 + sqrt(13))/2, am = (3 - sqrt(13))/2.
Let C = (3 + sqrt(13))/2 = exp arcsinh(3/2) = 3.3027756377... Then C^n, n > 0 = a(n)*(1/C) + a(n+1). Let X = the 2 X 2 matrix [0, 1; 1, 3]. Then X^n = [a(n-1), a(n); a(n), a(n+1)]. - Gary W. Adamson, Dec 21 2007
1/3 = 3/(1*10) + 3/(3*33) + 3/(10*109) + 3/(33*360) + 3/(109*1189) + ... . - Gary W. Adamson, Mar 16 2008
a(n) = ((3 + sqrt(13))^n - (3 - sqrt(13))^n)/(2^n*sqrt(13)). - Al Hakanson (hawkuu(AT)gmail.com), Jan 12 2009
a(p) == 13^((p-1)/2) mod p, for odd primes p. - Gary W. Adamson, Feb 22 2009
From Johannes W. Meijer, Jun 12 2010: (Start)
Limit_{k->oo} a(n+k)/a(k) = (A006497(n) + a(n)*sqrt(13))/2.
Limit_{n->oo} A006497(n)/a(n) = sqrt(13). (End)
Sum_{k>=1} (-1)^(k-1)/(a(k)*a(k+1)) = (sqrt(13)-3)/2. - Vladimir Shevelev, Feb 23 2013
From Vladimir Shevelev, Feb 24 2013: (Start)
(1) Expression a(n+1) via a(n): a(n+1) = (3*a(n) + sqrt(13*a(n)^2 + 4*(-1)^n))/2;
(2) a^2(n+1) - a(n)*a(n+2) = (-1)^n;
(3) Sum_{k=1..n} (-1)^(k-1)/(a(k)*a(k+1)) = a(n)/a(n+1);
(4) a(n)/a(n+1) = (sqrt(13)-3)/2 + r(n), where |r(n)| < 1/(a(n+1)*a(n+2)). (End)
a(n) = sqrt(13*(A006497(n))^2 + (-1)^(n-1)*52)/13. - Vladimir Shevelev, Mar 13 2013
Sum_{n >= 1} 1/( a(2*n) + 1/a(2*n) ) = 1/3; Sum_{n >= 1} 1/( a(2*n + 1) - 1/a(2*n + 1) ) = 1/9. - Peter Bala, Mar 26 2015
From Rogério Serôdio, Mar 30 2018: (Start)
Some properties:
(1) a(n)*a(n+1) = 3*Sum_{k=1..n} a(k)^2;
(2) a(n)^2 + a(n+1)^2 = a(2*n+1);
(3) a(n)^2 - a(n-2)^2 = 3*a(n-1)*(a(n) + a(n-2));
(4) a(m*(p+1)) = a(m*p)*a(m+1) + a(m*p-1)*a(m);
(5) a(n-k)*a(n+k) = a(n)^2 + (-1)^(n+k+1)*a(k)^2;
(6) a(2*n) = a(n)*(3*a(n) + 2*a(n-1));
(7) 3*Sum_{k=2..n+1} a(k)*a(k-1) is equal to a(n+1)^2 if n odd, and is equal to a(n+1)^2 - 1 if n is even;
(8) a(n) = alpha(k)*a(n-2*k+1) + a(n-4*k+2), where alpha(k) = ap^(2*k-1) + am^(2*k-1), with ap = (3 + sqrt(13))/2, am = (3 - sqrt(13))/2;
(9) 131|Sum_{k=n..n+9} a(k), for all positive n. (End)
From Kai Wang, Feb 10 2020: (Start)
a(n)^2 - a(n+r)*a(n-r) = (-1)^(n-r)*a(r)^2 - Catalan's identity.
arctan(1/a(2n)) - arctan(1/a(2n+2)) = arctan(a(2)/a(2n+1)).
arctan(1/a(2n)) = Sum_{m>=n} arctan(a(2)/a(2m+1)).
The same formula holds for Fibonacci numbers and Pell numbers. (End)
a(n+2) = 3^(n+1) + Sum_{k=0..n} a(k)*3^(n-k). - Greg Dresden and Gavron Campbell, Feb 22 2022
G.f. = 1/(1-Sum_{k>=1} P(k)*x^k), P(k)=A000129(k) (with a(0)=1). - Enrique Navarrete, Dec 17 2023
G.f.: x/(1 - 3*x - x^2) = Sum_{n >= 0} x^(n+1) *( Product_{k = 1..n} (m*k + 3 - m + x)/(1 + m*k*x) ) for arbitrary m (a telescoping series). - Peter Bala, May 08 2024
Sum_{n>=0} a(n)/phi^(3*n) = 1, where phi = A001622 is the golden ratio. - Diego Rattaggi, Apr 07 2025

Extensions

Second formula corrected by Johannes W. Meijer, Jun 02 2010

A099372 a(n) = A099371(n)^2.

Original entry on oeis.org

0, 1, 81, 6724, 558009, 46308025, 3843008064, 318923361289, 26466795978921, 2196425142889156, 182276820063821025, 15126779640154255921, 1255340433312739420416, 104178129185317217638609, 8645529381948016324584129, 717474760572500037722844100, 59541759598135555114671476169
Offset: 0

Views

Author

Wolfdieter Lang, Oct 18 2004

Keywords

Comments

See the comment in A099279. This is example a=9.
a(n+1) is the number of tilings of an n-board (a board with dimensions n X 1) using half-squares (1/2 X 1 pieces, always placed so that the shorter sides are horizontal) and (1/2,1/2)-fences if there are 9 kinds of half-square available. A (w,g)-fence is a tile composed of two w X 1 pieces separated horizontally by a gap of width g. a(n+1) also equals the number of tilings of an n-board using (1/4,1/4)-fences and (1/4,3/4)-fences if there are 9 kinds of (1/4,1/4)-fence available. - Michael A. Allen, Mar 21 2024

Crossrefs

Cf. other squares of k-metallonacci numbers (for k=1 to 10): A007598, A079291, A092936, A099279, A099365, A099366, A099367, A099369, this sequence, A099374.

Programs

  • Mathematica
    LinearRecurrence[{82,82,-1},{0,1,81},17] (* Stefano Spezia, Apr 06 2023 *)

Formula

a(n) = A099371(n)^2.
a(n) = 82*a(n-1) + 82*a(n-2) - a(n-3), n>=3; a(0)=0, a(1)=1, a(2)=81.
a(n) = 83*a(n-1) - a(n-2) - 2*(-1)^n, n>=2; a(0)=0, a(1)=1.
a(n) = 2*(T(n, 83/2)-(-1)^n)/85 with twice the Chebyshev polynomials of the first kind: 2*T(n, 83/2) = A099373(n).
G.f.: x*(1-x)/((1-83*x+x^2)*(1+x)) = x*(1-x)/(1-82*x-82*x^2+x^3).
E.g.f.: 2*exp(-x)*(exp(85*x/2)*cosh(9*sqrt(85)*x/2) - 1)/85. - Stefano Spezia, Apr 06 2023
a(n) = (1 - (-1)^n)/2 + 81*Sum_{r=1..n-1} r*a(n-r). - Michael A. Allen, Mar 21 2024
Product_{n>=2} (1 + (-1)^n/a(n)) = (9 + sqrt(85))/18 (Falcon, 2016, p. 189, eq. (3.1)). - Amiram Eldar, Dec 03 2024

A099279 Squares of A001076.

Original entry on oeis.org

0, 1, 16, 289, 5184, 93025, 1669264, 29953729, 537497856, 9645007681, 173072640400, 3105662519521, 55728852710976, 1000013686278049, 17944517500293904, 322001301319012225, 5778078906241926144, 103683419011035658369, 1860523463292399924496, 33385738920252162982561
Offset: 0

Views

Author

Wolfdieter Lang, Oct 18 2004

Keywords

Comments

For the generalized Fibonacci sequences U(n-1;a) = (ap(a)^n - am(a)^n)/(ap(a) - am(a)) with ap(a) = (a + sqrt(a^2+4))/2, am(a) = (a - sqrt(a^2+4))/2, a from the integers, one has for the squared sequences U(n-1;a)^2 = (2*T(n,(a^2+2)/2) - 2*(-1)^n)/(a^2+4). Here T(n,x) are Chebyshev's polynomials of the first kind (see A053120). Therefore the o.g.f. for the squared sequence is x*(1-x)/((1+x)*(1-(a^2+2)*x+x^2)) = x*(1-x)/(1 - (a^2+1)*x - (a^2+1)*x^2 + x^3). For this example a=4.
Unsigned member r=-16 of the family of Chebyshev sequences S_r(n) defined in A092184.
(-1)^(n+1)*a(n) = S_{-16}(n), n >= 0, defined in A092184.
a(n+1) is the number of tilings of an n-board (a board with dimensions n X 1) using half-squares (1/2 X 1 pieces, always placed so that the shorter sides are horizontal) and (1/2,1/2)-fences if there are 4 kinds of half-squares available. A (w,g)-fence is a tile composed of two w X 1 pieces separated horizontally by a gap of width g. a(n+1) also equals the number of tilings of an n-board using (1/4,1/4)-fences and (1/4,3/4)-fences if there are 4 kinds of (1/4,1/4)-fences available. - Michael A. Allen, Mar 12 2023

Crossrefs

Cf. other squares of k-metallonacci numbers (for k=1 to 10): A007598, A079291, A092936, this sequence, A099365, A099366, A099367, A099369, A099372, A099374.

Programs

  • Magma
    [Fibonacci(3*n)^2/4: n in [0..30]]; // G. C. Greubel, Aug 18 2022
  • Maple
    with (combinat):seq(fibonacci(n,4)^2,n=0..16); # Zerinvary Lajos, Apr 09 2008
    nmax:=48: with(combinat): for n from 0 to nmax do A001654(n):=fibonacci(n) * fibonacci(n+1) od: a(0):=0: for n from 1 to nmax/3 do a(n):=a(n-1)+A001654(3*n-2) od: seq(a(n),n=0..nmax/3); # Johannes W. Meijer, Sep 22 2010
  • Mathematica
    LinearRecurrence[{17,17,-1},{0,1,16},30] (* Harvey P. Dale, Mar 26 2012 *)
    Fibonacci[3*Range[0, 30]]^2/4 (* G. C. Greubel, Aug 18 2022 *)
  • MuPAD
    numlib::fibonacci(3*n)^2/4 $ n = 0..35; // Zerinvary Lajos, May 13 2008
    
  • PARI
    my(x='x+O('x^99)); concat([0], Vec(x*(1-x)/((1-18*x+x^2)*(1+x)))) \\ Altug Alkan, Dec 17 2017
    
  • Sage
    [(fibonacci(3*n))^2/4 for n in range(0, 17)] # Zerinvary Lajos, May 15 2009
    

Formula

a(n) = A001076(n)^2.
a(n) = 17*a(n-1) + 17*a(n-2) - a(n-3), n >= 3, a(0)=0, a(1)=1, a(2)=16.
a(n) = 18*a(n-1) - a(n-2) - 2*(-1)^n, n >= 2, a(0)=0, a(1)=1.
a(n) = (T(n, 9) - (-1)^n)/10 with Chebyshev's T(n, x) polynomials of the first kind. T(n, 9) = A023039(n).
G.f.: x*(1-x)/((1+x)*(1-18*x+x^2)) = x*(1-x)/(1-17*x-17*x^2+x^3).
a(n) = a(n-1) + A001654(3*n-2) with a(0)=0, where A001654 are the golden rectangle numbers. - Johannes W. Meijer, Sep 22 2010
a(n+1) = (1 + (-1)^n)/2 + 16*Sum_{r=1..n} ( r*a(n+1-r) ). - Michael A. Allen, Mar 12 2023
E.g.f.: exp(-x)*(exp(10*x)*cosh(4*sqrt(5)*x) - 1)/10. - Stefano Spezia, Apr 06 2023
Product_{n>=2} (1 + (-1)^n/a(n)) = (2 + sqrt(5))/4 (Falcon, 2016, p. 189, eq. (3.1)). - Amiram Eldar, Dec 03 2024

A099365 Squares of A052918(n-1) (generalized Fibonacci).

Original entry on oeis.org

0, 1, 25, 676, 18225, 491401, 13249600, 357247801, 9632441025, 259718659876, 7002771375625, 188815108482001, 5091005157638400, 137268324147754801, 3701153746831741225, 99793882840309258276, 2690733682941518232225
Offset: 0

Views

Author

Wolfdieter Lang, Oct 18 2004

Keywords

Comments

See the comment in A099279. This is example a=5.
a(n+1) is the number of tilings of an n-board (a board with dimensions n X 1) using half-squares (1/2 X 1 pieces, always placed so that the shorter sides are horizontal) and (1/2,1/2)-fences if there are 5 kinds of half-squares available. A (w,g)-fence is a tile composed of two w X 1 pieces separated horizontally by a gap of width g. a(n+1) also equals the number of tilings of an n-board using (1/4,1/4)-fences and (1/4,3/4)-fences if there are 5 kinds of (1/4,1/4)-fences available. - Michael A. Allen, Mar 30 2023

Crossrefs

Cf. other squares of k-metallonacci numbers (for k=1 to 10): A007598, A079291, A092936, A099279, this sequence, A099366, A099367, A099369, A099372, A099374.

Programs

  • Magma
    [(2/29)*(Evaluate(ChebyshevFirst(n), 27/2) -(-1)^n): n in [0..30]]; // G. C. Greubel, Aug 21 2022
    
  • Maple
    with (combinat):seq(fibonacci(n,5)^2,n=0..16); # Zerinvary Lajos, Apr 09 2008
  • Mathematica
    LinearRecurrence[{26,26,-1},{0,1,25},30] (* Harvey P. Dale, Sep 25 2019 *)
  • SageMath
    def A099365(n): return (2/29)*(chebyshev_T(n, 27/2) - (-1)^n)
    [A099365(n) for n in (0..30)] # G. C. Greubel, Aug 21 2022

Formula

a(n) = A052918(n-1)^2, n >= 1, a(0) = 0.
a(n) = 26*a(n-1) + 26*a(n-2) - a(n-3), n >= 3; a(0)=0, a(1)=1, a(2)=25.
a(n) = 27*a(n-1) - a(n-2) - 2*(-1)^n, n >= 2; a(0)=0, a(1)=1.
a(n) = 2*(T(n, 27/2) - (-1)^n)/29 with twice the Chebyshev's T(n, x) polynomials of the first kind. 2*T(n, 27/2) = A090248(n).
G.f.: x*(1-x)/((1-27*x+x^2)*(1+x)) = x*(1-x)/(1-26*x-26*x^2+x^3).
a(n) = (1 - (-1)^n)/2 + 25*Sum_{r=1..n-1} r*a(n-r). - Michael A. Allen, Mar 30 2023
Product_{n>=2} (1 + (-1)^n/a(n)) = (5 + sqrt(29))/10 (Falcon, 2016, p. 189, eq. (3.1)). - Amiram Eldar, Dec 03 2024

A099369 Squares of A041025(n-1), n>=1, (generalized Fibonacci).

Original entry on oeis.org

0, 1, 64, 4225, 278784, 18395521, 1213825600, 80094094081, 5284996383744, 348729667233025, 23010873040995904, 1518368891038496641, 100189335935499782400, 6610977802851947141761, 436224345652293011573824
Offset: 0

Views

Author

Wolfdieter Lang, Oct 18 2004

Keywords

Comments

See the comment in A099279. This is example a=8.
a(n+1) is the number of tilings of an n-board (a board with dimensions n X 1) using half-squares (1/2 X 1 pieces, always placed so that the shorter sides are horizontal) and (1/2,1/2)-fences if there are 8 kinds of half-squares available. A (w,g)-fence is a tile composed of two w X 1 pieces separated horizontally by a gap of width g. a(n+1) also equals the number of tilings of an n-board using (1/4,1/4)-fences and (1/4,3/4)-fences if there are 8 kinds of (1/4,1/4)-fences available. - Michael A. Allen, Apr 30 2023

Crossrefs

Cf. other squares of k-metallonacci numbers (for k=1 to 10): A007598, A079291, A092936, A099279, A099365, A099366, A099367, this sequence, A099372, A099374.

Programs

  • Mathematica
    LinearRecurrence[{65,65,-1},{0,1,64},20] (* Harvey P. Dale, Oct 05 2021 *)

Formula

a(n) = A041025(n-1)^2, n >= 1, a(0)=0.
a(n) = 65*a(n-1) + 65*a(n-2) - a(n-3), n >= 3; a(0)=0, a(1)=1, a(2)=64.
a(n) = 66*a(n-1) - a(n-2) - 2*(-1)^n, n >= 2; a(0)=0, a(1)=1.
a(n) = (T(n, 33) - (-1)^n)/34 with the Chebyshev polynomials of the first kind: T(n, 33) = A099370(n).
G.f.: x*(1-x)/((1-66*x+x^2)*(1+x)) = x*(1-x)/(1-65*x-65*x^2+x^3).
a(n) = (1 - (-1)^n)/2 + 64*Sum_{r=1..n-1} r*a(n-r). - Michael A. Allen, Apr 30 2023
Product_{n>=2} (1 + (-1)^n/a(n)) = (4 + sqrt(17))/8 (Falcon, 2016, p. 189, eq. (3.1)). - Amiram Eldar, Dec 03 2024

A099374 a(n) = A041041(n-1)^2, n >= 1.

Original entry on oeis.org

0, 1, 100, 10201, 1040400, 106110601, 10822240900, 1103762461201, 112572948801600, 11481337015302001, 1170983802612002500, 119428866529408953001, 12180573402197101203600
Offset: 0

Views

Author

Wolfdieter Lang, Oct 18 2004

Keywords

Comments

See the comment in A099279. This is example a=10.
a(n+1) is the number of tilings of an n-board (a board with dimensions n X 1) using half-squares (1/2 X 1 pieces, always placed so that the shorter sides are horizontal) and (1/2,1/2)-fences if there are 10 kinds of half-squares available. A (w,g)-fence is a tile composed of two w X 1 pieces separated horizontally by a gap of width g. a(n+1) also equals the number of tilings of an n-board using (1/4,1/4)-fences and (1/4,3/4)-fences if there are 10 kinds of (1/4,1/4)-fences available. - Michael A. Allen, Mar 21 2024

Crossrefs

Cf. other squares of k-metallonacci numbers (for k=1 to 10): A007598, A079291, A092936, A099279, A099365, A099366, A099367, A099369, A099372, this sequence.

Programs

  • Mathematica
    LinearRecurrence[{101,101,-1},{0,1,100},20] (* Harvey P. Dale, Nov 10 2021 *)

Formula

a(n) = A041041(n-1)^2, n >= 1, a(0)=0.
a(n) = 101*a(n-1) + 101*a(n-2) - a(n-3), n >= 3; a(0)=0, a(1)=1, a(2)=100.
a(n) = 102*a(n-1) - a(n-2) - 2*(-1)^n, n >= 2; a(0)=0, a(1)=1.
a(n) = (T(n, 51) - (-1)^n)/52 with the Chebyshev polynomials of the first kind: T(n, 51) = (n).
G.f.: x*(1-x)/((1-102*x+x^2)*(1+x)) = x*(1-x)/(1-101*x-101*x^2+x^3).
a(n) = (1 - (-1)^n)/2 + 100*Sum_{r=1..n-1} r*a(n-r). - Michael A. Allen, Mar 21 2024
Product_{n>=2} (1 + (-1)^n/a(n)) = (5 + sqrt(26))/10 (Falcon, 2016, p. 189, eq. (3.1)). - Amiram Eldar, Dec 03 2024

A099367 a(n) = A054413(n-1)^2, n >= 1.

Original entry on oeis.org

0, 1, 49, 2500, 127449, 6497401, 331240000, 16886742601, 860892632649, 43888637522500, 2237459621014849, 114066552034234801, 5815156694124960000, 296458924848338725201, 15113590010571150025249, 770496631614280312562500
Offset: 0

Views

Author

Wolfdieter Lang, Oct 18 2004

Keywords

Comments

See the comment in A099279. This is example a=7.

Crossrefs

Cf. A054413.
Cf. other squares of k-metallonacci numbers (for k=1 to 10): A007598, A079291, A092936, A099279, A099365, A099366, this sequence, A099369, A099372, A099374.

Programs

  • Mathematica
    LinearRecurrence[{50,50,-1},{0,1,49},20] (* Harvey P. Dale, Jul 27 2023 *)

Formula

a(n) = A054413(n-1)^2, n >= 1. a(0)=0.
a(n) = 50*a(n-1) + 50*a(n-2) - a(n-3), n >= 3; a(0)=0, a(1)=1, a(2)=49.
a(n) = 51*a(n-1) - a(n-2) - 2*(-1)^n, n >= 2; a(0)=0, a(1)=1.
a(n) = 2*(T(n, 51/2) - (-1)^n)/53 with twice the Chebyshev polynomials of the first kind: 2*T(n, 51/2) = A099368(n).
G.f.: x*(1-x)/((1-51*x+x^2)*(1+x)) = x*(1-x)/(1-50*x-50*x^2+x^3).
a(n+1) = (1 + (-1)^n)/2 + 49*Sum_{k=1..n} k*a(n+1-k). - Michael A. Allen, Feb 21 2023
Product_{n>=2} (1 + (-1)^n/a(n)) = (7 + sqrt(53))/14 (Falcon, 2016, p. 189, eq. (3.1)). - Amiram Eldar, Dec 03 2024

A099366 Squares of A005668.

Original entry on oeis.org

0, 1, 36, 1369, 51984, 1974025, 74960964, 2846542609, 108093658176, 4104712468081, 155870980128900, 5918992532430121, 224765845252215696, 8535183127051766329, 324112192982714904804, 12307728150216114616225
Offset: 0

Views

Author

Wolfdieter Lang, Oct 18 2004

Keywords

Comments

See the comment in A099279. This is example a=6.
a(n+1) is the number of tilings of an n-board (a board with dimensions n X 1) using half-squares (1/2 X 1 pieces, always placed so that the shorter sides are horizontal) and (1/2,1/2)-fences if there are 6 kinds of half-squares available. A (w,g)-fence is a tile composed of two w X 1 pieces separated horizontally by a gap of width g. a(n+1) also equals the number of tilings of an n-board using (1/4,1/4)-fences and (1/4,3/4)-fences if there are 6 kinds of (1/4,1/4)-fences available. - Michael A. Allen, Apr 21 2023

Crossrefs

Cf. other squares of k-metallonacci numbers (for k=1 to 10): A007598, A079291, A092936, A099279, A099365, this sequence, A099367, A099369, A099372, A099374.

Programs

  • Maple
    with (combinat):seq(fibonacci(n,6)^2,n=0..15); # Zerinvary Lajos, Apr 09 2008
  • Mathematica
    LinearRecurrence[{37,37,-1},{0,1,36},20] (* Harvey P. Dale, Sep 23 2018 *)

Formula

a(n) = A005668(n)^2.
a(n) = 37*a(n-1) + 37*a(n-2) - a(n-3), n >= 3; a(0)=0, a(1)=1, a(2)=36.
a(n) = 38*a(n-1) - a(n-2) - 2*(-1)^n, n >= 2; a(0)=0, a(1)=1.
a(n) = (T(n, 19) - (-1)^n)/20 with the Chebyshev polynomials of the first kind: T(n, 19) = A078986(n).
G.f.: x*(1-x)/((1 - 38*x + x^2)*(1+x)) = x*(1-x)/(1 - 37*x - 37*x^2 + x^3).
a(n) = (1 - (-1)^n)/2 + 36*Sum_{r=1..n-1} r*a(n-r). - Michael A. Allen, Apr 21 2023
Product_{n>=2} (1 + (-1)^n/a(n)) = (3 + sqrt(10))/6 (Falcon, 2016, p. 189, eq. (3.1)). - Amiram Eldar, Dec 03 2024

A176019 Decimal expansion of (3+sqrt(13))/6.

Original entry on oeis.org

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

Views

Author

Klaus Brockhaus, Apr 06 2010

Keywords

Comments

Continued fraction expansion of (3+sqrt(13))/6 is A010690.
Minimal polynomial: 9*x^2 - 9*x - 1. - Amiram Eldar, Dec 03 2024

Examples

			1.10092521257733154885320354457841599104188276230754...
		

Crossrefs

Cf. A010470 (decimal expansion of sqrt(13)), A010690 (repeat 1, 9), A092936.

Programs

  • Mathematica
    RealDigits[(3 + Sqrt[13])/6, 10, 120][[1]] (* Amiram Eldar, Dec 03 2024 *)

Formula

Equals Product_{k>=2} (1 + (-1)^k/A092936(k)). - Amiram Eldar, Dec 03 2024

A094083 Numerators of ratio of sides of n-th triple of rectangles of unit area sum around a triangle.

Original entry on oeis.org

1, 1, 1, 4, 9, 64, 25, 256, 1225, 16384, 3969, 65536, 53361, 1048576, 184041, 4194304, 41409225, 1073741824, 147744025, 4294967296, 2133423721, 68719476736, 7775536041, 274877906944, 457028729521, 17592186044416, 1690195005625
Offset: 1

Views

Author

Peter J. C. Moses, Apr 30 2004

Keywords

Comments

Page 13 of the link shows the type of configuration. When n is odd, the numerators 1,1,9,25,1225,3969,.. are A038534 and (A001790)^2, and the denominators 1,4,64,256,16384,65536,.. are A056982, A038533/2, and (A046161)^2. When n is even, the numerators 1,4,64,256,16384,65536,.. are A056982, A038533/2, and (A046161)^2, and the denominators 3,27,675,3675,297675,1440747,.. are 3*(A001803)^2. The limit of a(n+1)/a(n) as n(odd) tends to infinity = Pi^2/12, A072691. The limit of a(n+2)/a(n) as n tends to infinity = 1. a(n), for large odd n, tends to 2/(Pi*n). a(n), for large even n, tends to Pi/(6*n). The expansion of 2*x*EllipticK(x)/Pi gives the odd fractions. The expansion of 1/3*x*HypergeometricPFQ({1,1,1},{3/2,3/2},x) gives the even fractions.

Examples

			a(5) = a(5-2)*((5-2)/(5-1))^2 = 1/4*(3/4)^2 = 9/64
		

Crossrefs

Programs

  • Mathematica
    a[n_]:=If[OddQ[n], ((n/2-1)!)^2/(Pi*((n/2-1/2)!)^2), Pi*((n/2-1)!)^2/(12*((n/2-1/2)!)^2)] a[n_]:=If[OddQ[n], (2^(1-n)*(n-2)!!^2)/((n-1)/2)!^2, (2^(n-2)*((n-2)/2)!^2)/(3*(n-1)!!^2)] a[n_]:=((12+Pi^2+E^(I*n*Pi)*(Pi^2-12))*((n/2-1)!)^2)/(24*Pi*((n/2-1/2)!)^2) (CoefficientList[Series[(I*x*(6+Sqrt[3]*Pi)-2*x*Sqrt[3]*Log[x+Sqrt[x^2-1]])/(6*Sqrt[x^2-1]), {x, 0, 20}], x])^2

Formula

a(n)=a(n-2)*((n-2)/(n-1))^2, a(1)=1, a(2)=1/3. a(n)=((n/2-1)!)^2/(Pi*((n/2-1/2)!)^2) for n odd. a(n)=(2^(1-n)*(n-2)!!^2)/((n-1)/2)!^2 for n odd. a(n)=Pi*((n/2-1)!)^2/(12*((n/2-1/2)!)^2) for n even. a(n)=(2^(n-2)*((n-2)/2)!^2)/(3*(n-1)!!^2) for n even.
Showing 1-10 of 11 results. Next