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

A024495 a(n) = C(n,2) + C(n,5) + ... + C(n, 3*floor(n/3)+2).

Original entry on oeis.org

0, 0, 1, 3, 6, 11, 21, 42, 85, 171, 342, 683, 1365, 2730, 5461, 10923, 21846, 43691, 87381, 174762, 349525, 699051, 1398102, 2796203, 5592405, 11184810, 22369621, 44739243, 89478486, 178956971, 357913941, 715827882, 1431655765, 2863311531, 5726623062
Offset: 0

Views

Author

Keywords

Comments

Trisections give A082365, A132804, A132805. - Paul Curtz, Nov 18 2007
If the offset is changed to 1, this is the maximal number of closed regions bounded by straight lines after n straight line cuts in a plane: a(n) = a(n-1) + n - 3, a(1)=0; a(2)=0; a(3)=1; and so on. - Srikanth K S, Jan 23 2008
M^n * [1,0,0] = [A024493(n), a(n), A024494(n)]; where M = a 3x3 matrix [1,1,0; 0,1,1; 1,0,1]. Sum of terms = 2^n. Example: M^5 * [1,0,0] = [11, 11, 10], sum = 2^5 = 32. - Gary W. Adamson, Mar 13 2009
For n>=1, a(n-1) is the number of generalized compositions of n when there are i^2/2 - 3*i/2 + 1 different types of i, (i=1,2,...). - Milan Janjic, Sep 24 2010
Let M be any endomorphism on any vector space, such that M^3 = 1 (identity). Then (1+M)^n = A024493(n) + A024494(n)*M + a(n)*M^2. - Stanislav Sykora, Jun 10 2012
{A024493, A131708, A024495} is the difference analog of the hyperbolic functions {h_1(x), h_2(x), h_3(x)} of order 3. For the definitions of {h_i(x)} and the difference analog {H_i(n)} see [Erdelyi] and the Shevelev link respectively. - Vladimir Shevelev, Aug 01 2017
This is the p-INVERT of (1,1,1,1,1,...) for p(S) = 1 - S^3; see A291000. - Clark Kimberling, Aug 24 2017

References

  • A. Erdelyi, Higher Transcendental Functions, McGraw-Hill, 1955, Vol. 3, Chapter XVIII.
  • D. E. Knuth, The Art of Computer Programming. Addison-Wesley, Reading, MA, Vol. 1, 2nd. ed., Problem 38, p. 70.

Crossrefs

Sequences of the form 1/((1-x)^m - x^m): A000079 (m=1,2), this sequence (m=3), A000749 (m=4), A049016 (m=5), A192080 (m=6), A049017 (m=7), A290995 (m=8), A306939 (m=9).

Programs

  • Magma
    R:=PowerSeriesRing(Integers(), 30); [0,0] cat Coefficients(R!( x^2/((1-x)^3-x^3) )); // G. C. Greubel, Apr 11 2023
    
  • Maple
    a:= proc(n) option remember; `if`(n=0, 0, 2*a(n-1)+
          [-1, 0, 1, 1, 0, -1, -1][1+(n mod 6)])
        end:
    seq(a(n), n=0..33); # Paul Weisenhorn, May 17 2020
  • Mathematica
    LinearRecurrence[{3,-3,2},{0,0,1},40] (* Harvey P. Dale, Sep 20 2016 *)
  • PARI
    a(n) = sum(k=0,n\3,binomial(n,3*k+2)) /* Michael Somos, Feb 14 2006 */
    
  • PARI
    a(n)=if(n<0, 0, ([1,0,1;1,1,0;0,1,1]^n)[3,1]) /* Michael Somos, Feb 14 2006 */
    
  • SageMath
    def A024495(n): return (2^n - chebyshev_U(n, 1/2) - chebyshev_U(n-1, 1/2))/3
    [A024495(n) for n in range(41)] # G. C. Greubel, Apr 11 2023

Formula

