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 41-50 of 133 results. Next

A151548 When A160552 is regarded as a triangle with rows of lengths 1, 1, 2, 4, 8, 16, ..., this is what the rows converge to.

Original entry on oeis.org

1, 3, 5, 7, 5, 11, 17, 15, 5, 11, 17, 19, 21, 39, 49, 31, 5, 11, 17, 19, 21, 39, 49, 35, 21, 39, 53, 59, 81, 127, 129, 63, 5, 11, 17, 19, 21, 39, 49, 35, 21, 39, 53, 59, 81, 127, 129, 67, 21, 39, 53, 59, 81, 127, 133, 91, 81, 131, 165, 199, 289, 383, 321, 127, 5, 11, 17, 19, 21, 39
Offset: 0

Views

Author

David Applegate, May 18 2009

Keywords

Comments

When convolved with A151575: (1, 0, 2, -2, 6, -10, 22, -42, 86, -170, 342, ...) equals the toothpick sequence A139250: (1, 3, 7, 11, 15, 23, 35, 43, ...). - Gary W. Adamson, May 25 2009
Equals A160552: [1, 1, 3, 1, 3, 5, ...] convolved with [1, 2, 0, 0, 0, ...], equivalent to a(n) = 2*A160552(n) + A160552(n+1). - Gary W. Adamson, Jun 04 2009
Equals (1, 0, -2, 2, -2, 2, ...) convolved with the Toothpick sequence, A139250. - Gary W. Adamson, Mar 06 2012
It appears that the sums of two successive terms give A147646. - Omar E. Pol, Feb 18 2015

Examples

			From _Omar E. Pol_, Jul 24 2009: (Start)
When written as a triangle:
1;
3;
5,7;
5,11,17,15;
5,11,17,19,21,39,49,31;
5,11,17,19,21,39,49,35,21,39,53,59,81,127,129,63;
5,11,17,19,21,39,49,35,21,39,53,59,81,127,129,67,21,39,53,59,81,127,133,91,...
(End)
		

Crossrefs

