A099586 -id:A099586 - cp's OEIS frontend

A307039 Square array A(n,k), n >= 0, k >= 1, read by antidiagonals, where column k is the expansion of g.f. ((1-x)^(k-1))/((1-x)^k+x^k).

1, 1, 0, 1, 1, 0, 1, 1, 0, 0, 1, 1, 1, -2, 0, 1, 1, 1, 0, -4, 0, 1, 1, 1, 1, -3, -4, 0, 1, 1, 1, 1, 0, -9, 0, 0, 1, 1, 1, 1, 1, -4, -18, 8, 0, 1, 1, 1, 1, 1, 0, -14, -27, 16, 0, 1, 1, 1, 1, 1, 1, -5, -34, -27, 16, 0, 1, 1, 1, 1, 1, 1, 0, -20, -68, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, -6, -55, -116, 81, -32, 0

Offset: 0

Seiichi Manyama, Mar 21 2019

			Square array begins:
   1,  1,   1,    1,    1,   1,   1,  1, ...
   0,  1,   1,    1,    1,   1,   1,  1, ...
   0,  0,   1,    1,    1,   1,   1,  1, ...
   0, -2,   0,    1,    1,   1,   1,  1, ...
   0, -4,  -3,    0,    1,   1,   1,  1, ...
   0,  0, -18,  -14,   -5,   0,   1,  1, ...
   0,  8, -27,  -34,  -20,  -6,   0,  1, ...
   0, 16, -27,  -68,  -55, -27,  -7,  0, ...
   0, 16,   0, -116, -125, -83, -35, -8, ...

Seiichi Manyama, Antidiagonals n = 0..139, flattened

Columns 1-9 give A000007, A146559, A057681, A099586, A289306, A307040, A307041, A307044, A307045.

Cf. A306846, A306914.

T[n_, k_] := Sum[(-1)^j * Binomial[n, k*j], {j, 0, Floor[n/k]}]; Table[T[n-k, k], {n, 0, 13}, {k, n, 1, -1}] // Flatten (* Amiram Eldar, May 20 2021 *)

A(n,k) = Sum_{j=0..floor(n/k)} (-1)^j * binomial(n,k*j).

A099587 a(n) = coefficient of x in (1+x)^n mod (1+x^4).

Original entry on oeis.org

0, 1, 2, 3, 4, 4, 0, -14, -48, -116, -232, -396, -560, -560, 0, 1912, 6528, 15760, 31520, 53808, 76096, 76096, 0, -259808, -887040, -2141504, -4283008, -7311552, -10340096, -10340096, 0, 35303296, 120532992, 290992384, 581984768

Offset: 0

Ralf Stephan, Oct 24 2004

{A099586, A099587, A099588, A099589} is the difference analog of the trigonometric functions {k_1(x), k_2(x), k_3(x), k_4(x)} of order 4.

For the definition, see [Erdelyi] and the Shevelev link. - Vladimir Shevelev, Jul 03 2017

A. Erdelyi, Higher Transcendental Functions, McGraw-Hill, 1955, Vol. 3, Chapter XVIII.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Vladimir Shevelev, Combinatorial identities generated by difference analogs of hyperbolic and trigonometric functions of order n, arXiv:1706.01454 [math.CO], 2017.
Vladimir Shevelev, Coefficient of x^k in ((x+1)^n modulo x^N+1), seqfan, Thu Jul 20 2017.
G. Tollisen and T. Lengyel, A Congruential Identity and the 2-adic Order of Lacunary Sums of Binomial Coefficients, Integers 4 (2004), #A4.
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-2).

Cf. A099586, A099588, A099589.

Cf. A057681, A099586, A139398, A133476, A139714, A139748, A139761.