a(n) = ( 2^n + 2*cos((n-4)*Pi/3) )/3 = (2^n - A057079(n))/3.
a(n) = 2*a(n-1) + A010892(n-2) = a(n-1) + A024494(n-1). With initial zero, binomial transform of A011655 which is effectively A010892 unsigned. - Henry Bottomley, Jun 04 2001
a(2) = 1, a(3) = 3, a(n+2) = a(n+1) - a(n) + 2^n. - Benoit Cloitre, Sep 04 2002
a(n) = Sum_{k=0..n} 2^k*2*sin(Pi*(n-k)/3 + Pi/3)/sqrt(3) (offset 0). - Paul Barry, May 18 2004
G.f.: x^2/((1-x)^3 - x^3) = x^2 / ( (1-2*x)*(1-x+x^2) ).
a(n) = 3*a(n-1) - 3*a(n-2) + 2*a(n-3). - Paul Curtz, Nov 18 2007
a(n) + A024493(n-1) = A131577(n). - Paul Curtz, Jan 24 2008
From Paul Curtz, May 29 2011: (Start)
a(n) + a(n+3) = 3*2^n = A007283(n).
a(n+6) - a(n) = 21*2^n = A175805(n).
a(n) + a(n+9) = 171*2^n.
a(n+12) - a(n) = 1365*2^n. (End)
a(n) = A113405(n) + A113405(n+1). - Paul Curtz, Jun 05 2011
Start with x(0)=1, y(0)=0, z(0)=0 and set x(n+1) = x(n) + z(n), y(n+1) = y(n) + x(n), z(n+1) = z(n) + y(n). Then a(n) = z(n). - Stanislav Sykora, Jun 10 2012
G.f.: -x^2/( x^3 - 1 + 3*x/Q(0) ) where Q(k) = 1 + k*(x+1) + 3*x - x*(k+1)*(k+4)/Q(k+1) ; (continued fraction). - Sergei N. Gladkovskii, Mar 15 2013
a(n) = 1/18*(-4*(-1)^floor((n - 1)/3) - 6*(-1)^floor(n/3) - 3*(-1)^floor((n + 1)/3) + (-1)^(1 + floor((n + 2)/3)) + 3*2^(n + 1)). - John M. Campbell, Dec 23 2016
a(n) = (1/63)*(-40 + 21*2^n - 42*floor(n/6) + 32*floor((n+3)/6) + 16*floor((n+ 4)/6) - 24*floor((n+5)/6) - 22*floor((n+7)/6) + 21*floor((n+8)/6) + 10*floor((n+9)/6) + 5*floor((n+10)/6) + 3*floor((n+11)/6) + floor((n+ 13)/6)). - John M. Campbell, Dec 24 2016
a(n+m) = a(n)*A024493(m) + A131708(n)*A131708(m) + A024493(n)*a(m). - Vladimir Shevelev, Aug 01 2017
From Kevin Ryde, Sep 24 2020: (Start)
a(n) = (1/3)*2^n - (1/3)*cos((1/3)*Pi*n) - (1/sqrt(3))*sin((1/3)*Pi*n). [Cournot]
a(n) + A111927(n) + A131708(n) = 2^n - 1. [Cournot, page 96 last formula, but misprint should be 2^x - 1 rather than 2^p - 1] (End)
E.g.f.: (exp(2*x) - exp(x/2)*(cos(sqrt(3)*x/2) + sqrt(3)*sin(sqrt(3)*x/2)))/3. - Stefano Spezia, Feb 06 2025

A153130 Period 6: repeat [1, 2, 4, 8, 7, 5].

Original entry on oeis.org

1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, 5
Offset: 0

Views

Author

Paul Curtz, Dec 19 2008

Keywords

Comments

