A057079 -id:A057079 - cp's OEIS frontend

A024494 a(n) = C(n,1) + C(n,4) + ... + C(n, 3*floor(n/3) + 1).

1, 2, 3, 5, 10, 21, 43, 86, 171, 341, 682, 1365, 2731, 5462, 10923, 21845, 43690, 87381, 174763, 349526, 699051, 1398101, 2796202, 5592405, 11184811, 22369622, 44739243, 89478485, 178956970, 357913941, 715827883, 1431655766, 2863311531, 5726623061, 11453246122

Offset: 1

Clark Kimberling

M^n * [1,0,0] = [A024493(n), A024495(n), a(n)], where M is a 3 X 3 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

Let M be any endomorphism on any vector space such that M^3 = 1 (identity). Then (1+M)^n = A024493(n) + a(n)*M + A024495(n)*M^2. - Stanislav Sykora, Jun 10 2012

D. E. Knuth, The Art of Computer Programming. Addison-Wesley, Reading, MA, Vol. 1, 2nd. ed., Problem 38, p. 70.

G. C. Greubel, Table of n, a(n) for n = 1..1000
Index entries for linear recurrences with constant coefficients, signature (3,-3,2).

Cf. A010892, A024493, A024495, A057079.

See A131708 for another version.

[n le 3 select n else 3*Self(n-1) -3*Self(n-2) +2*Self(n-3): n in [1..40]]; // G. C. Greubel, Jan 23 2023

nn=20;a=1/(1-x);Drop[CoefficientList[Series[a x /(1-x-x^3 a^2),{x,0,nn}],x],1] (* Geoffrey Critzer, Dec 22 2013 *)
LinearRecurrence[{3,-3,2}, {1,2,3}, 40] (* G. C. Greubel, Jan 23 2023 *)

a(n) = sum(k=0,n\3,binomial(n,3*k+1)) /* Michael Somos, Feb 14 2006 */

a(n)=if(n<0, 0, ([1,0,1;1,1,0;0,1,1]^n)[2,1]) /* Michael Somos, Feb 14 2006 */

def A024494(n): return (1/3)*(2^n -chebyshev_U(n,1/2) +2*chebyshev_U(n-1,1/2))
[A024494(n) for n in range(1,41)] # G. C. Greubel, Jan 23 2023

3*a(n) = 2^n + 2*cos( (n-2)*Pi/3 ) = 2^n - A057079(n+2).
G.f.: x*(1-x)/((1-2*x)*(1-x+x^2)). - Paul Barry, Feb 11 2004
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*(1-x^2)*(1-x^3)/(1-x^6))/(1-2*x). - Michael Somos, Feb 14 2006
a(n+1) - 2*a(n) = A010892(n+1). - Michael Somos, Feb 14 2006
a(n) = 3*a(n-1) - 3*a(n-2) + 2*a(n-3). - Paul Curtz, Nov 20 2007
Equals binomial transform of (1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, ...). - Gary W. Adamson, Jul 03 2008
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) = y(n). - Stanislav Sykora, Jun 10 2012
E.g.f.: exp(x/2)*(exp(3*x/2) - cos(sqrt(3)*x/2) + sqrt(3)*sin(sqrt(3)*x/2))/3. - Stefano Spezia, Feb 06 2025

A294099 Rectangular array read by (upward) antidiagonals: A(n,k) = Sum_{j=0..k} (-1)^floor(j/2)binomial(k-floor((j+1)/2), floor(j/2))n^(k-j), n >= 1, k >= 0.

Original entry on oeis.org

1, 1, 2, 1, 3, 1, 1, 4, 5, -1, 1, 5, 11, 7, -2, 1, 6, 19, 29, 9, -1, 1, 7, 29, 71, 76, 11, 1, 1, 8, 41, 139, 265, 199, 13, 2, 1, 9, 55, 239, 666, 989, 521, 15, 1, 1, 10, 71, 377, 1393, 3191, 3691, 1364, 17, -1, 1, 11, 89, 559, 2584, 8119, 15289, 13775, 3571, 19, -2

Offset: 1

L. Edson Jeffery, Bob Selcoe and Andrew Hone, Oct 22 2017