RecurrenceTable[{a[1]=1, a[2]=2, a[3]=3, a[4]=4, a[n] = 4*a[n-1] - 6*a[n-2] + 4*a[n-3] - 2*a[n-4]}, a, {n, 1, 100}] (* G. C. Greubel, Nov 09 2015 *)
a[n_] := n*HypergeometricPFQ[{(1-n)/4, (2-n)/4, (3-n)/4, (4-n)/4}, {1/2, 3/4, 5/4}, -1]; Array[a, 40, 0] (* Jean-François Alcover, Jul 20 2017, from Vladimir Shevelev's first formula *)
LinearRecurrence[{4,-6,4,-2},{0,1,2,3},50] (* Harvey P. Dale, Mar 27 2022 *)

PARI
```
a(n) = polcoeff(((1+x)^n)%(x^4+1),1)
```

G.f.: x*(x-1)^2 / (2*x^4-4*x^3+6*x^2-4*x+1). - Colin Barker, Jul 15 2013
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - 2*a(n-4). - G. C. Greubel, Nov 09 2015
From Vladimir Shevelev, Jun 29 2017: (Start)
a(n) = Sum_{k >= 0}(-1)^k*binomial(n,4*k+1).
a(n) = round((2+sqrt(2))^(n/2)*cos(Pi*(n-2)/8)/2), where round(x) is the integer nearest to x.
a(n+m) = a(n)*K_1(m) + K_1(n)*a(m) - K_4(n)*K_3(m) - K_3(n)*K_4(m), where K_1 is A099586, K_3=A099588, and K_4=A099589.
(End)
a(n) = A099589(n+2)-2*A099589(n+1)+A099589(n). - R. J. Mathar, Jun 28 2025

a(0)=0 added by N. J. A. Sloane, Jun 30 2017

A099589 Expansion of x^3 / (1 - 4x + 6x^2 - 4x^3 + 2x^4).

Original entry on oeis.org

0, 0, 0, 1, 4, 10, 20, 34, 48, 48, 0, -164, -560, -1352, -2704, -4616, -6528, -6528, 0, 22288, 76096, 183712, 367424, 627232, 887040, 887040, 0, -3028544, -10340096, -24963200, -49926400, -85229696, -120532992, -120532992, 0, 411525376, 1405035520, 3392055808

Offset: 0

Ralf Stephan, Oct 24 2004

{A099586, A099587, A099588, A099589} is the difference analog of the trigonometric functions {k_1(x), k_2(x), k_3(x), k_4(x)} of order 4. For the definition, see [Erdelyi] and the Shevelev link. - Vladimir Shevelev, Jul 04 2017

A. Erdelyi, Higher Transcendental Functions, McGraw-Hill, 1955, Vol. 3, Chapter XVIII.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Vladimir Shevelev, Combinatorial identities generated by difference analogs of hyperbolic and trigonometric functions of order n, arXiv:1706.01454 [math.CO], 2017.
Vladimir Shevelev, Coefficient of x^k in ((x+1)^n modulo x^N+1), seqfan, Thu Jul 20 2017.
G. Tollisen and T. Lengyel, A congruential identity and the 2-adic order of lacunary sums of binomial coefficients, Integers 4 (2004), #A4.
Maran van Heesch, The multiplicative complexity of symmetric functions over a field with characteristic p, Thesis, 2014.
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-2).

Cf. A099586, A099587, A099588.