Digital root of 2^n.
A regular version of Pitoun's sequence: a(n) = A029898(n+1).
Also obtained from permutations of A141425, A020806, A070366, A153110, A153990, A154127, A154687, or A154815.
This sequence and its (again period 6) repeated differences produce the table:
1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, ...
1, 2, 4, -1, -2, -4, 1, 2, 4, -1, -2, ...
1, 2, -5, -1, -2, 5, 1, 2, -5, -1, -2, ...
1, -7, 4, -1, 7, -4, 1, -7, 4, -1, 7, ...
-8, 11, -5, 8,-11, 5, -8, 11, -5, 8,-11, ...
19,-16, 13,-19, 16,-13, 19,-16, 13,-19, 16, ...
-35, 29,-32, 35,-29, 32,-35, 29,-32, 35,-29, ...
64,-61, 67,-64, 61,-67, 64,-61, 67,-64, 61, ...
If each entry of this table is read modulo 9 we obtain the very regular table:
1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, ...
1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, ...
1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, ...
1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, ...
1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, ...
1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, ...
Also the decimal expansion of the constant 125/1001. - R. J. Mathar, Jan 23 2009
Digital root of the powers of any number congruent to 2 mod 9. - Alonso del Arte, Jan 26 2014

References

  • Cecil Balmond, Number 9: The Search for the Sigma Code. Munich, New York: Prestel (1998): 203.

Crossrefs

Cf. digital roots of powers of c mod 9: c = 4, A100402; c = 5, A070366; c = 7, A070403; c = 8, A010689.

Programs

Formula

a(n) + a(n+3) = 9 = A010734(n).
G.f.: (1+x+2x^2+5x^3)/((1-x)(1+x)(1-x+x^2)). - R. J. Mathar, Jan 23 2009
a(n) = A082365(n) mod 9. - Paul Curtz, Mar 31 2009
a(n) = -1/2*cos(Pi*n) - 3*cos(1/3*Pi*n) - 3^(1/2)*sin(1/3*Pi*n) + 9/2. - Leonid Bedratyuk, May 13 2012
a(n) = A010888(A004000(n+1)). - Ivan N. Ianakiev, Nov 27 2014
From Wesley Ivan Hurt, Apr 20 2015: (Start)
a(n) = a(n-6) for n>5.
a(n) = a(n-1) - a(n-3) + a(n-4) for n>3.
a(n) = (2+3*(n-1 mod 3))*(n mod 2) + (1+3*(-n mod 3))*(n-1 mod 2). (End)
a(n) = 2^n mod 9. - Nikita Sadkov, Oct 06 2018
From Stefano Spezia, Mar 20 2025: (Start)
E.g.f.: 4*cosh(x) - exp(x/2)*(3*cos(sqrt(3)*x/2) + sqrt(3)*sin(sqrt(3)*x/2)) + 5*sinh(x).
a(n) = A007953(2*a(n-1)) = A010888(2*a(n-1)). (End)

Extensions

Edited by R. J. Mathar, Apr 09 2009

A015565 a(n) = 7*a(n-1) + 8*a(n-2), a(0) = 0, a(1) = 1.

Original entry on oeis.org

0, 1, 7, 57, 455, 3641, 29127, 233017, 1864135, 14913081, 119304647, 954437177, 7635497415, 61083979321, 488671834567, 3909374676537, 31274997412295, 250199979298361, 2001599834386887, 16012798675095097, 128102389400760775, 1024819115206086201, 8198552921648689607
Offset: 0

Views

Author

Keywords

Comments

A linear 2nd order recurrence. A Jacobsthal number sequence.
Binomial transform of A053573 (preceded by zero). - Paul Barry, Apr 09 2003
Second binomial transform of A080424. Binomial transform of A053573, with leading zero. Binomial transform is 0,1,9,81,729,....(9^n - 0^n)/9. Second binomial transform is 0,1,11,111,1111,... (A002275: repunits). - Paul Barry, Mar 14 2004
Number of walks of length n between any two distinct nodes of the complete graph K_9. Example: a(2)=7 because the walks of length 2 between the nodes A and B of the complete graph ABCDEFGHI are: ACB, ADB, AEB, AFB, AGB, AHB and AIB. - Emeric Deutsch, Apr 01 2004
Unsigned version of A014990. - Philippe Deléham, Feb 13 2007
The ratio a(n+1)/a(n) converges to 8 as n approaches infinity. - Felix P. Muga II, Mar 09 2014

Examples

			G.f. = x + 7*x^2 + 57*x^3 + 455*x^4 + 3641*x^5 + 29127*x^6 + 233017*x^7 + ...
		