This array is used to compute A269254: A269254(n) = least k such that A(n,k) is a prime, or -1 if no such k exists.
For detailed theory, see [Hone]. - L. Edson Jeffery, Feb 09 2018
The array can be extended to k<0 with A(n, k) = -A(n, -k-1) for all k in Z. - Michael Somos, Jun 19 2023

			Array begins:
  1   2    1    -1     -2      -1        1         2          1          -1
  1   3    5     7      9      11       13        15         17          19
  1   4   11    29     76     199      521      1364       3571        9349
  1   5   19    71    265     989     3691     13775      51409      191861
  1   6   29   139    666    3191    15289     73254     350981     1681651
  1   7   41   239   1393    8119    47321    275807    1607521     9369319
  1   8   55   377   2584   17711   121393    832040    5702887    39088169
  1   9   71   559   4401   34649   272791   2147679   16908641   133121449
  1  10   89   791   7030   62479   555281   4935050   43860169   389806471
  1  11  109  1079  10681  105731  1046629  10360559  102558961  1015229051

Robert Price, Table of n, a(n) for n = 1..5050
Andrew N. W. Hone, et al., On a family of sequences related to Chebyshev polynomials, arXiv:1802.01793 [math.NT], 2018.

Cf. A285992, A299107, A299109, A088165, A117522, A299100, A299101, A113501, A269253, A269254, A294099, A298675, A298677, A298878, A299045, A299071.
Rows: A057079, A005408, A002878, A001834, A030221, A002315, A033890, A057080, A057081, A054320, ...
Columns: A000012, A000027, A028387, ...

(* Array: *)
Grid[Table[LinearRecurrence[{n, -1}, {1, 1 + n}, 10], {n, 10}]]
(* Array antidiagonals flattened (gives this sequence): *)
A294099[n_, k_] := Sum[(-1)^(Floor[j/2]) Binomial[k - Floor[(j + 1)/2], Floor[j/2]] n^(k - j), {j, 0, k}]; Flatten[Table[A294099[n - k, k], {n, 11}, {k, 0, n - 1}]]

{A(n, k) = sum(j=0, k, (-1)^(j\2)*binomial(k-(j+1)\2, j\2)*n^(k-j))}; /* Michael Somos, Jun 19 2023 */

A(n,0) = 1, A(n,1) = n + 1, A(n,k) = n*A(n,k-1) - A(n,k-2), n >= 1, k >= 2.
G.f. for row n: (1 + x)/(1 - n*x + x^2), n >= 1.
A(n, k) = B(-n, k) where B = A299045. - Michael Somos, Jun 19 2023

A174709 Partial sums of floor(n/6).

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 1, 2, 3, 4, 5, 6, 8, 10, 12, 14, 16, 18, 21, 24, 27, 30, 33, 36, 40, 44, 48, 52, 56, 60, 65, 70, 75, 80, 85, 90, 96, 102, 108, 114, 120, 126, 133, 140, 147, 154, 161, 168, 176, 184, 192

Offset: 0

Mircea Merca, Nov 30 2010

Partial sums of A152467.

			a(7) = floor(0/6) + floor(1/6) + floor(2/6) + floor(3/6) + floor(4/6) + floor(5/6) + floor(6/6) + floor(7/6) = 0 + 0 + 0 + 0 + 0 + 0 + 1 + 1 = 2.

Vincenzo Librandi, Table of n, a(n) for n = 0..10000
Mircea Merca, Inequalities and Identities Involving Sums of Integer Functions, J. Integer Sequences, Vol. 14 (2011), Article 11.9.1.
Index entries for linear recurrences with constant coefficients, signature (2,-1,0,0,0,1,-2,1).

Cf. A008724, A152467.

[Round(n*(n-4)/12): n in [0..60]]; // Vincenzo Librandi, Jun 22 2011

Maple
```
a(n):=round(n*(n-4)/12)
```

Table[Sum[Floor[k/6], {k,0,n}], {n,0,25}] (* G. C. Greubel, Dec 13 2016 *)

a(n)=(n-2)^2\12 \\ Charles R Greathouse IV, Oct 07 2015

a(n) = round(n*(n-4)/12) = round((2*n^2 - 8*n - 1)/24).
a(n) = floor((n-2)^2/12).
a(n) = ceiling((n+1)*(n-5)/12).
a(n) = a(n-6) + n - 5, n > 5.
From R. J. Mathar, Nov 30 2010: (Start)
a(n) = 2*a(n-1) - a(n-2) + a(n-6) - 2*a(n-7) + a(n-8).
G.f.: -x^6 / ( (1+x)*(x^2-x+1)*(1+x+x^2)*(x-1)^3 ).
a(n) = -n/3 + 5/72 + n^2/12 + (-1)^n/24 + A057079(n+5)/6 + A061347(n)/18. (End)
a(6n) = A000567(n), a(6n+1) = 2*A000326(n), a(6n+2) = A033428(n), a(6n+3) = A049451(n), a(6n+4) = A045944(n), a(6n+5) = A028896(n). - Philippe Deléham, Mar 26 2013
a(n) = A008724(n-2). - R. J. Mathar, Jul 10 2015
Sum_{n>=6} 1/a(n) = Pi^2/18 - Pi/(2*sqrt(3)) + 49/12. - Amiram Eldar, Aug 13 2022