Round@Table[(1/(2*Sqrt[2]))*((2-Sqrt[2])^(n/2)*(Cos[3*Pi*n/8] + Sin[3*Pi*n/8]) + (2+Sqrt[2])^(n/2)*(Sin[Pi*n/8] - Cos[Pi*n/8])), {n, 0, 40}] (* G. C. Greubel, Nov 07 2015 *)
RecurrenceTable[{a[n] == 4*a[n-1] - 6*a[n-2] + 4*a[n-3] - 2*a[n-4], a[0]==0, a[1]==0, a[2]==0, a[3]==1}, a, {n, 0, 40}] (* G. C. Greubel, Nov 10 2015 *)
Table[Sum[(-1)^k*Binomial[n, 4 k + 3], {k, 0, n}], {n, 0, 37}] (* Michael De Vlieger, Jun 30 2017 *)
a[n_] := n*(n-1)*(n-2)/6 HypergeometricPFQ[{(3-n)/4, (4-n)/4, (5-n)/4, (6-n)/4}, {5/4, 3/2, 7/4}, -1]; Array[a, 40, 0] (* Jean-François Alcover, Jul 20 2017, from Vladimir Shevelev's first formula *)

PARI
```
a(n) = polcoeff(((1+x)^n)%(x^4+1),3)
```

concat([0, 0], Vec(x^3/((1-x)^4+x^4) + O(x^50))) \\ Altug Alkan, Nov 08 2015

a(n) = sum(t=0, (n-3)\4, (-1)^t*binomial(n,4*t+3)); \\ Michel Marcus, Jun 30 2017

G.f.: x^3/((1-x)^4 + x^4), the binomial transform of x^3/(1+x^4). - Paul Barry, Apr 01 2005
Coefficient of x^3 in (1+x)^n mod (1 + x^4).
a(n) = (1/(2*sqrt(2)))*((2-sqrt(2))^(n/2)*(cos(3*Pi*n/8) + sin(3*Pi*n/8)) + (2+sqrt(2))^(n/2)*(sin(Pi*n/8) - cos(Pi*n/8))). - Paul Barry, Apr 01 2005
From Colin Barker, Nov 08 2015: (Start)
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - 2*a(n-4) for n > 4.
G.f.: x^3 / (2*x^4 - 4*x^3 + 6*x^2 - 4*x + 1). (End)
From Vladimir Shevelev, Jul 04 2017: (Start)
a(n) = Sum_{t >= 0} (-1)^t*binomial(n,4*t+3).
a(n) = round((2+sqrt(2))^(n/2)*cos(Pi*(n-6)/8)/2), where round(x) is the integer nearest to x.
a(n+m) = a(n)*K_1(m) + K_3(n)*K_2(m) + K_2(n)*K_3(m) + K_1(n)*a(m), where
K_1 is A099586, K_2 is A099587, K_3 is A099588. (End)

a(0)=0 added by N. J. A. Sloane, Jul 04 2017

A099588 Coefficient of x^2 in (1+x)^n mod 1+x^4.

Original entry on oeis.org

0, 0, 1, 3, 6, 10, 14, 14, 0, -48, -164, -396, -792, -1352, -1912, -1912, 0, 6528, 22288, 53808, 107616, 183712, 259808, 259808, 0, -887040, -3028544, -7311552, -14623104, -24963200, -35303296, -35303296, 0, 120532992, 411525376, 993510144, 1987020288, 3392055808, 4797091328, 4797091328

Offset: 0

Ralf Stephan, Oct 24 2004

{A099586, A099587, A099588, A099589} is the difference analog of the trigonometric functions {k_1(x), k_2(x), k_3(x), k_4(x)} of order 4. For the definition, see [Erdelyi] and the Shevelev link. - Vladimir Shevelev, Jul 04 2017

A. Erdelyi, Higher Transcendental Functions, McGraw-Hill, 1955, Vol. 3, Chapter XVIII.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Vladimir Shevelev, Combinatorial identities generated by difference analogs of hyperbolic and trigonometric functions of order n, arXiv:1706.01454 [math.CO], 2017.
Vladimir Shevelev, Coefficient of x^k in ((x+1)^n modulo x^N+1), seqfan, Thu Jul 20 2017.
G. Tollisen and T. Lengyel, A congruential identity and the 2-adic order of lacunary sums of binomial coefficients, Integers 4 (2004), #A4.
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-2).

Cf. A099586, A099587, A099589.

f:= rectoproc({rec, a(0)=0,a(1)=0,a(2)=1,a(3)=3},a(n),remember):
map(f, [$0..100]); # Robert Israel, Jun 30 2017

RecurrenceTable[{a[n]==4*a[n-1] - 6*a[n-2] + 4*a[n-3] - 2*a[n-4], a[1]=0, a[2]=1, a[3]=3, a[4]=6}, a, {n, 1, 200}] (* G. C. Greubel, Nov 10 2015 *)
Table[Sum[(-1)^k*Binomial[n, 4 k + 2], {k, 0, n}], {n, 0, 36}] (* Michael De Vlieger, Jun 30 2017 *)
a[n_] := n*(n-1)/2 HypergeometricPFQ[{(2-n)/4, (3-n)/4, (4-n)/4, (5-n)/4}, {3/4, 5/4, 3/2}, -1]; Array[a, 40, 0] (* Jean-François Alcover, Jul 20 2017, from Vladimir Shevelev's first formula *)

x='x+O('x^55); concat([0, 0], Vec(-x^2*(x-1)/(2*x^4-4*x^3+6*x^2-4*x+1))) \\ Altug Alkan, Nov 11 2015