Crossrefs

Programs

Formula

From Paul Barry, Apr 09 2003: (Start)
a(n) = (8^n - (-1)^n)/9.
a(n) = J(3*n)/3 = A001045(3*n)/3. (End)
From Emeric Deutsch, Apr 01 2004: (Start)
a(n) = 8^(n-1) - a(n-1).
G.f.: x/(1-7*x-8*x^2). (End)
a(n) = Sum_{k = 0..n} A106566(n,k)*A099322(k). - Philippe Deléham, Oct 30 2008
a(n) = round(8^n/9). - Mircea Merca, Dec 28 2010
From Peter Bala, May 31 2024: (Start)
G.f: A(x) = x/(1 - x^2) o x/(1 - x^2), where o denotes the black diamond product of power series as defined by Dukes and White. Cf. A054878.
The black diamond product A(x) o A(x) is the g.f. for the number of walks of length n between any two distinct nodes of the complete graph K_81.
Row 8 of A062160. (End)
E.g.f.: exp(-x)*(exp(9*x) - 1)/9. - Elmo R. Oliveira, Aug 17 2024

A097166 Expansion of g.f. (1+2*x)/((1-x)*(1-10*x)).

Original entry on oeis.org

1, 13, 133, 1333, 13333, 133333, 1333333, 13333333, 133333333, 1333333333, 13333333333, 133333333333, 1333333333333, 13333333333333, 133333333333333, 1333333333333333, 13333333333333333, 133333333333333333, 1333333333333333333, 13333333333333333333, 133333333333333333333
Offset: 0

Views

Author

Paul Barry, Jul 30 2004

Keywords

Comments

Partial sums of (1+2*x)/(1-10*x) = {1, 12, 120, 1200, ...}.
Second binomial transform of A082365.
These terms are the x's of A070152 and the corresponding y's are A350995 (see formula and examples). - Bernard Schott, Feb 15 2022

Examples

			a(0) = (4-1)/3 = 1 and Sum_{j=1..5} = 15.
a(1) = (40-1)/3 = 13 and Sum_{j=13..53} = 1353.
a(2) = (400-1)/3 = 133 and Sum_{j=133..533} = 133533.
		

Crossrefs

Cf. A056698 (index of primes), A082365, A097169, A309907 (squares of this).