A087204 Period 6: repeat [2, 1, -1, -2, -1, 1].

Original entry on oeis.org

2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1

Offset: 0

Nikolay V. Kosinov (kosinov(AT)unitron.com.ua), Oct 19 2003

Satisfies (a(n))^2 = a(2n) + 2. Shifted differences of itself.
Moebius transform is length 6 sequence [1, -2, -3, 0, 0, 6]. - Michael Somos, Oct 22 2006
Twice the real part of x^n, where x is either of the primitive 6th roots of unity. For the root with positive imaginary part, the imaginary part of x^n is i*A128834(n)*sqrt(3)/2. - Peter Munn, Apr 25 2022

			a(2) = -1 = a(1) - a(0) = 1 - 2 = ((1+sqrt(-3))/2)^2 + ((1-sqrt(-3))/2)^2 = -1 = -2/4 + 2*sqrt(-3)/4 - 2/4 -2 sqrt(-3)/4 = -1.
G.f. = 2 + x - x^2 - 2*x^3 - x^4 + x^5 + 2*x^6 + x^7 - x^8 - 2*x^9 - x^10 + ...

A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 176.
Paulo Ribenboim, My Numbers, My Friends: Popular Lectures on Number Theory, Springer-Verlag, NY, 2000, p. 6.

Tanya Khovanova, Recursive Sequences.
Wikipedia, Lucas sequence.
Index entries for linear recurrences with constant coefficients, signature (1,-1).
Index entries for recurrences a(n) = k*a(n - 1) +/- a(n - 2).

Essentially the same as A057079 and A100051. Pairwise sums of A010892.

Cf. A128834.

&cat[[2, 1, -1, -2, -1, 1]^^20]; // Wesley Ivan Hurt, Jun 19 2016

A087204:=n->[2, 1, -1, -2, -1, 1][(n mod 6)+1]: seq(A087204(n), n=0..100); # Wesley Ivan Hurt, Jun 19 2016

PadLeft[{}, 108, {2,1,-1,-2,-1,1}] (* Harvey P. Dale, Sep 11 2011 *)
a[ n_] := {1, -1, -2, -1, 1, 2}[[Mod[n, 6, 1]]]; (* Michael Somos, Jan 29 2015 *)
a[ n_] := 2 Re[ Exp[ Pi I n / 3]]; (* Michael Somos, Mar 29 2015 *)

{a(n) = [2, 1, -1, -2, -1, 1][n%6 + 1]}; /* Michael Somos, Oct 22 2006 */