a(n) = sum(t=0, (n-2)\4, (-1)^t*binomial(n,4*t+2)); \\ Michel Marcus, Jun 30 2017

a(n)=polcoeff(lift(Mod(1+x,1+x^4)^n),2); \\ Joerg Arndt, Feb 22 2018

G.f.: -x^2*(x-1) / (2*x^4-4*x^3+6*x^2-4*x+1). - Colin Barker, Jul 15 2013
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - 2*a(n-4). - G. C. Greubel, Nov 10 2015
a(n) = (1/2)*((2+sqrt(2))^(n/2) * sin(n*Pi/8) - (2-sqrt(2))^(n/2)*sin(3*n*Pi/8)). - G. C. Greubel, Nov 10 2015
From Vladimir Shevelev, Jul 04 2017 (Start)
a(n) = Sum_{k>=0}(-1)^k*binomial(n,4*k+2);
a(n) = round((2+sqrt(2))^(n/2)*cos(Pi*(n-4)/8)/2), where round(x) is the integer nearest to x;
a(n+m) = a(n)*K_1(m) + K_2(n)*K_2(m) + K_1(n)*a(m) - K_4(n)*K_4(m), where K_1 is A099586, K_2 is A099587, K_4 is A099589. (End)

a(0)=0 added by N. J. A. Sloane, Jul 04 2017

a(673) in b-file corrected by Andrew Howroyd, Feb 21 2018

A290286 Determinant of circulant matrix of order 4 with entries in the first row (-1)^jSum_{k>=0}(-1)^kbinomial(n, 4*k+j), j=0,1,2,3.

Original entry on oeis.org

1, 0, 0, 0, -1008, -37120, -473600, 0, 63996160, 702013440, 2893578240, 0, -393379835904, -12971004067840, -160377313820672, 0, 21792325059543040, 239501351489372160, 987061897553510400, 0, -134124249770961666048, -4422152303189489090560

Offset: 0

Vladimir Shevelev and Peter J. C. Moses, Jul 26 2017

sign

In the Shevelev link the author proved that, for odd N>=3 and every n>=1, the determinant of circulant matrix of order N with entries in the first row (-1)^j*Sum{k>=0}(-1)^k*binomial(n, N*k+j), j=0..N-1, is 0.

This sequence shows what happens for the first even N>3.

Vladimir Shevelev, Combinatorial identities generated by difference analogs of hyperbolic and trigonometric functions of order n, arXiv:1706.01454 [math.CO], 2017.
Wikipedia, Circulant matrix

Cf. A099586 (prefixed by a(0)=1), A099587, A099588, A099589, A290285.

seq(LinearAlgebra:-Determinant(Matrix(4,shape=Circulant[seq((-1)^j*
add((-1)^k*binomial(n, 4*k+j),k=0..n/4),j=0..3)])),n=0..50); # Robert Israel, Jul 26 2017

ro[n_] := Table[Sum[(-1)^(j+k) Binomial[n, 4k+j], {k, 0, n/4}], {j, 0, 3}];
M[n_] := Table[RotateRight[ro[n], m], {m, 0, 3}];
a[n_] := Det[M[n]];
Table[a[n], {n, 0, 21}] (* Jean-François Alcover, Aug 09 2018 *)