Programs

  • Magma
    [(4*10^n-1)/3 : n in [0..20]]; // Vincenzo Librandi, Nov 01 2011
    
  • Maple
    a:= n-> parse(cat(1, 3$n)):
    seq(a(n), n=0..18);  # Alois P. Heinz, Aug 23 2019
  • Mathematica
    NestList[10#+3&,1,20] (* Harvey P. Dale, Jan 22 2014 *)
  • Python
    [(4*10**n-1)//3 for n in range(25)] # Gennady Eremin, Mar 04 2022

Formula

a(n) = (4*10^n - 1)/3.
a(n) = A097169(2*n).
a(n) = 10*a(n-1) + 3, n>0. a(n) = 11*a(n-1) - 10*a(n-2), n>1. - Vincenzo Librandi, Nov 01 2011
A350994(n) = Sum_{j=a(n)..A350995(n)} = a(n).A350995(n) where "." means concatenation. - Bernard Schott, Jan 28 2022
From Elmo R. Oliveira, Apr 29 2025: (Start)
E.g.f.: exp(x)*(4*exp(9*x) - 1)/3.
a(n) = A198970(n)/3. (End)

A082311 A Jacobsthal sequence trisection.

Original entry on oeis.org

1, 5, 43, 341, 2731, 21845, 174763, 1398101, 11184811, 89478485, 715827883, 5726623061, 45812984491, 366503875925, 2932031007403, 23456248059221, 187649984473771, 1501199875790165, 12009599006321323, 96076792050570581, 768614336404564651, 6148914691236517205
Offset: 0

Views

Author

Paul Barry, Apr 09 2003

Keywords

Crossrefs

Programs

  • Magma
    [2*8^n/3+(-1)^n/3 : n in [0..30]]; // Vincenzo Librandi, Aug 13 2011
    
  • Mathematica
    f[n_] := (2*8^n + (-1)^n)/3; Array[f, 25, 0] (* Robert G. Wilson v, Aug 13 2011 *)
  • PARI
    x='x+O('x^30); Vec((1-2*x)/((1+x)*(1-8*x))) \\ G. C. Greubel, Sep 16 2018

Formula

a(n) = (2*8^n + (-1)^n)/3 = A001045(3*n+1).
From R. J. Mathar, Feb 23 2009: (Start)
a(n) = 7*a(n-1) + 8*a(n-2).
G.f.: (1-2*x)/((1+x)*(1-8*x)). (End)
a(n) = A024494(3*n+1). a(n) = 8*a(n-1) + 3*(-1)^n. Sum of digits = A070366. - Paul Curtz, Nov 20 2007
a(n)= A007613(n) + A132805(n) = A081374(1+3*n). - Paul Curtz, Jun 06 2011
E.g.f.: (cosh(x) + 2*cosh(8*x) - sinh(x) + 2*sinh(8*x))/3. - Stefano Spezia, Jul 15 2024

A083322 a(n) = 2^n - A081374(n).

Original entry on oeis.org

1, 2, 6, 11, 22, 42, 85, 170, 342, 683, 1366, 2730, 5461, 10922, 21846, 43691, 87382, 174762, 349525, 699050, 1398102, 2796203, 5592406, 11184810, 22369621, 44739242, 89478486, 178956971, 357913942, 715827882, 1431655765, 2863311530, 5726623062
Offset: 1

Views

Author

David Applegate, Aug 22 2003

Keywords

Crossrefs

Cf. A081374.
Trisections: A082365, A007613, A132804.

Programs

  • Magma
    I:=[1,2,6,11]; [n le 4 select I[n] else 2*Self(n-1)-Self(n-3)+2*Self(n-4): n in [1..40]]; // Vincenzo Librandi, Jul 08 2016
  • Mathematica
    CoefficientList[Series[(1 + 2 x^2) / ((1 - 2 x) (1 + x) (1 - x + x^2)), {x, 0, 40}], x] (* Vincenzo Librandi, Jul 08 2016 *)
    LinearRecurrence[{2,0,-1,2},{1,2,6,11},40] (* Harvey P. Dale, Jan 30 2024 *)

Formula

G.f.: x*(1+2*x^2) / ( (1-2*x)*(1+x)*(1-x+x^2) ). - R. J. Mathar, May 27 2011
From Paul Curtz, May 27 2011: (Start)
a(n) = 2*a(n-1) - a(n-3) + 2*a(n-4).
a(n)+a(n+3) = 3*2^(n+1) = A007283(n+1).
a(n+6)-a(n) = 21*2^(n+1) = A175805(n+1).
(End)

A087462 Generalized mod 3 multiplicative Jacobsthal sequence.

Original entry on oeis.org

1, 1, 1, 8, 5, 11, 64, 43, 85, 512, 341, 683, 4096, 2731, 5461, 32768, 21845, 43691, 262144, 174763, 349525, 2097152, 1398101, 2796203, 16777216, 11184811, 22369621, 134217728, 89478485, 178956971, 1073741824, 715827883, 1431655765, 8589934592, 5726623061
Offset: 0

Views

Author

Paul Barry, Sep 08 2003

Keywords

Comments

2^n = a(n) + A087463(n) + A087464(n) provides a decomposition of Pascal's triangle.
Multiplicative analog of A078008.

Crossrefs

Cf. A001045, A001018 (trisection), A082311 (trisection), A082365 (trisection).

Programs

  • PARI
    Vec(-(4*x^5-2*x^4+x^3+x^2+x+1)/((x+1)*(2*x-1)*(x^2-x+1)*(4*x^2+2*x+1)) + O(x^100)) \\ Colin Barker, Nov 02 2015

Formula

a(n) = Sum_{k=0..n} if (mod(n*k, 3)=0, 1, 0) * C(n, k).
a(n) = 2^n-2/3*(1-cos(2*Pi*n/3))*(A001045(n)+2*A001045(n-1)+0^n).
From Colin Barker, Nov 02 2015: (Start)
a(n) = 7*a(n-3)+8*a(n-6) for n>5.
G.f.: -(4*x^5-2*x^4+x^3+x^2+x+1) / ((x+1)*(2*x-1)*(x^2-x+1)*(4*x^2+2*x+1)).
(End)

A090409 a(n) = (7*8^n + 2*(-1)^n)/9.

Original entry on oeis.org

1, 6, 50, 398, 3186, 25486, 203890, 1631118, 13048946, 104391566, 835132530, 6681060238, 53448481906, 427587855246, 3420702841970, 27365622735758, 218924981886066, 1751399855088526, 14011198840708210, 112089590725665678, 896716725805325426, 7173733806442603406
Offset: 0

Views

Author

Paul Barry, Nov 29 2003

Keywords

Crossrefs

First differences of A015565.

Programs

  • Mathematica
    LinearRecurrence[{7,8},{1,6},20] (* Harvey P. Dale, Aug 15 2016 *)

Formula

a(n) = Sum_{j=0..2} Sum_{k=0..n} C(3*n+j, 3*k)/3.
a(n) = (A007613(n) + A082311(n) + A082365(n))/3.
G.f.: (-1+x)/((1+x)*(8*x-1)). - R. J. Mathar, Dec 10 2014
From Elmo R. Oliveira, Aug 18 2024: (Start)
E.g.f.: exp(-x)*(7*exp(9*x) + 2)/9.
a(n) = 7*a(n-1) + 8*a(n-2) for n > 1. (End)

Extensions

a(20)-a(21) from Elmo R. Oliveira, Aug 18 2024

A141355 The Jacobsthal sequence, dropping each third term.

Original entry on oeis.org

1, 1, 5, 11, 43, 85, 341, 683, 2731, 5461, 21845, 43691, 174763, 349525, 1398101, 2796203, 11184811, 22369621, 89478485, 178956971, 715827883, 1431655765, 5726623061, 11453246123, 45812984491, 91625968981, 366503875925, 733007751851
Offset: 0

Views

Author

Paul Curtz, Aug 03 2008

Keywords

Comments

A001045 after removal of the subsequence A132805.

Programs

  • Python
    def A141355(n): return ((1<<(n+1<<1)-(n+1>>1)-1)|1)//3 # Chai Wah Wu, Apr 19 2025

Formula

a(2n+1)-a(2n) = 6*A015565(n).
a(4n+1)=2a(4n)-1. a(4n+2)=4a(4n+1)+1. a(4n+3)=2a(4n+2)+1. a(4n+4)=4a(4n+3)-1.
a(2n)= A082311(n). a(2n+1) = A082365(n). - R. J. Mathar, Feb 23 2009
a(n)=7*a(n-2)+8*a(n-4). G.f.: (1+x-2*x^2+4*x^3)/((1-8*x^2)*(1+x^2)). - R. J. Mathar, Feb 23 2009

Extensions

Edited and extended by R. J. Mathar, Feb 23 2009

A352692 a(n) + a(n+1) = 2^n for n >= 0 with a(0) = 4.

Original entry on oeis.org

4, -3, 5, -1, 9, 7, 25, 39, 89, 167, 345, 679, 1369, 2727, 5465, 10919, 21849, 43687, 87385, 174759, 349529, 699047, 1398105, 2796199, 5592409, 11184807, 22369625, 44739239, 89478489, 178956967, 357913945, 715827879, 1431655769, 2863311527, 5726623065, 11453246119, 22906492249
Offset: 0

Views

Author

Paul Curtz, Mar 29 2022

Keywords

Comments

Difference table D(n,k) = D(n-1,k+1) - D(n-1,k), D(0,k) = a(k):
4, -3, 5, -1, 9, 7, 25, ...
-7, 8, -6, 10, -2, 18, 14, 50, ...
15, -14, 16, -12, 20, -4, 36, 28, 100, ...
-29, 30, -28, 32, -24, 40, -8, 72, 56, 200, ...
59, -58, 60, -56, 64, -48, 80, -16, 144, 112, 400, ...
...
The diagonals are given by D(n,n+k) = a(k)*2^n.
D(n,1) = -(-1)^n* A340627(n).
a(n) - a(n) = 0, 0, 0, 0, 0, ... (trivially)
a(n+1) + a(n) = 1, 2, 4, 8, 16, ... = 2^n (by definition)
a(n+2) - a(n) = 1, 2, 4, 8, 16, ... = 2^n
a(n+3) + a(n) = 3, 6, 12, 24, 48, ... = 2^n*3
a(n+4) - a(n) = 5, 10, 20, 40, 80, ... = 2^n*5
a(n+5) + a(n) = 11, 22, 44, 88, 176, ... = 2^n*11
(...)
This table is given by T(r,n) = A001045(r)*2^n with r, n >= 0.
Sums of antidiagonals are A045883(n).
Main diagonal: A192382(n).
First upper diagonal: A054881(n+1).
First subdiagonal: A003683(n+1).
Second subdiagonal: A246036(n).
Now consider the array from c(n) = (-1)^n*a(n) with its difference table:
4, 3, 5, 1, 9, -7, 25, -39, ... = c(n)
-1, 2, -4, 8, -16, 32, -64, 128, ... = -A122803(n)
3, -6, 12, -24, 48, -96, 192, -384, ... =
-9, 18, -36, 72, -144, 288, -576, 1152, ...
27, -54, 108, -216, 432, -864, 1728, -3456, ...
...
The first subdiagonal is -A000400(n). The second is A169604(n).

Crossrefs

If a(0) = k then A001045 (k=0), A078008 (k=1), A140966 (k=2), A154879 (k=3), this sequence (k=4).
Essentially the same as A115335.

Programs

  • Maple
    a := proc(n) option remember; ifelse(n = 0, 4, 2^(n-1) - a(n-1)) end: # Peter Luschny, Mar 29 2022
    A352691 := proc(n)
        (11*(-1)^n + 2^n)/3
    end proc: # R. J. Mathar, Apr 26 2022
  • Mathematica
    LinearRecurrence[{1, 2}, {4, -3}, 40] (* Amiram Eldar, Mar 29 2022 *)
  • PARI
    a(n) = (11*(-1)^n + 2^n)/3; \\ Thomas Scheuerle, Mar 29 2022

Formula

abs(a(n)) = A115335(n-1) for n >= 1.
a(3*n) - (-1)^n*4 = A132805(n).
a(3*n+1) + (-1)^n*4 = A082311(n).
a(3*n+2) - (-1)^n*4 = A082365(n).
From Thomas Scheuerle, Mar 29 2022: (Start)
G.f.: (-4 + 7*x)/(-1 + x + 2*x^2).
Sum_{k=0..n} (-1)^(n-k)*binomial(n, k)*a(m + 2*n-k) = a(m)*2^n.
Sum_{k=0..n} (-1)^(n-k)*binomial(n, k)*a(1 + n-k) = -(-1)^n*A340627(n).
a(n) = (11*(-1)^n + 2^n)/3.
a(n + 2*m) = a(n) + A002450(m)*2^n.
a(2*n) = A192382(n+1) + (-1)^n*a(n).
a(n) = ( A045883(n) - Sum_{k=0..n-1}(-1)^k*a(k) )/n, for n > 0. (End)
a(n) = A001045(n) + 4*(-1)^n.
a(n+1) = 2*a(n) -11*(-1)^n.
a(n+2) = a(n) + 2^n.
a(n+4) = a(n) + A020714(n).
a(n+6) = a(n) + A175805(n).
a(2*n) = A163868(n).
a(2*n+1) = (2^(2*n+1) - 11)/3.

Extensions

Warning: The DATA is correct, but there may be errors in the COMMENTS, which should be rechecked. - Editors of OEIS, Apr 26 2022
Edited by M. F. Hasler, Apr 26 2022.
Showing 1-10 of 11 results. Next