A087204(n) = if(0==n, 2, my(f = factor(n)); prod(k=1, #f~, if(f[k, 1]<=3, 1-f[k, 1], 1))); \\ (After David W. Wilson's multiplicative formula) - Antti Karttunen, Apr 28 2022

[lucas_number2(n,1,1) for n in range(0, 102)] # Zerinvary Lajos, Apr 30 2009

a(n) = a(n-1) - a(n-2), starting with a(0) = 2 and a(1) = 1.
G.f.: (2-x)/(1-x+x^2).
a(n) = Sum_{k>=0} (-1)^k*n/(n-k)*C(n-k, k).
a(n) = (1/2)*((-1)^floor(n/3) + 2*(-1)^floor((n+1)/3) + (-1)^floor((n+2)/3)).
Multiplicative with a(2^e) = -1, a(3^e) = -2, a(p^e) = 1 otherwise. - David W. Wilson, Jun 12 2005
a(n) = a(-n) = -a(n-3) for all n in Z. - Michael Somos, Oct 22 2006
E.g.f.: 2*exp(x/2)*cos(sqrt(3)*x/2). - Sergei N. Gladkovskii, Aug 12 2012
a(n) = r^n + s^n, with r=(1+i*sqrt(3))/2 and s=(1-i*sqrt(3))/2 the roots of 1-x+x^2. - Ralf Stephan, Jul 19 2013
a(n) = 2*cos(n*Pi/3). - Wesley Ivan Hurt, Jun 19 2016
Dirichlet g.f.: zeta(s)*(1-2^(1-s)-3^(1-s)+6^(1-s)). - Amiram Eldar, Jan 01 2023

Edited by Ralf Stephan, Feb 04 2005

A100051 A Chebyshev transform of 1,1,1,...

Original entry on oeis.org

1, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1

Offset: 0

Paul Barry, Oct 31 2004

1, followed by period 6: repeat [1, -1, -2, -1, 1, 2]. - Joerg Arndt, Aug 28 2024
A Chebyshev transform of 1/(1-x): if A(x) is the g.f. of a sequence, map it to ((1-x^2)/(1+x^2))A(x/(1+x^2)).
Transform of 1/(1+x) under the mapping g(x)->((1+x)/(1-x))g(x/(1-x)^2). - Paul Barry, Dec 01 2004
Multiplicative with a(p^e) = -1 if p = 2; -2 if p = 3; 1 otherwise. - David W. Wilson, Jun 10 2005

			G.f. = 1 + x - x^2 - 2*x^3 - x^4 + x^5 + 2*x^6 + x^7 - x^8 - 2*x^9 - x^10 + ...

G. C. Greubel, Table of n, a(n) for n = 0..1000
Michael Somos, Rational Function Multiplicative Coefficients
Index entries for linear recurrences with constant coefficients, signature (1,-1).

Cf. A011655, A057079, A087204, A099837, A099443, A100047, A100048, A100050.

Row sums of array A127677.

CoefficientList[Series[(1 - x^2)/(1 - x + x^2), {x,0,50}], x] (* G. C. Greubel, May 03 2017 *)
LinearRecurrence[{1,-1},{1,1,-1},80] (* Harvey P. Dale, Mar 25 2019 *)

{a(n) = - (n == 0) + [2, 1, -1, -2, -1, 1][n%6 + 1]}; /* Michael Somos, Mar 21 2011 */

From Paul Barry, Dec 01 2004: (Start)
G.f.: (1-x^2)/(1-x+x^2).
a(n) = a(n-1) - a(n-2), n>2.
a(n) = n*Sum_{k=0..floor(n/2)} (-1)^k*binomial(n-k, k)/(n-k).
a(n) = Sum_{k=0..n} binomial(n+k, 2k)*(2n/(n+k))*(-1)^k, n>1. (End)
Moebius transform is length 6 sequence [1, -2, -3, 0, 0, 6].
Euler transform of length 6 sequence [1, -2, -1, 0, 0, 1].
a(n) = a(-n). a(n) = c_6(n) if n>1 where c_k(n) is Ramanujan's sum. - Michael Somos, Mar 21 2011
a(n) = A087204(n), n>0. - R. J. Mathar, Sep 02 2008
a(n) = A057079(n+1), n>0. Dirichlet g.f. zeta(s) *(1-2^(1-s)-3^(1-s)+6^(1-s)). - R. J. Mathar, Apr 11 2011

A135356 Triangle T(n,k) read by rows: coefficients in the recurrence of sequences which equal their n-th differences.

Original entry on oeis.org

2, 2, 0, 3, -3, 2, 4, -6, 4, 0, 5, -10, 10, -5, 2, 6, -15, 20, -15, 6, 0, 7, -21, 35, -35, 21, -7, 2, 8, -28, 56, -70, 56, -28, 8, 0, 9, -36, 84, -126, 126, -84, 36, -9, 2, 10, -45, 120, -210, 252, -210, 120, -45, 10, 0, 11, -55, 165, -330, 462, -462, 330, -165, 55, -11, 2

Offset: 1

Paul Curtz, Dec 08 2007, Mar 25 2008, Apr 28 2008

Sequences which equal their p-th differences obey recurrences a(n) = Sum_{s=1..p} T(p,s)*a(n-s).
This defines T(p,s) as essentially a signed version of a chopped Pascal triangle A014410, see A130785.
For cases like p=2, 4, 6, 8, 10, 12, 14, the denominator of the rational generating function of a(n) contains a factor 1-x; depending on the first terms in the sequences a(n), additional, simpler recurrences may exist if this cancels with a factor in the numerator. - R. J. Mathar, Jun 10 2008

			Triangle begins with row n=1:
  2;
  2,   0;
  3,  -3,  2;
  4,  -6,  4,    0;
  5, -10, 10,   -5,   2;
  6, -15, 20,  -15,   6,   0;
  7, -21, 35,  -35,  21,  -7,  2;
  8, -28, 56,  -70,  56, -28,  8,  0;
  9, -36, 84, -126, 126, -84, 36, -9, 2;

Alois P. Heinz, Rows n = 1..141, flattened

Related sequences: A000079 (n=1), A131577 (n=2), (A131708 , A130785, A131562, A057079) (n=3), (A000749, A038503, A009545, A038505) (n=4), A133476 (n=5), A140343 (n=6), A140342 (n=7).

Cf. A000225, A000984, A051049, A130785.

A135356:= func< n,k | k eq n select 1-(-1)^n else (-1)^(k+1)*Binomial(n,k) >;
[A135356(n,k): k in [1..n], n in [1..12]]; // G. C. Greubel, Apr 09 2023

T:= (p, s)->  `if`(p=s, 2*irem(p, 2), (-1)^(s+1) *binomial(p, s)):
seq(seq(T(p, s), s=1..p), p=1..11);  # Alois P. Heinz, Aug 26 2011

T[p_, s_]:= If[p==s, 2*Mod[s, 2], (-1)^(s+1)*Binomial[p, s]];
Table[T[p, s], {p, 12}, {s, p}]//Flatten (* Jean-François Alcover, Feb 19 2015, after Alois P. Heinz *)

def A135356(n,k):
    if (k==n): return 2*(n%2)
    else: return (-1)^(k+1)*binomial(n,k)
flatten([[A135356(n,k) for k in range(1,n+1)] for n in range(1,13)]) # G. C. Greubel, Apr 09 2023

T(n,k) = (-1)^(k+1)*A007318(n, k). T(n,n) = 1 - (-1)^n.
Sum_{k=1..n} T(n, k) = 2.
From G. C. Greubel, Apr 09 2023: (Start)
Sum_{k=1..n} (-1)^(k-1)*T(n, k) = 2*A051049(n-1).
Sum_{k=1..n-1} T(n, k) = (1 + (-1)^n).
Sum_{k=1..n-1} (-1)^(k-1)*T(n, k) = A000225(n-1).
T(2*n, n) = (-1)^(n-1)*A000984(n), n >= 1. (End)

Edited by R. J. Mathar, Jun 10 2008

A129818 Riordan array (1/(1+x), x/(1+x)^2), inverse array is A039599.

Original entry on oeis.org

1, -1, 1, 1, -3, 1, -1, 6, -5, 1, 1, -10, 15, -7, 1, -1, 15, -35, 28, -9, 1, 1, -21, 70, -84, 45, -11, 1, -1, 28, -126, 210, -165, 66, -13, 1, 1, -36, 210, -462, 495, -286, 91, -15, 1, -1, 45, -330, 924, -1287, 1001, -455, 120, -17, 1, 1, -55, 495, -1716, 3003, -3003, 1820, -680, 153, -19, 1

Offset: 0

Philippe Deléham, Jun 09 2007

This sequence is up to sign the same as A129818. - T. D. Noe, Sep 30 2011
Row sums: A057078. - Philippe Deléham, Jun 11 2007
Subtriangle of the triangle given by (0, -1, 0, -1, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (1, 0, 1, -1, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Mar 19 2012
This triangle provides the coefficients of powers of x^2 for the even-indexed Chebyshev S polynomials (see A049310): S(2*n,x) = Sum_{k=0..n} T(n,k)*x^(2*k), n >= 0. - Wolfdieter Lang, Dec 17 2012
If L(x^n) := C(n) = A000108(n) (Catalan numbers), then the polynomials P_n(x) := Sum_{k=0..n} T(n,k)*x^k are orthogonal with respect to the inner product given by (f(x),g(x)) := L(f(x)*g(x)). - Michael Somos, Jan 03 2019

			Triangle T(n,k) begins:
  n\k  0   1    2     3     4     5    6    7    8   9 10 ...
   0:  1
   1: -1   1
   2:  1  -3    1
   3: -1   6   -5     1
   4:  1 -10   15    -7     1
   5: -1  15  -35    28    -9     1
   6:  1 -21   70   -84    45   -11    1
   7: -1  28 -126   210  -165    66  -13    1
   8:  1 -36  210  -462   495  -286   91  -15    1
   9: -1  45 -330   924 -1287  1001 -455  120  -17   1
  10:  1 -55  495 -1716  3003 -3003 1820 -680  153 -19  1
  ... Reformatted by _Wolfdieter Lang_, Dec 17 2012
Recurrence from the A-sequence A115141:
15 = T(4,2) = 1*6 + (-2)*(-5) + (-1)*1.
(0, -1, 0, -1, 0, 0, ...) DELTA (1, 0, 1, -1, 0, 0, ...) begins:
  1
  0,  1
  0, -1,   1
  0,  1,  -3,   1
  0, -1,   6,  -5,  1
  0,  1, -10,  15, -7,  1
  0, -1,  15, -35, 28, -9, 1. - _Philippe Deléham_, Mar 19 2012
Row polynomial for n=3 in terms of x^2: S(6,x) = -1 + 6*x^2 -5*x^4 + 1*x^6, with Chebyshev's S polynomial. See a comment above. - _Wolfdieter Lang_, Dec 17 2012
Boas-Buck type recurrence: -35 = T(5,2) = (5/3)*(-1*1 +1*(-5) - 1*15) = -3*7 = -35. - _Wolfdieter Lang_, Jun 03 2020

Vincenzo Librandi, Rows n = 1..101, flattened
Paul Barry, A Catalan Transform and Related Transformations on Integer Sequences, Journal of Integer Sequences, Vol. 8 (2005), Article 05.4.5.
Paul Barry and Aoife Hennessy, The Euler-Seidel Matrix, Hankel Matrices and Moment Sequences, J. Int. Seq. 13 (2010) # 10.8.2, example 15.
Aoife Hennessy, A Study of Riordan Arrays with Applications to Continued Fractions, Orthogonal Polynomials and Lattice Paths, Ph. D. Thesis, Waterford Institute of Technology, Oct. 2011.

Cf. A039599, A085478, A123970, A321620.

# The function RiordanSquare is defined in A321620.
RiordanSquare((1 - sqrt(1 - 4*x))/(2*x), 10):
LinearAlgebra[MatrixInverse](%); # Peter Luschny, Jan 04 2019

max = 10; Flatten[ CoefficientList[#, y] & /@ CoefficientList[ Series[ (1 + x)/(1 + (2 - y)*x + x^2), {x, 0, max}], x]] (* Jean-François Alcover, Sep 29 2011, after Wolfdieter Lang *)

@CachedFunction
def A129818(n,k):
    if n< 0: return 0
    if n==0: return 1 if k == 0 else 0
    h = A129818(n-1,k) if n==1 else 2*A129818(n-1,k)
    return A129818(n-1,k-1) - A129818(n-2,k) - h
for n in (0..9): [A129818(n,k) for k in (0..n)] # Peter Luschny, Nov 20 2012

T(n,k) = (-1)^(n-k)*A085478(n,k) = (-1)^(n-k)*binomial(n+k,2*k).
Sum_{k=0..n} T(n,k)*A000531(k) = n^2, with A000531(0)=0. - Philippe Deléham, Jun 11 2007
Sum_{k=0..n} T(n,k)*x^k = A033999(n), A057078(n), A057077(n), A057079(n), A005408(n), A002878(n), A001834(n), A030221(n), A002315(n), A033890(n), A057080(n), A057081(n), A054320(n), A097783(n), A077416(n), A126866(n), A028230(n+1) for x = 0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, respectively. - Philippe Deléham, Nov 19 2009
O.g.f.: (1+x)/(1+(2-y)*x+x^2). - Wolfdieter Lang, Dec 15 2010
O.g.f. column k with leading zeros (Riordan array, see NAME): (1/(1+x))*(x/(1+x)^2)^k, k >= 0. - Wolfdieter Lang, Dec 15 2010
From Wolfdieter Lang, Dec 20 2010: (Start)
Recurrences from the Z- and A-sequences for Riordan arrays. See the W. Lang link under A006232 for details and references.
T(n,0) = -1*T(n-1,0), n >= 1, from the o.g.f. -1 for the Z-sequence (trivial result).
T(n,k) = Sum_{j=0..n-k} A(j)*T(n-1,k-1+j), n >= k >= 1, with A(j):= A115141(j) = [1,-2,-1,-2,-5,-14,...], j >= 0 (o.g.f. 1/c(x)^2 with the A000108 (Catalan) o.g.f. c(x)). (End)
T(n,k) = (-1)^n*A123970(n,k). - Philippe Deléham, Feb 18 2012
T(n,k) = -2*T(n-1,k) + T(n-1,k-1) - T(n-2,k), T(0,0) = T(1,1) = 1, T(1,0) = -1, T(n,k) = 0 if k < 0 or if k > n. - Philippe Deléham, Mar 19 2012
A039599(m,n) = Sum_{k=0..n} T(n,k) * C(k+m) where C(n) are the Catalan numbers. - Michael Somos, Jan 03 2019
Equals the matrix inverse of the Riordan square (cf. A321620) of the Catalan numbers. - Peter Luschny, Jan 04 2019
Boas-Buck type recurrence for column k >= 0 (see Aug 10 2017 comment in A046521 with references): T(n,k) = ((1 + 2*k)/(n - k))*Sum_{j = k..n-1} (-1)^(n-j)*T(j,k), with input T(n,n) = 1, and T(n,k) = 0 for n < k. - Wolfdieter Lang, Jun 03 2020

A128018 Expansion of (1-4x)/(1-2x+4*x^2).

Original entry on oeis.org

1, -2, -8, -8, 16, 64, 64, -128, -512, -512, 1024, 4096, 4096, -8192, -32768, -32768, 65536, 262144, 262144, -524288, -2097152, -2097152, 4194304, 16777216, 16777216, -33554432, -134217728, -134217728, 268435456, 1073741824, 1073741824, -2147483648, -8589934592

Offset: 0

Paul Barry, Feb 11 2007

Hankel transform of A128014(n+1). Binomial transform of A128019.
Hankel transform of A002426(n+1). - Paul Barry, Mar 15 2008
Hankel transform of A007971(n+1). - Paul Barry, Sep 30 2009
Hankel transform of A103970 is a(n)/4^C(n+1,2). - Paul Barry, Nov 20 2009
The real part of Q^(n+1), where Q is the quaternion 1+i+j+k. - Stanislav Sykora, Jun 11 2012.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (2,-4).

Cf. A002426, A007971, A128014, A103970, A128019, A138340.

CoefficientList[Series[(1 - 4*x)/(1 - 2*x + 4*x^2), {x,0,50}], x] (* or *) LinearRecurrence[{2,-4},{1,-2},50] (* G. C. Greubel, Feb 28 2017 *)

x='x+O('x^50); Vec((1-4*x)/(1-2*x+4*x^2)) \\ G. C. Greubel, Feb 28 2017

a(n) = A138340(n)/2^n. - Philippe Deléham, Nov 14 2008
a(n) = 2^(n+1)*cos(Pi*(n+1)/3). - Richard Choulet, Nov 19 2008
From Paul Barry, Oct 21 2009: (Start)
a(n) = Sum_{k=0..floor((n+1)/2)} C(n+1,2*k)*(-3)^k.
a(n) = ((1+i*sqrt(3))^(n+1) + (1-i*sqrt(3))^(n+1))/2, i=sqrt(-1). (End)
G.f.: G(0)/(2*x)-1/x, where G(k)= 1 + 1/(1 - x*(3*k+1)/(x*(3*k+4) + 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 27 2013
a(n) = 2^n*A057079(n+2). - R. J. Mathar, Mar 04 2018
Sum_{n>=0} 1/a(n) = 1/3. - Amiram Eldar, Feb 14 2023

A130755 Binomial transform of periodic sequence (3, 1, 2).

Original entry on oeis.org

3, 4, 7, 15, 32, 65, 129, 256, 511, 1023, 2048, 4097, 8193, 16384, 32767, 65535, 131072, 262145, 524289, 1048576, 2097151, 4194303, 8388608, 16777217, 33554433, 67108864, 134217727, 268435455, 536870912, 1073741825, 2147483649

Offset: 0

Paul Curtz, Jul 13 2007

The third sequence of "less twisted numbers"; this sequence, A130750 and A130752 form a "suite en trio" (cf. reference, p. 130).
First differences of A130752, second differences of A130750.
Sequence equals its third differences:
  3     4     7    15    32    65   129   256   511  1023
     1     3     8    17    33    64   127   255   512
        2     5     9    16    31    63   128   257
           3     4     7    15    32    65   129

P. Curtz, Exercise Book, manuscript, 1995.

G. C. Greubel, Table of n, a(n) for n = 0..1000 (terms 0..300 from Vincenzo Librandi)
Index entries for linear recurrences with constant coefficients, signature (3,-3,2).

Cf. A010882 (periodic (1, 2, 3)), A128834 (periodic (0, 1, 1, 0, -1, -1)), A057079 (periodic (1, 2, 1, -1, -2, -1)), A130750 (first differences), A130752 (second differences).

m:=31; S:=[ [3, 1, 2][(n-1) mod 3 +1]: n in [1..m] ]; [ &+[ Binomial(i-1, k-1)*S[k]: k in [1..i] ]: i in [1..m] ]; // Klaus Brockhaus, Aug 03 2007

I:=[3,4,7]; [n le 3 select I[n] else 3*Self(n-1) - 3*Self(n-2) + 2*Self(n-3): n in [1..30]]; // G. C. Greubel, Jan 15 2018

CoefficientList[Series[(3-5*x+4*x^2)/((1-2*x)*(1-x+x^2)), {x, 0, 50}], x] (* or *) LinearRecurrence[{3,-3,2}, {3,4,7}, 30] (* G. C. Greubel, Jan 15 2018 *)

{m=31; v=vector(m); v[1]=3; v[2]=4; v[3]=7; for(n=4, m, v[n]=3*v[n-1]-3*v[n-2]+2*v[n-3]); v} \\ Klaus Brockhaus, Aug 03 2007

{for(n=0, 30, print1(2^(n+1)+[1, 0, -1, -1, 0, 1][n%6+1], ","))} \\ Klaus Brockhaus, Aug 03 2007

G.f.: (3-5*x+4*x^2)/((1-2*x)*(1-x+x^2)).
a(0) = 3; a(1) = 4; a(2) = 7; for n > 2, a(n) = 3*a(n-1) - 3*a(n-2) + 2*a(n-3).
a(n) = 2^(n+1) + A128834(n+2).
a(0) = 3; for n > 0, a(n) = 2*a(n-1) + A057079(n+3).

Edited and extended by Klaus Brockhaus, Aug 03 2007

A007039 Number of cyclic binary n-bit strings with no alternating substring of length > 2.

Original entry on oeis.org

2, 2, 2, 6, 12, 20, 30, 46, 74, 122, 200, 324, 522, 842, 1362, 2206, 3572, 5780, 9350, 15126, 24474, 39602, 64080, 103684, 167762, 271442, 439202, 710646, 1149852, 1860500, 3010350, 4870846, 7881194, 12752042, 20633240, 33385284, 54018522, 87403802, 141422322

Offset: 1

N. J. A. Sloane

John W. Layman observes that the second differences give the sequence shifted to the right.

			G.f. = 2*x + 2*x^2 + 2*x^3 + 6*x^4 + 12*x^5 + 20*x^6 + 30*x^7 + 46*x^8 + ...

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Vincenzo Librandi, Table of n, a(n) for n = 1000
Z. Agur, A. S. Fraenkel, and S. T. Klein, The number of fixed points of the majority rule, Discr. Math., 70 (1988), 295-302.
Matthew Macauley, Jon McCammond, and Henning S. Mortveit, Dynamics groups of asynchronous cellular automata, Journal of Algebraic Combinatorics, Vol 33, No 1 (2011), pp. 11-35.
A. McLeod and W. O. J. Moser, Counting cyclic binary strings, Math. Mag., 80 (No. 1, 2007), 29-37.
W. O. J. Moser, On cyclic binary n-bit strings, Report from the Department of Mathematics and Statistics, McGill University, 1991. (Annotated scanned copy)
W. O. J. Moser, Cyclic binary strings without long runs of like (alternating) bits, Fibonacci Quart. 31 (1993), no. 1, 2-6.
E. Munarini and N. Z. Salvi, Circular Binary Strings without Zigzags, Integers: Electronic Journal of Combinatorial Number Theory 3 (2003), #A19. See relation (8) on page 5.
Index entries for linear recurrences with constant coefficients, signature (2,-1,0,1).

Cf. A000032, A007040, A057079, A111734.

CoefficientList[Series[2*(1+x)*(1-2*x+2*x^2)/((1-x+x^2)*(1-x-x^2)),{x,0,40}],x] (* Vincenzo Librandi, Apr 16 2012 *)
a[n_ /; n<4] = 2; a[4] = 6; a[n_] := a[n] = 2*a[n-1] - a[n-2] + a[n-4]; Array[a, 39] (* Jean-François Alcover, Oct 08 2017 *)

Vec(2*x*(1-x+2*x^3)/((1-x-x^2)*(1-x+x^2))+O(x^66)) \\ Joerg Arndt, Oct 27 2015

For n >= 5, a(n) = 2a(n-1) - a(n-2) + a(n-4). - David W. Wilson
G.f.: 2*x*(1+x)*(1-2*x+2*x^2)/((1-x+x^2)*(1-x-x^2)). - Colin Barker, Mar 28 2012
a(n) = A000032(n) + A057079(n + 1). - John M. Campbell, Dec 29 2016
a(n) = abs(A111734(n)). - Alois P. Heinz, Oct 08 2017
E.g.f.: 2*exp(x/2)*(cos(sqrt(3)*x/2) + cosh(sqrt(5)*x/2)) - 4. - Stefano Spezia, Mar 09 2025

cp's OEIS Frontend

A024494 a(n) = C(n,1) + C(n,4) + ... + C(n, 3*floor(n/3) + 1).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

SageMath

Formula

A294099 Rectangular array read by (upward) antidiagonals: A(n,k) = Sum_{j=0..k} (-1)^floor(j/2)*binomial(k-floor((j+1)/2), floor(j/2))*n^(k-j), n >= 1, k >= 0.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A174709 Partial sums of floor(n/6).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

A087204 Period 6: repeat [2, 1, -1, -2, -1, 1].

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

Sage

Formula

Extensions

A100051 A Chebyshev transform of 1,1,1,...

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A135356 Triangle T(n,k) read by rows: coefficients in the recurrence of sequences which equal their n-th differences.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

A294099 Rectangular array read by (upward) antidiagonals: A(n,k) = Sum_{j=0..k} (-1)^floor(j/2)binomial(k-floor((j+1)/2), floor(j/2))n^(k-j), n >= 1, k >= 0.

A128018 Expansion of (1-4x)/(1-2x+4*x^2).