Programs

  • Maple
    G := 1/(1+x) + 4*x*mul(1+x^(2^k-1)+2*x^(2^k),k=1..20); # N. J. A. Sloane, May 23 2009
    S2:=proc(n) option remember; local i,j;
    if n <= 1 then RETURN(2*n+1); fi;
    i:=floor(log(n)/log(2));
    j:=n-2^i;
    if j=0 then 5 elif j=2^i-1 then 2*n+1
    else 2*S2(j)+S2(j+1); fi;
    end; # - N. J. A. Sloane, May 22 2009
  • Mathematica
    terms = 70; CoefficientList[1/(1 + x) + 4*x*Product[1 + x^(2^k - 1) + 2*x^(2^k), {k, 1, Log[2, terms] // Ceiling}] + O[x]^terms, x] (* Jean-François Alcover, Nov 14 2017, after N. J. A. Sloane *)

Formula

a(2^k-1) = 2^(k+1)-1 for k >= 0; otherwise a(2^k) = 5 for k >= 1; otherwise a(2^i+j) = 2a(j)+a(j+1) for i >= 2, 1 <= j <= 2^i-2. - N. J. A. Sloane, May 22 2009
G.f.: 1/(1+x) + 4*x*mul(1+x^(2^k-1)+2*x^(2^k),k=1..oo). - N. J. A. Sloane, May 23 2009
a(n) = A147646(n) - a(n-1), n >= 1. - Omar E. Pol, Feb 19 2015

A151575 G.f.: (1+x)/(1+x-2*x^2).

Original entry on oeis.org

1, 0, 2, -2, 6, -10, 22, -42, 86, -170, 342, -682, 1366, -2730, 5462, -10922, 21846, -43690, 87382, -174762, 349526, -699050, 1398102, -2796202, 5592406, -11184810, 22369622, -44739242, 89478486, -178956970, 357913942, -715827882, 1431655766, -2863311530, 5726623062
Offset: 0

Views

Author

N. J. A. Sloane, May 25 2009, based on a suggestion from Gary W. Adamson

Keywords

Comments

Or, g.f. = (1+x)/((1-x)*(1-2*x)).
A signed version of A078008, which is the main entry.
[1, 0, 2, -2, 6, -10, 22, -42, 86, ...] = an operator for toothpick sequences. The sequence convolved with A151548 = toothpick sequence A139250. The sequence convolved with A151555 = toothpick sequence A153006. - Gary W. Adamson, May 25 2009

Crossrefs

Programs

  • Mathematica
    CoefficientList[Series[(1+x)/(1+x-2x^2),{x,0,40}],x] (* or *) LinearRecurrence[{-1,2},{1,0},40] (* Harvey P. Dale, May 31 2023 *)

Formula

From R. J. Mathar, Jul 08 2009: (Start)
a(n) = (2 + (-2)^n)/3 = (-1)^n*A078008(n), n>=0.
a(n) = 2*A077925(n-2), n>1. (End)
a(n) = A084247(n+1)/2. - Philippe Deléham, Sep 21 2009
G.f.: 1 + x - x*Q(0), where Q(k) = 1 + 2*x^2 - (2*k+3)*x + x*(2*k+1 - 2*x)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Oct 05 2013

A062092 a(n) = 2*a(n-1) - (-1)^n for n > 0, a(0)=2.

Original entry on oeis.org

2, 5, 9, 19, 37, 75, 149, 299, 597, 1195, 2389, 4779, 9557, 19115, 38229, 76459, 152917, 305835, 611669, 1223339, 2446677, 4893355, 9786709, 19573419, 39146837, 78293675, 156587349, 313174699, 626349397, 1252698795, 2505397589
Offset: 0

Views

Author

Amarnath Murthy, Jun 16 2001

Keywords

Comments

Let A be the Hessenberg matrix of order n, defined by: A[1,j] = A[i,i] = 1, A[i,i-1] = -1, and A[i,j] = 0 otherwise. Then, for n>=1, a(n-1) = charpoly(A,3). - Milan Janjic, Jan 24 2010

References

  • T. Koshy, Fibonacci and Lucas numbers with applications, Wiley, 2001, p. 98.

Crossrefs

Cf. A171160 (first differences).

Programs

Formula

a(n) = a(n-1) + 2*a(n-2).
a(n) = (7*2^n - (-1)^n)/3.
a(n) = 2^(n+1) + A001045(n).
A002487(a(n)) = A000032(n+1).
G.f.: (2+3*x)/(1-x-2*x^2).
E.g.f.: (7*exp(2*x) - exp(-x))/3.
a(n) = Sum_{j=0..2} A001045(n-j) (sum of 3 consecutive elements of the Jacobsthal sequence). - Alexander Adamchuk, May 16 2006
From Paul Curtz, Jun 03 2022: (Start)
a(n) = A001045(n+3) - A078008(n).
a(n) = A078008(n+3) - A001045(n).
a(n) = A005009(n-1) - a(n-1) for n >= 1.
a(n) = a(n-2) + A005009(n-2) for n >= 2.
a(n) = A154879(n-2) + 3*A201630(n-2) for n >= 2. (End)

Extensions

More terms from Jason Earls, Jun 18 2001
Additional comments from Michael Somos, Jun 24 2002

A140966 a(n) = (5 + (-2)^n)/3.

Original entry on oeis.org

2, 1, 3, -1, 7, -9, 23, -41, 87, -169, 343, -681, 1367, -2729, 5463, -10921, 21847, -43689, 87383, -174761, 349527, -699049, 1398103, -2796201, 5592407, -11184809, 22369623, -44739241, 89478487, -178956969, 357913943, -715827881, 1431655767, -2863311529, 5726623063
Offset: 0

Views

Author

Paul Curtz, Jul 27 2008

Keywords

Comments

Inverse binomial transform of A048573.
This is an example of the case k=-1 of sequences with recurrences a(n) = k*a(n-1) + (k+3)*a(n-2) - (2*k+2)*a(n-3).
The case k=1 is covered, for example, by A097163, A135520, A136326, A136336, or A137208.
Sequences with k=2 are A094554 and A094555.
Sequences with k=3 are A084175, A108924, and A139818.

Crossrefs

Programs

Formula

a(n) = -a(n-1) + 2*a(n-2).
G.f.: (2+3*x)/((1-x)*(1+2*x)).
a(n+1) - a(n) = (-1)^(n+1)*A000079(n).
a(n+3) = (-1)^n*A083582(n).
a(n+1) - 2*a(n) = -a(n+2).
a(n+1) - 3*a(n) = 5*(-1)^(n+1)*A078008(n) = (-1)^(n+1)*A001045(n-1).
a(2n+3) = -A083584(n), a(2n) = A163834(n). - Philippe Deléham, Feb 24 2014
E.g.f.: (5*exp(x) + exp(-2*x))/3. - Stefano Spezia, Jul 27 2024

Extensions

Definition simplified by R. J. Mathar, Sep 11 2009

A015552 a(n) = 6*a(n-1) + 7*a(n-2), a(0) = 0, a(1) = 1.

Original entry on oeis.org

0, 1, 6, 43, 300, 2101, 14706, 102943, 720600, 5044201, 35309406, 247165843, 1730160900, 12111126301, 84777884106, 593445188743, 4154116321200, 29078814248401, 203551699738806, 1424861898171643, 9974033287201500, 69818233010410501, 488727631072873506
Offset: 0

Views

Author

Keywords

Comments

Number of walks of length n between any two distinct nodes of the complete graph K_8. Example: a(2)=6 because the walks of length 2 between the nodes A and B of the complete graph ABCDEFGH are: ACB, ADB, AEB, AFB, AGB and AHB. - Emeric Deutsch, Apr 01 2004
The ratio a(n+1)/a(n) converges to 7 as n approaches infinity. - Felix P. Muga II, Mar 09 2014

Examples

			G.f. = x + 6*x^2 + 43*x^3 + 300*x^4 + 2101*x^5 + 14706*x^6 + 102943*x^7 + ...
		

Crossrefs

Programs

  • Magma
    [Round(7^n/8): n in [0..30]]; // Vincenzo Librandi, Jun 24 2011
  • Maple
    seq(round(7^n/8),n=0..25); # Mircea Merca, Dec 28 2010
  • Mathematica
    k=0;lst={k};Do[k=7^n-k;AppendTo[lst, k], {n, 0, 5!}];lst (* Vladimir Joseph Stephan Orlovsky, Dec 11 2008 *)
    Table[(7^n - (-1)^n)/8, {n,0,30}] (* G. C. Greubel, Dec 30 2017 *)
  • PARI
    {a(n) = if ( n<0, 0, (7^n - (-1)^n) / 8)};
    
  • Sage
    [lucas_number1(n,6,-7) for n in range(0, 21)] # Zerinvary Lajos, Apr 24 2009
    

Formula

a(n) = 6*a(n-1) + 7*a(n-2).
From Emeric Deutsch, Apr 01 2004: (Start)
G.f.: x/(1-6*x-7*x^2).
a(n) = 7^(n-1) - a(n-1). (End)
a(n) = (7^n - (-1)^n)/8. - Rolf Pleisch, Jul 06 2009
a(n) = round(7^n/8). - Mircea Merca, Dec 28 2010
E.g.f. exp(3*x)*sinh(4*x)/4. - Elmo R. Oliveira, Aug 17 2024

A015540 a(n) = 5*a(n-1) + 6*a(n-2), a(0) = 0, a(1) = 1.

Original entry on oeis.org

0, 1, 5, 31, 185, 1111, 6665, 39991, 239945, 1439671, 8638025, 51828151, 310968905, 1865813431, 11194880585, 67169283511, 403015701065, 2418094206391, 14508565238345, 87051391430071, 522308348580425, 3133850091482551, 18803100548895305, 112818603293371831
Offset: 0

Views

Author

Keywords

Comments

Number of walks of length n between any two distinct vertices of the complete graph K_7. Example: a(2)=5 because the walks of length 2 between the vertices A and B of the complete graph ABCDEFG are ACB, ADB, AEB, AFB and AGB. - Emeric Deutsch, Apr 01 2004
Pisano period lengths: 1, 1, 2, 2, 2, 2, 14, 2, 2, 2, 10, 2, 12, 14, 2, 2, 16, 2, 18, 2, ... - R. J. Mathar, Aug 10 2012
Sum_{i=0..m} (-1)^(m+i)*6^i, for m >= 0, gives all terms after 0. - Bruno Berselli, Aug 28 2013
The ratio a(n+1)/a(n) converges to 6 as n approaches infinity. Also A053524, A080424, A051958. - Felix P. Muga II, Mar 09 2014

Examples

			G.f. = x + 5*x^2 + 31*x^3 + 185*x^4 + 1111*x^5 + 6665*x^6 + 39991*x^7 + ...
		

Crossrefs

Partial sums are in A033116. Cf. A014987.

Programs

  • Magma
    [Floor(6^n/7-(-1)^n/7): n in [0..30]]; // Vincenzo Librandi, Jun 24 2011
    
  • Maple
    seq(round(6^n/7),n=0..25); # Mircea Merca, Dec 28 2010
  • Mathematica
    k=0; lst={k}; Do[k = 6^n-k; AppendTo[lst, k], {n, 0, 5!}];lst (* Vladimir Joseph Stephan Orlovsky, Dec 11 2008 *)
    CoefficientList[Series[x / ((1 - 6 x) (1 + x)), {x, 0, 50}], x] (* Vincenzo Librandi, Mar 26 2014 *)
    LinearRecurrence[{5,6},{0,1},30] (* Harvey P. Dale, May 12 2015 *)
  • PARI
    my(x='x+O('x^30)); concat([0], Vec(x/((1-6*x)*(1+x)))) \\ G. C. Greubel, Jan 24 2018
    
  • PARI
    a(n) = round(6^n/7); \\ Altug Alkan, Sep 05 2018
  • Sage
    [lucas_number1(n,5,-6) for n in range(21)] # Zerinvary Lajos, Apr 24 2009
    

Formula

a(n) = 5*a(n-1) + 6*a(n-2).
From Paul Barry, Apr 20 2003: (Start)
a(n) = (6^n - (-1)^n)/7.
G.f.: x/((1-6*x)*(1+x)).
E.g.f.: (exp(6*x) - exp(-x))/7. (End)
a(n) = 6^(n-1) - a(n-1). - Emeric Deutsch, Apr 01 2004
a(n+1) = Sum_{k=0..n} binomial(n-k, k)*5^(n-2*k)*6^k. - Paul Barry, Jul 29 2004
a(n) = round(6^n/7). - Mircea Merca, Dec 28 2010
a(n) = (-1)^(n-1)*Sum_{k=0..n-1} A135278(n-1,k)*(-7)^k = (6^n - (-1)^n)/7 = (-1)^(n-1)*Sum_{k=0..n-1} (-6)^k. Equals (-1)^(n-1)*Phi(n,-6), where Phi is the cyclotomic polynomial when n is an odd prime. (For n > 0.) - Tom Copeland, Apr 14 2014

A053088 a(n) = 3*a(n-2) + 2*a(n-3) for n > 2, a(0)=1, a(1)=0, a(2)=3.

Original entry on oeis.org

1, 0, 3, 2, 9, 12, 31, 54, 117, 224, 459, 906, 1825, 3636, 7287, 14558, 29133, 58248, 116515, 233010, 466041, 932060, 1864143, 3728262, 7456549, 14913072, 29826171, 59652314, 119304657, 238609284, 477218599, 954437166, 1908874365
Offset: 0

Views

Author

Pauline Gorman (pauline(AT)gorman65.freeserve.co.uk), Feb 26 2000

Keywords

Comments

Growth of happy bug population in GCSE math course work assignment.
The generalized (3,2)-Padovan sequence p(3,2;n). See the W. Lang link under A000931. - Wolfdieter Lang, Jun 25 2010
With offset 1: a(n) = -2^n*Sum_{k=0..n} k^p*q^k for p=1, q=-1/2. See also A232603 (p=2, q=-1/2), A232604 (p=3, q=-1/2). - Stanislav Sykora, Nov 27 2013
From Paul Curtz, Nov 02 2021 (Start)
a(n-2) difference table (from 0, 0, a(n)):
0 0 1 0 3 2 9 12 31 54 ...
0 1 -1 3 -1 7 3 19 23 63 ...
1 -2 4 -4 8 -4 16 4 40 44 ...
-3 6 -8 12 -12 20 -12 36 4 84 ...
9 -14 20 -24 32 -32 48 -32 80 0 ...
-23 34 -44 56 -64 80 -80 112 -80 176 ...
57 -78 100 -120 144 -160 192 -192 256 -192 ...
... .
The signature is valid for every row.
a(n-2) + a(n-1) = A001045(n).
a(n-2) + a(n+1) = A062510(n) = 3*A001045(n).
a(n-2) + a(n+3) = see A144472(n+1).
Second subdiagonal: 1, 6, 20, 56, 144, 352, ... = A014480(n).
First subdiagonal: -A036895(n) = -2*A001787(n).
Main diagonal: A001787(n) = -first and -third upper diagonals.
Second, fourth and fifth upper diagonals: A001792(n), A045891(n+2) and A172160(n+1). (End)

Crossrefs

Programs

  • Mathematica
    CoefficientList[Series[1/(1 - 3 x^2 - 2 x^3), {x, 0, 32}], x] (* Michael De Vlieger, Sep 30 2019 *)
  • PARI
    c(n)=(2^(n+1)-(-1)^n*(3*n+2))/9; a(n)=c(n+1); \\ Stanislav Sykora, Nov 27 2013

Formula

G.f.: 1 / (1-3*x^2-2*x^3).
With offset 1: a(1)=1; a(n) = 2*a(n-1) - (-1)^n*n; a(n) = (1/9)*(2^(n+1) - (-1)^n*(3*n+2)). - Benoit Cloitre, Nov 02 2002
a(n) = Sum_{k=0..floor(n/2)} A078008(n-2k). - Paul Barry, Nov 24 2003
a(n) = Sum_{k=0..floor(n/2)} binomial(k, n-2k)*3^k*(2/3)^(n-2k). - Paul Barry, Oct 16 2004
a(n) = Sum_{k=0..n} A078008(k)*(1 - (-1)^(n+k-1))/2. - Paul Barry, Apr 16 2005
a(n) = ( 2^(n+2) + (-1)^n*(3*n+5) )/9 (see also the B. Cloitre comment above). From the o.g.f. 1/(1-3*x^2-2*x^3) = 1/((1-2*x)*(1+x)^2) = (3/(1+x)^2 + 2/(1+x) + 4/(1-2*x))/9. - Wolfdieter Lang, Jun 25 2010
From Wolfdieter Lang, Aug 26 2010: (Start)
a(n) = a(n-1) + 2*a(n-2) + (-1)^n for n > 1, a(0)=1, a(1)=0.
Due to the identity for the o.g.f. A(x): A(x) = x*(1+2*x)*A(x) + 1/(1+x).
(This recurrence was observed by Gary Detlefs in a 08/25/10 e-mail to the author.) (End)
G.f.: Sum_{n>=0} binomial(3*n,n)*x^n / (1+x)^(3*n+3). - Paul D. Hanna, Mar 03 2012
E.g.f.: 1 + (1/9)*(exp(-x)*(3*x - 2) + 2*exp(2*x)). - Stefano Spezia, Sep 27 2019

Extensions

More terms from James Sellers, Feb 28 2000 and Christian G. Bower, Feb 29 2000

A109500 Number of closed walks of length n on the complete graph on 6 nodes from a given node.

Original entry on oeis.org

1, 0, 5, 20, 105, 520, 2605, 13020, 65105, 325520, 1627605, 8138020, 40690105, 203450520, 1017252605, 5086263020, 25431315105, 127156575520, 635782877605, 3178914388020, 15894571940105, 79472859700520
Offset: 0

Views

Author

Mitch Harris, Jun 30 2005

Keywords

Crossrefs

Cf. A109502.
Cf. sequences with the same recurrence form: A001045, A078008, A097073, A115341, A015518, A054878, A015521, A109499, A015531. - Vladimir Joseph Stephan Orlovsky, Dec 11 2008

Programs

  • Magma
    [(5^n + 5*(-1)^n)/6: n in [0..30]]; // G. C. Greubel, Dec 30 2017
  • Mathematica
    k=0;lst={k};Do[k=5^n-k;AppendTo[lst, k], {n, 1, 5!}];lst (* Vladimir Joseph Stephan Orlovsky, Dec 11 2008 *)
    CoefficientList[Series[(1 - 4*x)/(1 - 4*x - 5*x^2), {x, 0, 50}], x] (* or *) Table[(5^n + 5*(-1)^n)/6, {n,0,30}] (* G. C. Greubel, Dec 30 2017 *)
  • PARI
    for(n=0, 30, print1((5^n + 5*(-1)^n)/6, ", ")) \\ G. C. Greubel, Dec 30 2017
    

Formula

G.f.: (1 - 4*x)/(1 - 4*x - 5*x^2).
a(n) = (5^n + 5*(-1)^n)/6.
a(n) = 5^(n-1) - a(n-1), a(0) = 1. - Jon E. Schoenfield, Feb 08 2015

Extensions

Corrected by Franklin T. Adams-Watters, Sep 18 2006
Edited by Jon E. Schoenfield, Feb 08 2015

A154879 Third differences of the Jacobsthal sequence A001045.

Original entry on oeis.org

3, -2, 4, 0, 8, 8, 24, 40, 88, 168, 344, 680, 1368, 2728, 5464, 10920, 21848, 43688, 87384, 174760, 349528, 699048, 1398104, 2796200, 5592408, 11184808, 22369624, 44739240, 89478488, 178956968, 357913944, 715827880, 1431655768, 2863311528, 5726623064
Offset: 0

Views

Author

Paul Curtz, Jan 16 2009

Keywords

Comments

Second differences of A078008. First differences of the sequence (-1)^(n+1)*A084247(n).

Crossrefs

Cf. A115341.

Programs

  • Magma
    [(1/3)*(8*(-1)^n+2^n): n in [0..35]]; // Vincenzo Librandi, Jul 24 2011
    
  • Mathematica
    Differences[LinearRecurrence[{1,2},{0,1},40],3] (* or *) LinearRecurrence[ {1,2},{3,-2},40] (* Harvey P. Dale, Apr 20 2018 *)
  • Python
    def A154879(n): return ((1<2 else (3,-2,4)[n] # Chai Wah Wu, Apr 18 2025

Formula

a(n) + a(n+1) = A000079(n), n > 1.
a(n+3) = 8*A001045(n) = 4*A078008(n+1) = 2*A097073(n+1).
From R. J. Mathar, Jan 23 2009: (Start)
a(n) = a(n-1) + 2*a(n-2).
G.f.: (3-5*x)/((1+x)*(1-2*x)). (End)

Extensions

Edited and extended by R. J. Mathar, Jan 23 2009
Typo in A-number in formula corrected by R. J. Mathar, Feb 23 2009

A003665 a(n) = 2^(n-1)*( 2^n + (-1)^n ).

Original entry on oeis.org

1, 1, 10, 28, 136, 496, 2080, 8128, 32896, 130816, 524800, 2096128, 8390656, 33550336, 134225920, 536854528, 2147516416, 8589869056, 34359869440, 137438691328, 549756338176, 2199022206976, 8796095119360, 35184367894528, 140737496743936, 562949936644096, 2251799847239680
Offset: 0

Views

Author

Keywords

Comments

Binomial transform of expansion of cosh(3*x), the aerated version of A001019, 1,0,9,0,81,0,729,... - Paul Barry, Apr 05 2003
Alternatively: start with the fraction 1/1, take the numerators of fractions built according to the rule: add top and bottom to get the new bottom, add top and 9 times the bottom to get the new top. The limit of the sequence of fractions used to generate this sequence is sqrt(9). - Cino Hilliard, Sep 25 2005
This sequence also gives the number of ordered pairs of subsets (A, B) of {1, 2, ..., n} such that |A u B| is even. (Here "u" stands for the set-theoretic union.) The special case n = 13 can be found as in CRUX Problem 3257. - Walther Janous (walther.janous(AT)tirol.com), Mar 01 2008
For n > 0, a(n) is term (1,1) in the n-th power of the 2 X 2 matrix [1,3; 3,1]. - Gary W. Adamson, Aug 06 2010
a(n) is the number of compositions of n when there are 1 type of 1 and 9 types of other natural numbers. - Milan Janjic, Aug 13 2010
a(n) = ((1+3)^n+(1-3)^n)/2. In general, if b(0),b(1),... is the k-th binomial transform of the sequence ((3^n+(-3)^n)/2), then b(n) = ((k+3)^n+(k-3)^n)/2. More generally, if b(0),b(1),... is the k-th binomial transform of the sequence ((m^n+(-m)^n)/2), then b(n) = ((k+m)^n+(k-m)^n)/2. See A034494, A081340-A081342, A034659. - Charlie Marion, Jun 25 2011
Pisano period lengths: 1, 1, 1, 1, 4, 1, 6, 1, 1, 4, 5, 1, 12, 6, 4, 1, 8, 1, 9, 4, ... - R. J. Mathar, Aug 10 2012
Starting with offset 1 the sequence is the INVERT transform of (1, 9, 9, 9, ...). - Gary W. Adamson, Aug 06 2016

References

  • John Derbyshire, Prime Obsession, Joseph Henry Press, April 2004, p. 16.
  • M. Gardner, Riddles of Sphinx, M.A.A., 1987, p. 145.

Crossrefs

Programs

  • GAP
    List([0..30], n-> 2^(n-1)*(2^n +(-1)^n)); # G. C. Greubel, Aug 02 2019
  • Magma
    [2^(n-1)*( 2^n + (-1)^n ): n in [0..30]]; // Vincenzo Librandi, Aug 19 2011
    
  • Maple
    A003665:=n->2^(n-1)*( 2^n + (-1)^n ): seq(A003665(n), n=0..30); # Wesley Ivan Hurt, Apr 28 2017
  • Mathematica
    CoefficientList[Series[(1+8x)/(1-2x-8x^2), {x,0,30}], x] (* or *)
    LinearRecurrence[{2,8}, {1,1}, 30] (* Robert G. Wilson v, Sep 18 2013 *)
  • PARI
    a(n)=2^(n-1)*( 2^n + (-1)^n );
    
  • Sage
    [2^(n-1)*(2^n +(-1)^n) for n in (0..30)] # G. C. Greubel, Aug 02 2019
    

Formula

From Paul Barry, Mar 01 2003: (Start)
a(n) = 2*a(n-1) + 8*a(n-2), a(0)=a(1)=1.
a(n) = (4^n + (-2)^n)/2.
G.f.: (1-x)/((1+2*x)*(1-4*x)). (End)
From Paul Barry, Apr 05 2003: (Start)
a(n) = Sum_{k=0..floor(n/2)} binomial(n, 2*k)*9^k.
E.g.f. exp(x)*cosh(3*x). (End)
a(n) = (A078008(n) + A001045(n+1))*2^(n-1) = A014551(n)*2^(n-1). - Paul Barry, Feb 12 2004
Given a(0)=1, b(0)=1 then for i=1, 2, ..., a(i)/b(i) = (a(i-1) + 9*b(i-1)) / (a(i-1) + b(i-1)). - Cino Hilliard, Sep 25 2005
a(n) = Sum_{k=0..n} A098158(n,k)*9^(n-k). - Philippe Deléham, Dec 26 2007
a(n) = ((1+sqrt(9))^n + (1-sqrt(9))^n)/2. - Al Hakanson (hawkuu(AT)gmail.com), Dec 08 2008
If p[1]=1, and p[i]=9, (i>1), and if A is Hessenberg matrix of order n defined by: A[i,j]=p[j-i+1], (i<=j), A[i,j]=-1, (i=j+1), and A[i,j]=0 otherwise. Then, for n>=1, a(n)=det(A). - Milan Janjic, Apr 29 2010
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - x*(9*k-1)/(x*(9*k+8) - 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 28 2013

Extensions

Entry revised by N. J. A. Sloane, Nov 22 2006
Previous Showing 41-50 of 133 results. Next