from sympy.matrices import Matrix
from sympy import binomial
def mj(j, n): return (-1)**j*sum((-1)**k*binomial(n, 4*k + j) for k in range(n//4 + 1))
def a(n):
    m=Matrix(4, 4, lambda i,j: mj((i-j)%4,n))
    return m.det()
print([a(n) for n in range(22)]) # Indranil Ghosh, Jul 31 2017

a(n) = 0 for n == 3 (mod 4).

G.f. (empirical): (1/8)*(68*x^2+1)/(16*x^4+136*x^2+1)+(1/4)*(68*x^2-8*x+1)/(16*x^4+64*x^3+128*x^2-16*x+1)+(1/2)*(12*x^2+1)/(16*x^4+24*x^2+1)+3/(8*(4*x^2+1))-(1/4)*(12*x^2-4*x+1)/(16*x^4-32*x^3+32*x^2-8*x+1)-(1/4)*(4*x^2+1)/(16*x^4+1)+(1/4)*(12*x^2+4*x+1)/(16*x^4+32*x^3+32*x^2+8*x+1). - Robert Israel, Jul 26 2017

A348309 a(n) = Sum_{k=0..floor(n/8)} (-1)^k * binomial(n-4k,4k).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 0, -4, -14, -34, -69, -125, -209, -329, -493, -705, -955, -1199, -1324, -1092, -56, 2560, 8025, 18313, 36353, 66273, 113525, 184653, 286257, 422377, 589028, 763912, 888378, 837502, 372835, -928725, -3776537, -9302337, -19226889, -36034869, -63099331, -104630831, -165212760

Offset: 0

Seiichi Manyama, Oct 11 2021

Seiichi Manyama, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1,0,0,0,-1).

Cf. A348308, A348310.

Cf. A099586, A348289.

LinearRecurrence[{4, -6, 4, -1, 0, 0, 0, -1}, {1, 1, 1, 1, 1, 1, 1, 1}, 45] (* Amiram Eldar, Oct 11 2021 *)

a(n) = sum(k=0, n\8, (-1)^k*binomial(n-4*k, 4*k));

my(N=66, x='x+O('x^N)); Vec((1-x)^3/((1-x)^4+x^8))

G.f.: (1-x)^3/((1-x)^4 + x^8).

a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4) - a(n-8).

A099590 2^(n-1) times coefficient of x in (1+x)^n mod U(n,x), U the Chebyshev polynomials.

Original entry on oeis.org

0, 4, 14, 32, 77, 192, 452, 1024, 2299, 5120, 11270, 24576, 53241, 114688, 245768, 524288, 1114103, 2359296, 4980746, 10485760, 22020085, 46137344, 96469004, 201326592, 419430387, 872415232, 1811939342, 3758096384, 7784628209

Offset: 1

Ralf Stephan, Oct 24 2004

nonn

Cf. A099586.

Apparently, a(n+1) = (n+1)2^n - (n+2)/4 * (I^n + (-I)^n).

Empirical g.f.: x^2*(x+2)*(x^2-2*x+2) / ((2*x-1)^2*(x^2+1)^2). - Colin Barker, Jul 16 2013

cp's OEIS Frontend

A307039 Square array A(n,k), n >= 0, k >= 1, read by antidiagonals, where column k is the expansion of g.f. ((1-x)^(k-1))/((1-x)^k+x^k).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A099587 a(n) = coefficient of x in (1+x)^n mod (1+x^4).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A099589 Expansion of x^3 / (1 - 4*x + 6*x^2 - 4*x^3 + 2*x^4).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Mathematica

PARI

PARI

PARI

Formula

Extensions

A099588 Coefficient of x^2 in (1+x)^n mod 1+x^4.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

PARI

Formula

Extensions

A290286 Determinant of circulant matrix of order 4 with entries in the first row (-1)^j*Sum_{k>=0}(-1)^k*binomial(n, 4*k+j), j=0,1,2,3.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

Python

Formula

A348309 a(n) = Sum_{k=0..floor(n/8)} (-1)^k * binomial(n-4*k,4*k).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

A099590 2^(n-1) times coefficient of x in (1+x)^n mod U(n,x), U the Chebyshev polynomials.

Views

A099589 Expansion of x^3 / (1 - 4x + 6x^2 - 4x^3 + 2x^4).

A290286 Determinant of circulant matrix of order 4 with entries in the first row (-1)^jSum_{k>=0}(-1)^kbinomial(n, 4*k+j), j=0,1,2,3.

A348309 a(n) = Sum_{k=0..floor(n/8)} (-1)^k * binomial(n-4k,4k).