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

A049310 Triangle of coefficients of Chebyshev's S(n,x) := U(n,x/2) polynomials (exponents in increasing order).

Original entry on oeis.org

1, 0, 1, -1, 0, 1, 0, -2, 0, 1, 1, 0, -3, 0, 1, 0, 3, 0, -4, 0, 1, -1, 0, 6, 0, -5, 0, 1, 0, -4, 0, 10, 0, -6, 0, 1, 1, 0, -10, 0, 15, 0, -7, 0, 1, 0, 5, 0, -20, 0, 21, 0, -8, 0, 1, -1, 0, 15, 0, -35, 0, 28, 0, -9, 0, 1, 0, -6, 0, 35, 0, -56, 0, 36, 0, -10, 0, 1, 1, 0, -21, 0, 70, 0, -84, 0
Offset: 0

Views

Author

Keywords

Comments

G.f. for row polynomials S(n,x) (signed triangle): 1/(1-x*z+z^2). Unsigned triangle |a(n,m)| has Fibonacci polynomials F(n+1,x) as row polynomials with g.f. 1/(1-x*z-z^2). |a(n,m)| triangle has rows of Pascal's triangle A007318 in the even-numbered diagonals (odd-numbered ones have only 0's).
Row sums (unsigned triangle) A000045(n+1) (Fibonacci). Row sums (signed triangle) S(n,1) sequence = periodic(1,1,0,-1,-1,0) = A010892.
Alternating row sums A049347(n) = S(n,-1) = periodic(1,-1,0). - Wolfdieter Lang, Nov 04 2011
S(n,x) is the characteristic polynomial of the adjacency matrix of the n-path. - Michael Somos, Jun 24 2002
S(n,x) is also the matching polynomial of the n-path. - Eric W. Weisstein, Apr 10 2017
|T(n,k)| = number of compositions of n+1 into k+1 odd parts. Example: |T(7,3)| = 10 because we have (1,1,3,3), (1,3,1,3), (1,3,3,1), (3,1,1,3), (3,1,3,1), (3,3,1,1), (1,1,1,5), (1,1,5,1), (1,5,1,1) and (5,1,1,1). - Emeric Deutsch, Apr 09 2005
S(n,x)= R(n,x) + S(n-2,x), n >= 2, S(-1,x)=0, S(0,x)=1, R(n,x):=2*T(n,x/2) = Sum_{m=0..n} A127672(n,m)*x^m (monic integer Chebyshev T-Polynomials). This is the rewritten so-called trace of the transfer matrix formula for the T-polynomials. - Wolfdieter Lang, Dec 02 2010
In a regular N-gon inscribed in a unit circle, the side length is d(N,1) = 2*sin(Pi/N). The length ratio R(N,k):=d(N,k)/d(N,1) for the (k-1)-th diagonal, with k from {2,3,...,floor(N/2)}, N >= 4, equals S(k-1,x) = sin(k*Pi/N)/sin(Pi/N) with x=rho(N):=R(N,2) = 2*cos(Pi/N). Example: N=7 (heptagon), rho=R(7,2), sigma:=R(N,3) = S(2,rho) = rho^2 - 1. Motivated by the quoted paper by P. Steinbach. - Wolfdieter Lang, Dec 02 2010
From Wolfdieter Lang, Jul 12 2011: (Start)
In q- or basic analysis, q-numbers are [n]_q := S(n-1,q+1/q) = (q^n-(1/q)^n)/(q-1/q), with the row polynomials S(n,x), n >= 0.
The zeros of the row polynomials S(n-1,x) are (from those of Chebyshev U-polynomials):
x(n-1;k) = +- t(k,rho(n)), k = 1..ceiling((n-1)/2), n >= 2, with t(n,x) the row polynomials of A127672 and rho(n):= 2*cos(Pi/n). The simple vanishing zero for even n appears here as +0 and -0.
Factorization of the row polynomials S(n-1,x), x >= 1, in terms of the minimal polynomials of cos(2 Pi/2), called Psi(n,x), with coefficients given by A181875/A181876:
S(n-1,x) = (2^(n-1))*Product_{n>=1}(Psi(d,x/2), 2 < d | 2n).
(From the rewritten eq. (3) of the Watkins and Zeitlin reference, given under A181872.) [See the W. Lang ArXiv link, Proposition 9, eq. (62). - Wolfdieter Lang, Apr 14 2018]
(End)
The discriminants of the S(n,x) polynomials are found in A127670. - Wolfdieter Lang, Aug 03 2011
This is an example for a subclass of Riordan convolution arrays (lower triangular matrices) called Bell arrays. See the L. W. Shapiro et al. reference under A007318. If a Riordan array is named (G(z),F(z)) with F(z)=z*Fhat(z), the o.g.f. for the row polynomials is G(z)/(1-x*z*Fhat(z)), and it becomes a Bell array if G(z)=Fhat(z). For the present Bell type triangle G(z)=1/(1+z^2) (see the o.g.f. comment above). This leads to the o.g.f. for the column no. k, k >= 0, x^k/(1+x^2)^(k+1) (see the formula section), the one for the row sums and for the alternating row sums (see comments above). The Riordan (Bell) A- and Z-sequences (defined in a W. Lang link under A006232, with references) have o.g.f.s 1-x*c(x^2) and -x*c(x^2), with the o.g.f. of the Catalan numbers A000108. Together they lead to a recurrence given in the formula section. - Wolfdieter Lang, Nov 04 2011
The determinant of the N x N matrix S(N,[x[1], ..., x[N]]) with elements S(m-1,x[n]), for n, m = 1, 2, ..., N, and for any x[n], is identical with the determinant of V(N,[x[1], ..., x[N]]) with elements x[n]^(m-1) (a Vandermondian, which equals Product_{1 <= i < j<= N} (x[j] - x[i])). This is a special instance of a theorem valid for any N >= 1 and any monic polynomial system p(m,x), m>=0, with p(0,x) = 1. For this theorem see the Vein-Dale reference, p. 59. Thanks to L. Edson Jeffery for an email asking for a proof of the non-singularity of the matrix S(N,[x[1], ...., x[N]]) if and only if the x[j], j = 1..N, are pairwise distinct. - Wolfdieter Lang, Aug 26 2013
These S polynomials also appear in the context of modular forms. The rescaled Hecke operator T*n = n^((1-k)/2)*T_n acting on modular forms of weight k satisfies T*(p^n) = S(n, T*p), for each prime p and positive integer n. See the Koecher-Krieg reference, p. 223. - _Wolfdieter Lang, Jan 22 2016
For a shifted o.g.f. (mod signs), its compositional inverse, and connections to Motzkin and Fibonacci polynomials, non-crossing partitions and other combinatorial structures, see A097610. - Tom Copeland, Jan 23 2016
From M. Sinan Kul, Jan 30 2016; edited by Wolfdieter Lang, Jan 31 2016 and Feb 01 2016: (Start)
Solutions of the Diophantine equation u^2 + v^2 - k*u*v = 1 for integer k given by (u(k,n), v(k,n)) = (S(n,k), S(n-1,k)) because of the Cassini-Simson identity: S(n,x)^2 - S(n+1,x)*S(n-1, x) = 1, after use of the S-recurrence. Note that S(-n, x) = -S(-n-2, x), n >= 1, and the periodicity of some S(n, k) sequences.
Hence another way to obtain the row polynomials would be to take powers of the matrix [x, -1; 1,0]: S(n, x) = (([x, -1; 1, 0])^n)[1,1], n >= 0.
See also a Feb 01 2016 comment on A115139 for a well-known S(n, x) sum formula.
Then we have with the present T triangle
A039834(n) = -i^(n+1)*T(n-1, k) where i is the imaginary unit and n >= 0.
A051286(n) = Sum_{i=0..n} T(n,i)^2 (see the Philippe Deléham, Nov 21 2005 formula),
A181545(n) = Sum_{i=0..n+1} abs(T(n,i)^3),
A181546(n) = Sum_{i=0..n+1} T(n,i)^4,
A181547(n) = Sum_{i=0..n+1} abs(T(n,i)^5).
S(n, 0) = A056594(n), and for k = 1..10 the sequences S(n-1, k) with offset n = 0 are A128834, A001477, A001906, A001353, A004254, A001109, A004187, A001090, A018913, A004189.
(End)
For more on the Diophantine equation presented by Kul, see the Ismail paper. - Tom Copeland, Jan 31 2016
The o.g.f. for the Legendre polynomials L(n,x) is 1 / sqrt(1- 2x*z + z^2), and squaring it gives the o.g.f. of U(n,x), A053117, so Sum_{k=0..n} L(k,x/2) L(n-k,x/2) = S(n,x). This gives S(n,x) = L(n/2,x/2)^2 + 2*Sum_{k=0..n/2-1} L(k,x/2) L(n-k,x/2) for n even and S(n,x) = 2*Sum_{k=0..(n-1)/2} L(k,x/2) L(n-k,x/2) for odd n. For a connection to elliptic curves and modular forms, see A053117. For the normalized Legendre polynomials, see A100258. For other properties and relations to other polynomials, see Allouche et al. - Tom Copeland, Feb 04 2016
LG(x,h1,h2) = -log(1 - h1*x + h2*x^2) = Sum_{n>0} F(n,-h1,h2,0,..,0) x^n/n is a log series generator of the bivariate row polynomials of A127672 with A127672(0,0) = 0 and where F(n,b1,b2,..,bn) are the Faber polynomials of A263916. Exp(LG(x,h1,h2)) = 1 / (1 - h1*x + h2*x^2 ) is the o.g.f. of the bivariate row polynomials of this entry. - Tom Copeland, Feb 15 2016 (Instances of the bivariate o.g.f. for this entry are on pp. 5 and 18 of Sunada. - Tom Copeland, Jan 18 2021)
For distinct odd primes p and q the Legendre symbol can be written as Legendre(q,p) = Product_{k=1..P} S(q-1, 2*cos(2*Pi*k/p)), with P = (p-1)/2. See the Lemmermeyer reference, eq. (8.1) on p. 236. Using the zeros of S(q-1, x) (see above) one has S(q-1, x) = Product_{l=1..Q} (x^2 - (2*cos(Pi*l/q))^2), with Q = (q-1)/2. Thus S(q-1, 2*cos(2*Pi*k/p)) = ((-4)^Q)*Product_{l=1..Q} (sin^2(2*Pi*k/p) - sin^2(Pi*l/q)) = ((-4)^Q)*Product_{m=1..Q} (sin^2(2*Pi*k/p) - sin^2(2*Pi*m/q)). For the proof of the last equality see a W. Lang comment on the triangle A057059 for n = Q and an obvious function f. This leads to Eisenstein's proof of the quadratic reciprocity law Legendre(q,p) = ((-1)^(P*Q)) * Legendre(p,q), See the Lemmermeyer reference, pp. 236-237. - Wolfdieter Lang, Aug 28 2016
For connections to generalized Fibonacci polynomials, compare their generating function on p. 5 of the Amdeberhan et al. link with the o.g.f. given above for the bivariate row polynomials of this entry. - Tom Copeland, Jan 08 2017
The formula for Ramanujan's tau function (see A000594) for prime powers is tau(p^k) = p^(11*k/2)*S(k, p^(-11/2)*tau(p)) for k >= 1, and p = A000040(n), n >= 1. See the Hardy reference, p. 164, eqs. (10.3.4) and (10.3.6) rewritten in terms of S. - Wolfdieter Lang, Jan 27 2017
From Wolfdieter Lang, May 08 2017: (Start)
The number of zeros Z(n) of the S(n, x) polynomials in the open interval (-1,+1) is 2*b(n) for even n >= 0 and 1 + 2*b(n) for odd n >= 1, where b(n) = floor(n/2) - floor((n+1)/3). This b(n) is the number of integers k in the interval (n+1)/3 < k <= floor(n/2). See a comment on the zeros of S(n, x) above, and b(n) = A008615(n-2), n >= 0. The numbers Z(n) have been proposed (with a conjecture related to A008611) by Michel Lagneau, as the number of zeros of Fibonacci polynomials on the imaginary axis (-I,+I), with I=sqrt(-1). They are Z(n) = A008611(n-1), n >= 0, with A008611(-1) = 0. Also Z(n) = A194960(n-4), n >= 0. Proof using the A008611 version. A194960 follows from this.
In general the number of zeros Z(a;n) of S(n, x) for n >= 0 in the open interval (-a,+a) for a from the interval (0,2) (x >= 2 never has zeros, and a=0 is trivial: Z(0;n) = 0) is with b(a;n) = floor(n//2) - floor((n+1)*arccos(a/2)/Pi), as above Z(a;n) = 2*b(a;n) for even n >= 0 and 1 + 2*b(a;n) for odd n >= 1. For the closed interval [-a,+a] Z(0;n) = 1 and for a from (0,1) one uses for Z(a;n) the values b(a;n) = floor(n/2) - ceiling((n+1)*arccos(a/2)/Pi) + 1. (End)
The Riordan row polynomials S(n, x) (Chebyshev S) belong to the Boas-Buck class (see a comment and references in A046521), hence they satisfy the Boas-Buck identity: (E_x - n*1)*S(n, x) = (E_x + 1)*Sum_{p=0..n-1} (1 - (-1)^p)*(-1)^((p+1)/2)*S(n-1-p, x), for n >= 0, where E_x = x*d/dx (Euler operator). For the triangle T(n, k) this entails a recurrence for the sequence of column k, given in the formula section. - Wolfdieter Lang, Aug 11 2017
The e.g.f. E(x,t) := Sum_{n>=0} (t^n/n!)*S(n,x) for the row polynomials is obtained via inverse Laplace transformation from the above given o.g.f. as E(x,t) = ((1/xm)*exp(t/xm) - (1/xp)*exp(t/xp) )/(xp - xm) with xp = (x + sqrt(x^2-4))/2 and xm = (x - sqrt(x^2-4))/2. - Wolfdieter Lang, Nov 08 2017
From Wolfdieter Lang, Apr 12 2018: (Start)
Factorization of row polynomials S(n, x), for n >= 1, in terms of C polynomials (not Chebyshev C) with coefficients given in A187360. This is obtained from the factorization into Psi polynomials (see the Jul 12 2011 comment above) but written in terms of minimal polynomials of 2*cos(2*Pi/n) with coefficients in A232624:
S(2*k, x) = Product_{2 <= d | (2*k+1)} C(d, x)*(-1)^deg(d)*C(d, -x), with deg(d) = A055034(d) the degree of C(d, x).
S(2*k+1, x) = Product_{2 <= d | 2*(k+1)} C(d, x) * Product_{3 <= 2*d + 1 | (k+1)} (-1)^(deg(2*d+1))*C(2*d+1, -x).
Note that (-1)^(deg(2*d+1))*C(2*d+1, -x)*C(2*d+1, x) pairs always appear.
The number of C factors of S(2*k, x), for k >= 0, is 2*(tau(2*k+1) - 1) = 2*(A099774(k+1) - 1) = 2*A095374(k), and for S(2*k+1, x), for k >= 0, it is tau(2*(k+1)) + tau_{odd}(k+1) - 2 = A302707(k), with tau(2*k+1) = A099774(k+1), tau(n) = A000005 and tau(2*(k+1)) = A099777(k+1).
For the reverse problem, the factorization of C polynomials into S polynomials, see A255237. (End)
The S polynomials with general initial conditions S(a,b;n,x) = x*S(a,b;n-1,x) - S(a,b;n-2,x), for n >= 1, with S(a,b;-1,x) = a and S(a,b;0,x) = b are S(a,b;n,x) = b*S(n, x) - a*S(n-1, x), for n >= -1. Recall that S(-2, x) = -1 and S(-1, x) = 0. The o.g.f. is G(a,b;z,x) = (b - a*z)/(1 - x*z + z^2). - Wolfdieter Lang, Oct 18 2019
Also the convolution triangle of A101455. - Peter Luschny, Oct 06 2022
From Wolfdieter Lang, Apr 26 2023: (Start)
Multi-section of S-polynomials: S(m*n+k, x) = S(m+k, x)*S(n-1, R(m, x)) - S(k, x)*S(n-2, R(m, x)), with R(n, x) = S(n, x) - S(n-2, x) (see A127672), S(-2, x) = -1, and S(-1, x) = 0, for n >= 0, m >= 1, and k = 0, 1, ..., m-1.
O.g.f. of {S(m*n+k, y)}_{n>=0}: G(m,k,y,x) = (S(k, y) - (S(k, y)*R(m, y) - S(m+k, y))*x)/(1 - R(m,y)*x + x^2).
See eqs. (40) and (49), with r = x or y and s =-1, of the G. Detlefs and W. Lang link at A034807. (End)
S(n, x) for complex n and complex x: S(n, x) = ((-i/2)/sqrt(1 - (x/2)^2))*(q(x/2)*exp(+n*log(q(x/2))) - (1/q(x/2))*exp(-n*log(q(x/2)))), with q(x) = x + sqrt(1 - x^2)*i. Here log(z) = |z| + Arg(z)*i, with Arg(z) from [-Pi,+Pi) (principal branch). This satisfies the recurrence relation for S because it is derived from the Binet - de Moivre formula for S. Examples: S(n/m, 0) = cos((n/m)*Pi/4), for n >= 0 and m >= 1. S(n*i, 0) = (1/2)*(1 + exp(n*Pi))*exp(-(n/2)*Pi), for n >= 0. S(1+i, 2+i) = 0.6397424847... + 1.0355669490...*i. Thanks to Roberto Alfano for asking a question leading to this formula. - Wolfdieter Lang, Jun 05 2023
Lim_{n->oo} S(n, x)/S(n-1, x) = r(x) = (x - sqrt(x^2 -4))/2, for |x| >= 2. For x = +-2, this limit is +-1. - Wolfdieter Lang, Nov 15 2023

Examples

			The triangle T(n, k) begins:
  n\k  0  1   2   3   4   5   6    7   8   9  10  11
  0:   1
  1:   0  1
  2:  -1  0   1
  3:   0 -2   0   1
  4:   1  0  -3   0   1
  5:   0  3   0  -4   0   1
  6:  -1  0   6   0  -5   0   1
  7:   0 -4   0  10   0  -6   0    1
  8:   1  0 -10   0  15   0  -7    0   1
  9:   0  5   0 -20   0  21   0   -8   0   1
  10: -1  0  15   0 -35   0  28    0  -9   0   1
  11:  0 -6   0  35   0 -56   0   36   0 -10   0   1
  ... Reformatted and extended by _Wolfdieter Lang_, Oct 24 2012
For more rows see the link.
E.g., fourth row {0,-2,0,1} corresponds to polynomial S(3,x)= -2*x + x^3.
From _Wolfdieter Lang_, Jul 12 2011: (Start)
Zeros of S(3,x) with rho(4)= 2*cos(Pi/4) = sqrt(2):
  +- t(1,sqrt(2)) = +- sqrt(2) and
  +- t(2,sqrt(2)) = +- 0.
Factorization of S(3,x) in terms of Psi polynomials:
S(3,x) = (2^3)*Psi(4,x/2)*Psi(8,x/2) = x*(x^2-2).
(End)
From _Wolfdieter Lang_, Nov 04 2011: (Start)
A- and Z- sequence recurrence:
T(4,0) = - (C(0)*T(3,1) + C(1)*T(3,3)) = -(-2 + 1) = +1,
T(5,3) = -3 - 1*1 = -4.
(End)
Boas-Buck recurrence for column k = 2, n = 6: S(6, 2) = (3/4)*(0 - 2* S(4 ,2) + 0 + 2*S(2, 2)) = (3/4)*(-2*(-3) + 2) = 6. - _Wolfdieter Lang_, Aug 11 2017
From _Wolfdieter Lang_, Apr 12 2018: (Start)
Factorization into C polynomials (see the Apr 12 2018 comment):
S(4, x) = 1 - 3*x^2 + x^4 = (-1 + x + x^2)*(-1 - x + x^2) = (-C(5, -x)) * C(5, x); the number of factors is 2 = 2*A095374(2).
S(5, x) = 3*x - 4*x^3 + x^5 = x*(-1 + x)*(1 + x)*(-3 + x^2) = C(2, x)*C(3, x)*(-C(3, -x))*C(6, x); the number of factors is 4 = A302707(2). (End)
		

References

  • G. H. Hardy, Ramanujan: twelve lectures on subjects suggested by his life and work, AMS Chelsea Publishing, Providence, Rhode Island, 2002, p. 164.
  • Max Koecher and Aloys Krieg, Elliptische Funktionen und Modulformen, 2. Auflage, Springer, 2007, p. 223.
  • Franz Lemmermeyer, Reciprocity Laws. From Euler to Eisenstein, Springer, 2000.
  • D. S. Mitrinovic, Analytic Inequalities, Springer-Verlag, 1970; p. 232, Sect. 3.3.38.
  • Theodore J. Rivlin, Chebyshev polynomials: from approximation theory to algebra and number theory, 2. ed., Wiley, New York, 1990, pp. 60 - 61.
  • R. Vein and P. Dale, Determinants and Their Applications in Mathematical Physics, Springer, 1999.

Crossrefs

Cf. A000005, A000217, A000292, A000332, A000389, A001227, A007318, A008611, A008615, A101455, A010892, A011973, A053112 (without zeros), A053117, A053119 (reflection), A053121 (inverse triangle), A055034, A097610, A099774, A099777, A100258, A112552 (first column clipped), A127672, A168561 (absolute values), A187360. A194960, A232624, A255237.
Triangles of coefficients of Chebyshev's S(n,x+k) for k = 5, 4, 3, 2, 1, 0, -1, -2, -3, -4, -5: A207824, A207823, A125662, A078812, A101950, A049310, A104562, A053122, A207815, A159764, A123967.

Programs

  • Magma
    A049310:= func< n,k | ((n+k) mod 2) eq 0 select (-1)^(Floor((n+k)/2)+k)*Binomial(Floor((n+k)/2), k) else 0 >;
    [A049310(n,k): k in [0..n], n in [0..15]]; // G. C. Greubel, Jul 25 2022
  • Maple
    A049310 := proc(n,k): binomial((n+k)/2,(n-k)/2)*cos(Pi*(n-k)/2)*(1+(-1)^(n-k))/2 end: seq(seq(A049310(n,k), k=0..n),n=0..11); # Johannes W. Meijer, Aug 08 2011
    # Uses function PMatrix from A357368. Adds a row above and a column to the left.
    PMatrix(10, n -> ifelse(irem(n, 2) = 0, 0, (-1)^iquo(n-1, 2))); # Peter Luschny, Oct 06 2022
  • Mathematica
    t[n_, k_] /; EvenQ[n+k] = ((-1)^((n+k)/2+k))*Binomial[(n+k)/2, k]; t[n_, k_] /; OddQ[n+k] = 0; Flatten[Table[t[n, k], {n, 0, 12}, {k, 0, n}]][[;; 86]] (* Jean-François Alcover, Jul 05 2011 *)
    Table[Coefficient[(-I)^n Fibonacci[n + 1, - I x], x, k], {n, 0, 10}, {k, 0, n}] //Flatten (* Clark Kimberling, Aug 02 2011; corrected by Eric W. Weisstein, Apr 06 2017 *)
    CoefficientList[ChebyshevU[Range[0, 10], -x/2], x] // Flatten (* Eric W. Weisstein, Apr 06 2017 *)
    CoefficientList[Table[(-I)^n Fibonacci[n + 1, -I x], {n, 0, 10}], x] // Flatten (* Eric W. Weisstein, Apr 06 2017 *)
  • PARI
    {T(n, k) = if( k<0 || k>n || (n + k)%2, 0, (-1)^((n + k)/2 + k) * binomial((n + k)/2, k))} /* Michael Somos, Jun 24 2002 */
    
  • SageMath
    @CachedFunction
    def A049310(n,k):
        if n< 0: return 0
        if n==0: return 1 if k == 0 else 0
        return A049310(n-1,k-1) - A049310(n-2,k)
    for n in (0..9): [A049310(n,k) for k in (0..n)] # Peter Luschny, Nov 20 2012
    

Formula

T(n,k) := 0 if n < k or n+k odd, otherwise ((-1)^((n+k)/2+k))*binomial((n+k)/2, k); T(n, k) = -T(n-2, k)+T(n-1, k-1), T(n, -1) := 0 =: T(-1, k), T(0, 0)=1, T(n, k)= 0 if n < k or n+k odd; g.f. k-th column: (1 / (1 + x^2)^(k + 1)) * x^k. - Michael Somos, Jun 24 2002
T(n,k) = binomial((n+k)/2, (n-k)/2)*cos(Pi*(n-k)/2)*(1+(-1)^(n-k))/2. - Paul Barry, Aug 28 2005
Sum_{k=0..n} T(n,k)^2 = A051286(n). - Philippe Deléham, Nov 21 2005
Recurrence for the (unsigned) Fibonacci polynomials: F(1)=1, F(2)=x; for n > 2, F(n) = x*F(n-1) + F(n-2).
From Wolfdieter Lang, Nov 04 2011: (Start)
The Riordan A- and Z-sequences, given in a comment above, lead together to the recurrence:
T(n,k) = 0 if n < k, if k=0 then T(0,0)=1 and
T(n,0)= -Sum_{i=0..floor((n-1)/2)} C(i)*T(n-1,2*i+1), otherwise T(n,k) = T(n-1,k-1) - Sum_{i=1..floor((n-k)/2)} C(i)*T(n-1,k-1+2*i), with the Catalan numbers C(n)=A000108(n).
(End)
The row polynomials satisfy also S(n,x) = 2*(T(n+2, x/2) - T(n, x/2))/(x^2-4) with the Chebyshev T-polynomials. Proof: Use the trace formula 2*T(n, x/2) = S(n, x) - S(n-2, x) (see the Dec 02 2010 comment above) and the S-recurrence several times. This is a formula which expresses the S- in terms of the T-polynomials. - Wolfdieter Lang, Aug 07 2014
From Tom Copeland, Dec 06 2015: (Start)
The non-vanishing, unsigned subdiagonals Diag_(2n) contain the elements D(n,k) = Sum_{j=0..k} D(n-1,j) = (k+1) (k+2) ... (k+n) / n! = binomial(n+k,n), so the o.g.f. for the subdiagonal is (1-x)^(-(n+1)). E.g., Diag_4 contains D(2,3) = D(1,0) + D(1,1) + D(1,2) + D(1,3) = 1 + 2 + 3 + 4 = 10 = binomial(5,2). Diag_4 is shifted A000217; Diag_6, shifted A000292: Diag_8, shifted A000332; and Diag_10, A000389.
The non-vanishing antidiagonals are signed rows of the Pascal triangle A007318.
For a reversed, unsigned version with the zeros removed, see A011973. (End)
The Boas-Buck recurrence (see a comment above) for the sequence of column k is: S(n, k) = ((k+1)/(n-k))*Sum_{p=0..n-1-k} (1 - (-1)^p)*(-1)^((p+1)/2) * S(n-1-p, k), for n > k >= 0 and input S(k, k) = 1. - Wolfdieter Lang, Aug 11 2017
The m-th row consecutive nonzero entries in order are (-1)^c*(c+b)!/c!b! with c = m/2, m/2-1, ..., 0 and b = m-2c if m is even and with c = (m-1)/2, (m-1)/2-1, ..., 0 with b = m-2c if m is odd. For the 8th row starting at a(36) the 5 consecutive nonzero entries in order are 1,-10,15,-7,1 given by c = 4,3,2,1,0 and b = 0,2,4,6,8. - Richard Turk, Aug 20 2017
O.g.f.: exp( Sum_{n >= 0} 2*T(n,x/2)*t^n/n ) = 1 + x*t + (-1 + x^2)*t^2 + (-2*x + x^3)*t^3 + (1 - 3*x^2 + x^4)*t^4 + ..., where T(n,x) denotes the n-th Chebyshev polynomial of the first kind. - Peter Bala, Aug 15 2022

A008836 Liouville's function lambda(n) = (-1)^k, where k is number of primes dividing n (counted with multiplicity).

Original entry on oeis.org

1, -1, -1, 1, -1, 1, -1, -1, 1, 1, -1, -1, -1, 1, 1, 1, -1, -1, -1, -1, 1, 1, -1, 1, 1, 1, -1, -1, -1, -1, -1, -1, 1, 1, 1, 1, -1, 1, 1, 1, -1, -1, -1, -1, -1, 1, -1, -1, 1, -1, 1, -1, -1, 1, 1, 1, 1, 1, -1, 1, -1, 1, -1, 1, 1, -1, -1, -1, 1, -1, -1, -1, -1, 1, -1, -1, 1, -1, -1, -1, 1, 1, -1, 1, 1, 1, 1, 1, -1, 1, 1, -1, 1, 1, 1, 1, -1, -1, -1, 1, -1
Offset: 1

Views

Author

Keywords

Comments

Coons and Borwein: "We give a new proof of Fatou's theorem: if an algebraic function has a power series expansion with bounded integer coefficients, then it must be a rational function. This result is applied to show that for any non-trivial completely multiplicative function from N to {-1,1}, the series sum_{n=1..infinity} f(n)z^n is transcendental over {Z}[z]; in particular, sum_{n=1..infinity} lambda(n)z^n is transcendental, where lambda is Liouville's function. The transcendence of sum_{n=1..infinity} mu(n)z^n is also proved." - Jonathan Vos Post, Jun 11 2008
Coons proves that a(n) is not k-automatic for any k > 2. - Jonathan Vos Post, Oct 22 2008
The Riemann hypothesis is equivalent to the statement that for every fixed epsilon > 0, lim_{n -> infinity} (a(1) + a(2) + ... + a(n))/n^(1/2 + epsilon) = 0 (Borwein et al., theorem 1.2). - Arkadiusz Wesolowski, Oct 08 2013

Examples

			a(4) = 1 because since bigomega(4) = 2 (the prime divisor 2 is counted twice), then (-1)^2 = 1.
a(5) = -1 because 5 is prime and therefore bigomega(5) = 1 and (-1)^1 = -1.
		

References

  • T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 37.
  • P. Borwein, S. Choi, B. Rooney and A. Weirathmueller, The Riemann Hypothesis: A Resource for the Aficionado and Virtuoso Alike, Springer, Berlin, 2008, pp. 1-11.
  • H. Gupta, On a table of values of L(n), Proceedings of the Indian Academy of Sciences. Section A, 12 (1940), 407-409.
  • H. Gupta, A table of values of Liouville's function L(n), Research Bulletin of East Panjab University, No. 3 (Feb. 1950), 45-55.
  • P. Ribenboim, Algebraic Numbers, p. 44.
  • J. Roberts, Lure of the Integers, Math. Assoc. America, 1992, p. 279.
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, Exercise 3.3.5 on page 99.
  • J. V. Uspensky and M. A. Heaslet, Elementary Number Theory, McGraw-Hill, NY, 1939, p. 112.

Crossrefs

Möbius transform of A010052.
Cf. A182448 (Dgf at s=2), A347328 (Dgf at s=3), A347329 (Dgf at s=4).

Programs

  • Haskell
    a008836 = (1 -) . (* 2) . a066829  -- Reinhard Zumkeller, Nov 19 2011
    
  • Maple
    A008836 := n -> (-1)^numtheory[bigomega](n); # Peter Luschny, Sep 15 2011
    with(numtheory): A008836 := proc(n) local i,it,s; it := ifactors(n): s := (-1)^add(it[2][i][2], i=1..nops(it[2])): RETURN(s) end:
  • Mathematica
    Table[LiouvilleLambda[n], {n, 100}] (* Enrique Pérez Herrero, Dec 28 2009 *)
    Table[If[OddQ[PrimeOmega[n]],-1,1],{n,110}] (* Harvey P. Dale, Sep 10 2014 *)
  • PARI
    {a(n) = if( n<1, 0, n=factor(n); (-1)^sum(i=1, matsize(n)[1], n[i,2]))}; /* Michael Somos, Jan 01 2006 */
    
  • PARI
    a(n)=(-1)^bigomega(n) \\ Charles R Greathouse IV, Jan 09 2013
    
  • Python
    from sympy import factorint
    def A008836(n): return -1 if sum(factorint(n).values()) % 2 else 1 # Chai Wah Wu, May 24 2022

Formula

Dirichlet g.f.: zeta(2s)/zeta(s); Dirichlet inverse of A008966.
Sum_{ d divides n } lambda(d) = 1 if n is a square, otherwise 0.
Completely multiplicative with a(p) = -1, p prime.
a(n) = (-1)^A001222(n) = (-1)^bigomega(n). - Jonathan Vos Post, Apr 16 2006
a(n) = A033999(A001222(n)). - Jaroslav Krizek, Sep 28 2009
Sum_{d|n} a(d) *(A000005(d))^2 = a(n) *Sum{d|n} A000005(d). - Vladimir Shevelev, May 22 2010
a(n) = 1 - 2*A066829(n). - Reinhard Zumkeller, Nov 19 2011
a(n) = i^(tau(n^2)-1) where tau(n) is A000005 and i is the imaginary unit. - Anthony Browne, May 11 2016
a(n) = A106400(A156552(n)). - Antti Karttunen, May 30 2017
Recurrence: a(1)=1, n > 1: a(n) = sign(1/2 - Sum_{dMats Granvik, Oct 11 2017
a(n) = Sum_{ d | n } A008683(d)*A010052(n/d). - Jinyuan Wang, Apr 20 2019
a(1) = 1; a(n) = -Sum_{d|n, d < n} mu(n/d)^2 * a(d). - Ilya Gutkovskiy, Mar 10 2021
a(n) = (-1)^A349905(n). - Antti Karttunen, Apr 26 2022
From Ridouane Oudra, Jun 02 2024: (Start)
a(n) = (-1)^A066829(n);
a(n) = (-1)^A063647(n);
a(n) = A101455(A048691(n));
a(n) = sin(tau(n^2)*Pi/2). (End)

A056594 Period 4: repeat [1,0,-1,0]; expansion of 1/(1 + x^2).

Original entry on oeis.org

1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0
Offset: 0

Views

Author

Wolfdieter Lang, Aug 04 2000

Keywords

Comments

G.f. is inverse of cyclotomic(4,x). Unsigned: A000035(n+1).
Real part of i^n and imaginary part of i^(n+1), i=sqrt(-1). - Reinhard Zumkeller, Jul 22 2007
The BINOMIAL transform generates A009116(n); the inverse BINOMIAL transform generates (-1)^n*A009116(n). - R. J. Mathar, Apr 07 2008
a(n-1), n >= 1, is the nontrivial Dirichlet character modulo 4, called Chi_2(4;n) (the trivial one is Chi_1(4;n) given by periodic(1,0) = A000035(n)). See the Apostol reference, p. 139, the k = 4, phi(k) = 2 table. - Wolfdieter Lang, Jun 21 2011
a(n-1), n >= 1, is the character of the Dirichlet beta function. - Daniel Forgues, Sep 15 2012
a(n-1), n >= 1, is also the (strongly) multiplicative function h(n) of Theorem 5.12, p. 150, of the Niven-Zuckerman reference. See the formula section. This function h(n) can be employed to count the integer solutions to n = x^2 + y^2. See A002654 for a comment with the formula. - Wolfdieter Lang, Apr 19 2013
This sequence is duplicated in A101455 but with offset 1. - Gary Detlefs, Oct 04 2013
For n >= 2 this gives the determinant of the bipartite graph with 2*n nodes and the adjacency matrix A(n) with elements A(n;1,2) = 1 = A(n;n,n-1), and for 1 < i < n A(n;i,i+1) = 1 = A(n;i,i-1), otherwise 0. - Wolfdieter Lang, Jun 25 2023

Examples

			With a(n-1) = h(n) of Niven-Zuckerman: a(62) = h(63) = h(3^2*7^1) = (-1)^(2*1)*(-1)^(1*3) = -1 = h(3)^2*h(7) = a(2)^2*a(6) = (-1)^2*(-1) = -1. - _Wolfdieter Lang_, Apr 19 2013
		

References

  • T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1986.
  • I. S. Gradstein and I. M. Ryshik, Tables of series, products, and integrals, Volume 1, Verlag Harri Deutsch, 1981.
  • Ivan Niven and Herbert S. Zuckerman, An Introduction to the Theory of Numbers, New York: John Wiley (1980), p. 150.
  • Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 32, equation 32:6:1 at page 300.

Crossrefs

Cf. A049310, A074661, A131852, A002654, A146559 (binomial transform).

Programs

  • Magma
    &cat[ [1, 0, -1, 0]: n in [0..23] ]; // Bruno Berselli, Feb 08 2011
    
  • Maple
    A056594 := n->(1-irem(n,2))*(-1)^iquo(n,2); # Peter Luschny, Jul 27 2011
  • Mathematica
    CoefficientList[Series[1/(1 + x^2), {x, 0, 50}], x]
    a[n_]:= KroneckerSymbol[-4,n+1];Table[a[n],{n,0,93}] (* Thanks to Jean-François Alcover. - Wolfdieter Lang, May 31 2013 *)
    CoefficientList[Series[1/Cyclotomic[4, x], {x, 0, 100}], x] (* Vincenzo Librandi, Apr 03 2014 *)
  • Maxima
    A056594(n) := block(
            [1,0,-1,0][1+mod(n,4)]
    )$ /* R. J. Mathar, Mar 19 2012 */
    
  • PARI
    {a(n) = real( I^n )}
    
  • PARI
    {a(n) = kronecker(-4, n+1) }
    
  • Python
    def A056594(n): return (1,0,-1,0)[n&3] # Chai Wah Wu, Sep 23 2023

Formula

G.f.: 1/(1+x^2).
E.g.f.: cos(x).
a(n) = (1/2)*((-i)^n + i^n), where i = sqrt(-1). - Mitch Harris, Apr 19 2005
a(n) = (1/2)*((-1)^(n+floor(n/2)) + (-1)^floor(n/2)).
Recurrence: a(n)=a(n-4), a(0)=1, a(1)=0, a(2)=-1, a(3)=0.
a(n) = T(n, 0) = A053120(n, 0); T(n, x) Chebyshev polynomials of the first kind. - Wolfdieter Lang, Aug 21 2009
a(n) = S(n, 0) = A049310(n, 0); S(n, x) := U(n, x/2), Chebyshev polynomials of 2nd kind.
Sum_{k>=0} a(k)/(k+1) = Pi/4. - Jaume Oliver Lafont, Mar 30 2010
a(n) = Sum_{k=0..n} A101950(n,k)*(-1)^k. - Philippe Deléham, Feb 10 2012
a(n) = (1/2)*(1 + (-1)^n)*(-1)^(n/2). - Bruno Berselli, Mar 13 2012
a(0) = 1, a(n-1) = 0 if n is even, a(n-1) = Product_{j=1..m} (-1)^(e_j*(p_j-1)/2) if the odd n-1 = p_1^(e_1)*p_2^(e_2)*...*p_m^(e_m) with distinct odd primes p_j, j=1..m. See the function h(n) of Theorem 5.12 of the Niven-Zuckerman reference. - Wolfdieter Lang, Apr 19 2013
a(n) = (-4/(n+1)), n >= 0, where (k/n) is the Kronecker symbol. See the Eric Weisstein and Wikipedia links. Thanks to Wesley Ivan Hurt. - Wolfdieter Lang, May 31 2013
a(n) = R(n,0)/2 with the row polynomials R of A127672. This follows from the product of the zeros of R, and the formula Product_{k=0..n-1} 2*cos((2*k+1)*Pi/(2*n)) = (1 + (-1)^n)*(-1)^(n/2), n >= 1 (see the Gradstein and Ryshik reference, p. 63, 1.396 4., with x = sqrt(-1)). - Wolfdieter Lang, Oct 21 2013
a(n) = Sum_{k=0..n} i^(k*(k+1)), where i=sqrt(-1). - Bruno Berselli, Mar 11 2015
Dirichlet g.f. of a(n) shifted right: L(chi_2(4),s) = beta(s) = (1-2^(-s))*(d.g.f. of A034947), see comments by Lang and Forgues. - Ralf Stephan, Mar 27 2015
a(n) = cos(n*Pi/2). - Ridouane Oudra, Sep 29 2024

A003881 Decimal expansion of Pi/4.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

Also the ratio of the area of a circle to the circumscribed square. More generally, the ratio of the area of an ellipse to the circumscribed rectangle. Also the ratio of the volume of a cylinder to the circumscribed cube. - Omar E. Pol, Sep 25 2013
Also the surface area of a quarter-sphere of diameter 1. - Omar E. Pol, Oct 03 2013
Least positive solution to sin(x) = cos(x). - Franklin T. Adams-Watters, Jun 17 2014
Dirichlet L-series of the non-principal character modulo 4 (A101455) at 1. See e.g. Table 22 of arXiv:1008.2547. - R. J. Mathar, May 27 2016
This constant is also equal to the infinite sum of the arctangent functions with nested radicals consisting of square roots of two. Specifically, one of the Viete-like formulas for Pi is given by Pi/4 = Sum_{k = 2..oo} arctan(sqrt(2 - a_{k - 1})/a_k), where the nested radicals are defined by recurrence relations a_k = sqrt(2 + a_{k - 1}) and a_1 = sqrt(2) (see the article [Abrarov and Quine]). - Sanjar Abrarov, Jan 09 2017
Pi/4 is the area enclosed between circumcircle and incircle of a regular polygon of unit side. - Mohammed Yaseen, Nov 29 2023

Examples

			0.785398163397448309615660845819875721049292349843776455243736148...
N = 2, m = 6: Pi/4 = 4!*3^4 Sum_{k >= 0} (-1)^k/((2*k - 11)*(2*k - 5)*(2*k + 1)*(2*k + 7)*(2*k + 13)). - _Peter Bala_, Nov 15 2016
		

References

  • Jörg Arndt and Christoph Haenel, Pi: Algorithmen, Computer, Arithmetik, Springer 2000, p. 150.
  • Florian Cajori, A History of Mathematical Notations, Dover edition (2012), par. 437.
  • Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, 2003, Sections 6.3 and 8.4, pp. 429 and 492.
  • Douglas R. Hofstadter, Gödel, Escher, Bach: An Eternal Golden Braid, Basic Books, p. 408.
  • J. Rivaud, Analyse, Séries, équations différentielles, Mathématiques supérieures et spéciales, Premier cycle universitaire, Vuibert, 1981, Exercice 3, p. 136.
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 119.

Crossrefs

Cf. A006752 (beta(2)=Catalan), A153071 (beta(3)), A175572 (beta(4)), A175571 (beta(5)), A175570 (beta(6)), A258814 (beta(7)), A258815 (beta(8)), A258816 (beta(9)).
Cf. A001622.

Programs

  • Haskell
    -- see link: Literate Programs
    import Data.Char (digitToInt)
    a003881_list len = map digitToInt $ show $ machin `div` (10 ^ 10) where
       machin = 4 * arccot 5 unity - arccot 239 unity
       unity = 10 ^ (len + 10)
       arccot x unity = arccot' x unity 0 (unity `div` x) 1 1 where
         arccot' x unity summa xpow n sign
        | term == 0 = summa
        | otherwise = arccot'
          x unity (summa + sign * term) (xpow `div` x ^ 2) (n + 2) (- sign)
        where term = xpow `div` n
    -- Reinhard Zumkeller, Nov 20 2012
    
  • Magma
    R:= RealField(100); Pi(R)/4; // G. C. Greubel, Mar 08 2018
  • Maple
    evalf(Pi/4) ;
  • Mathematica
    RealDigits[N[Pi/4,6! ]]  (* Vladimir Joseph Stephan Orlovsky, Dec 02 2009 *)
    (* PROGRAM STARTS *)
    (* Define the nested radicals a_k by recurrence *)
    a[k_] := Nest[Sqrt[2 + #1] & , 0, k]
    (* Example of Pi/4 approximation at K = 100 *)
    Print["The actual value of Pi/4 is"]
    N[Pi/4, 40]
    Print["At K = 100 the approximated value of Pi/4 is"]
    K := 100;  (* the truncating integer *)
    N[Sum[ArcTan[Sqrt[2 - a[k - 1]]/a[k]], {k, 2, K}], 40] (* equation (8) *)
    (* Error terms for Pi/4 approximations *)
    Print["Error terms for Pi/4"]
    k := 1; (* initial value of the index k *)
    K := 10; (* initial value of the truncating integer K *)
    sqn := {}; (* initiate the sequence *)
    AppendTo[sqn, {"Truncating integer K ", " Error term in Pi/4"}];
    While[K <= 30,
    AppendTo[sqn, {K,
       N[Pi/4 - Sum[ArcTan[Sqrt[2 - a[k - 1]]/a[k]], {k, 2, K}], 1000] //
        N}]; K++]
    Print[MatrixForm[sqn]]
    (* Sanjar Abrarov, Jan 09 2017 *)
  • PARI
    Pi/4 \\ Charles R Greathouse IV, Jul 07 2014
    
  • SageMath
    # Leibniz/Cohen/Villegas/Zagier/Arndt/Haenel
    def FastLeibniz(n):
        b = 2^(2*n-1); c = b; s = 0
        for k in range(n-1,-1,-1):
            t = 2*k+1
            s = s + c/t if is_even(k) else s - c/t
            b *= (t*(k+1))/(2*(n-k)*(n+k))
            c += b
        return s/c
    A003881 = RealField(3333)(FastLeibniz(1330))
    print(A003881)  # Peter Luschny, Nov 20 2012
    

Formula

Equals Integral_{x=0..oo} sin(2x)/(2x) dx.
Equals lim_{n->oo} n*A001586(n-1)/A001586(n) (conjecture). - Mats Granvik, Feb 23 2011
Equals Integral_{x=0..1} 1/(1+x^2) dx. - Gary W. Adamson, Jun 22 2003
Equals Integral_{x=0..Pi/2} sin(x)^2 dx, or Integral_{x=0..Pi/2} cos(x)^2 dx. - Jean-François Alcover, Mar 26 2013
Equals (Sum_{x=0..oo} sin(x)*cos(x)/x) - 1/2. - Bruno Berselli, May 13 2013
Equals (-digamma(1/4) + digamma(3/4))/4. - Jean-François Alcover, May 31 2013
Equals Sum_{n>=0} (-1)^n/(2*n+1). - Geoffrey Critzer, Nov 03 2013
Equals Integral_{x=0..1} Product_{k>=1} (1-x^(8*k))^3 dx [cf. A258414]. - Vaclav Kotesovec, May 30 2015
Equals Product_{k in A071904} (if k mod 4 = 1 then (k-1)/(k+1)) else (if k mod 4 = 3 then (k+1)/(k-1)). - Dimitris Valianatos, Oct 05 2016
From Peter Bala, Nov 15 2016: (Start)
For N even: 2*(Pi/4 - Sum_{k = 1..N/2} (-1)^(k-1)/(2*k - 1)) ~ (-1)^(N/2)*(1/N - 1/N^3 + 5/N^5 - 61/N^7 + 1385/N^9 - ...), where the sequence of unsigned coefficients [1, 1, 5, 61, 1385, ...] is A000364. See Borwein et al., Theorem 1 (a).
For N odd: 2*(Pi/4 - Sum_{k = 1..(N-1)/2} (-1)^(k-1)/(2*k - 1)) ~ (-1)^((N-1)/2)*(1/N - 1/N^2 + 2/N^4 - 16/N^6 + 272/N^8 - ...), where the sequence of unsigned coefficients [1, 1, 2, 16, 272, ...] is A000182 with an extra initial term of 1.
For N = 0,1,2,... and m = 1,3,5,... there holds Pi/4 = (2*N)! * m^(2*N) * Sum_{k >= 0} ( (-1)^(N+k) * 1/Product_{j = -N..N} (2*k + 1 + 2*m*j) ); when N = 0 we get the Madhava-Gregory-Leibniz series for Pi/4.
For examples of asymptotic expansions for the tails of these series representations for Pi/4 see A024235 (case N = 1, m = 1), A278080 (case N = 2, m = 1) and A278195 (case N = 3, m = 1).
For N = 0,1,2,..., Pi/4 = 4^(N-1)*N!/(2*N)! * Sum_{k >= 0} 2^(k+1)*(k + N)!* (k + 2*N)!/(2*k + 2*N + 1)!, follows by applying Euler's series transformation to the above series representation for Pi/4 in the case m = 1. (End)
From Peter Bala, Nov 05 2019: (Start)
For k = 0,1,2,..., Pi/4 = k!*Sum_{n = -oo..oo} 1/((4*n+1)*(4*n+3)* ...*(4*n+2*k+1)), where Sum_{n = -oo..oo} f(n) is understood as lim_{j -> oo} Sum_{n = -j..j} f(n).
Equals Integral_{x = 0..oo} sin(x)^4/x^2 dx = Sum_{n >= 1} sin(n)^4/n^2, by the Abel-Plana formula.
Equals Integral_{x = 0..oo} sin(x)^3/x dx = Sum_{n >= 1} sin(n)^3/n, by the Abel-Plana formula. (End)
From Amiram Eldar, Aug 19 2020: (Start)
Equals arcsin(1/sqrt(2)).
Equals Product_{k>=1} (1 - 1/(2*k+1)^2).
Equals Integral_{x=0..oo} x/(x^4 + 1) dx.
Equals Integral_{x=0..oo} 1/(x^2 + 4) dx. (End)
With offset 1, equals 5 * Pi / 2. - Sean A. Irvine, Aug 19 2021
Equals (1/2)!^2 = Gamma(3/2)^2. - Gary W. Adamson, Aug 23 2021
Equals Integral_{x = 0..oo} exp(-x)*sin(x)/x dx (see Rivaud reference). - Bernard Schott, Jan 28 2022
From Amiram Eldar, Nov 06 2023: (Start)
Equals beta(1), where beta is the Dirichlet beta function.
Equals Product_{p prime >= 3} (1 - (-1)^((p-1)/2)/p)^(-1). (End)
Equals arctan( F(1)/F(4) ) + arctan( F(2)/F(3) ), where F(1), F(2), F(3), and F(4) are any four consecutive Fibonacci numbers. - Gary W. Adamson, Mar 03 2024
Pi/4 = Sum_{n >= 1} i/(n*P(n, i)*P(n-1, i)) = (1/2)*Sum_{n >= 1} (-1)^(n+1)*4^n/(n*A006139(n)*A006139(n-1)), where i = sqrt(-1) and P(n, x) denotes the n-th Legendre polynomial. The n-th summand of the series is O( 1/(3 + 2*sqrt(3))^n ). - Peter Bala, Mar 16 2024
Equals arctan( phi^(-3) ) + arctan(phi^(-1) ). - Gary W. Adamson, Mar 27 2024
Equals Sum_{n>=1} eta(n)/2^n, where eta(n) is the Dirichlet eta function. - Antonio Graciá Llorente, Oct 04 2024
Equals Product_{k>=2} ((k + 1)^(k*(2*k + 1))*(k - 1)^(k*(2*k - 1)))/k^(4*k^2). - Antonio Graciá Llorente, Apr 12 2025
Equals Integral_{x=sqrt(2)..oo} dx/(x*sqrt(x^2 - 1)). - Kritsada Moomuang, May 29 2025

Extensions

a(98) and a(99) corrected by Reinhard Zumkeller, Nov 20 2012

A002654 Number of ways of writing n as a sum of at most two nonzero squares, where order matters; also (number of divisors of n of form 4m+1) - (number of divisors of form 4m+3).

Original entry on oeis.org

1, 1, 0, 1, 2, 0, 0, 1, 1, 2, 0, 0, 2, 0, 0, 1, 2, 1, 0, 2, 0, 0, 0, 0, 3, 2, 0, 0, 2, 0, 0, 1, 0, 2, 0, 1, 2, 0, 0, 2, 2, 0, 0, 0, 2, 0, 0, 0, 1, 3, 0, 2, 2, 0, 0, 0, 0, 2, 0, 0, 2, 0, 0, 1, 4, 0, 0, 2, 0, 0, 0, 1, 2, 2, 0, 0, 0, 0, 0, 2, 1, 2, 0, 0, 4, 0, 0, 0, 2, 2, 0, 0, 0, 0, 0, 0, 2, 1, 0, 3, 2, 0, 0, 2, 0
Offset: 1

Views

Author

Keywords

Comments

Glaisher calls this E(n) or E_0(n). - N. J. A. Sloane, Nov 24 2018
Number of sublattices of Z X Z of index n that are similar to Z X Z; number of (principal) ideals of Z[i] of norm n.
a(n) is also one fourth of the number of integer solutions of n = x^2 + y^2 (order and signs matter, and 0 (without signs) is allowed). a(n) = N(n)/4, with N(n) from p. 147 of the Niven-Zuckermann reference. See also Theorem 5.12, p. 150, which defines a (strongly) multiplicative function h(n) which coincides with A056594(n-1), n >= 1, and N(n)/4 = sum(h(d), d divides n). - Wolfdieter Lang, Apr 19 2013
a(2+8*N) = A008441(N) gives the number of ways of writing N as the sum of 2 (nonnegative) triangular numbers for N >= 0. - Wolfdieter Lang, Jan 12 2017
Coefficients of Dedekind zeta function for the quadratic number field of discriminant -4. See A002324 for formula and Maple code. - N. J. A. Sloane, Mar 22 2022

Examples

			4 = 2^2, so a(4) = 1; 5 = 1^2 + 2^2 = 2^2 + 1^2, so a(5) = 2.
x + x^2 + x^4 + 2*x^5 + x^8 + x^9 + 2*x^10 + 2*x^13 + x^16 + 2*x^17 + x^18 + ...
2 = (+1)^2 + (+1)^2 = (+1)^2 + (-1)^2  = (-1)^2 + (+1)^2 = (-1)^2 + (-1)^2. Hence there are 4 integer solutions, called N(2) in the Niven-Zuckerman reference, and a(2) = N(2)/4 = 1.  4 = 0^1 + (+2)^2 = (+2)^2 + 0^2 = 0^2 + (-2)^2 = (-2)^2 + 0^2. Hence N(4) = 4 and a(4) = N(4)/4 = 1. N(5) = 8, a(5) = 2. - _Wolfdieter Lang_, Apr 19 2013
		

References

  • J. M. Borwein, D. H. Bailey and R. Girgensohn, Experimentation in Mathematics, A K Peters, Ltd., Natick, MA, 2004. x+357 pp. See p. 194.
  • George Chrystal, Algebra: An elementary text-book for the higher classes of secondary schools and for colleges, 6th ed., Chelsea Publishing Co., New York, 1959, Part II, p. 346 Exercise XXI(17). MR0121327 (22 #12066)
  • Emil Grosswald, Representations of Integers as Sums of Squares. Springer-Verlag, NY, 1985, p. 15.
  • Ivan Niven and Herbert S. Zuckerman, An Introduction to the Theory of Numbers, New York: John Wiley, 1980, pp. 147 and 150.
  • Günter Scheja and Uwe Storch, Lehrbuch der Algebra, Tuebner, 1988, p. 251.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 89.
  • J. V. Uspensky and M. A. Heaslet, Elementary Number Theory, McGraw-Hill, NY, 1939, p. 340.

Crossrefs

Equals 1/4 of A004018. Partial sums give A014200.
Cf. A002175, A008441, A121444, A122856, A122865, A022544, A143574, A000265, A027748, A124010, A025426 (two squares, order does not matter), A120630 (Dirichlet inverse), A101455 (Mobius transform), A000089, A241011.
If one simply reads the table in Glaisher, PLMS 1884, which omits the zero entries, one gets A213408.
Dedekind zeta functions for imaginary quadratic number fields of discriminants -3, -4, -7, -8, -11, -15, -19, -20 are A002324, A002654, A035182, A002325, A035179, A035175, A035171, A035170, respectively.
Dedekind zeta functions for real quadratic number fields of discriminants 5, 8, 12, 13, 17, 21, 24, 28, 29, 33, 37, 40 are A035187, A035185, A035194, A035195, A035199, A035203, A035188, A035210, A035211, A035215, A035219, A035192, respectively.

Programs

  • Haskell
    a002654 n = product $ zipWith f (a027748_row m) (a124010_row m) where
       f p e | p `mod` 4 == 1 = e + 1
             | otherwise      = (e + 1) `mod` 2
       m = a000265 n
    -- Reinhard Zumkeller, Mar 18 2013
    
  • Maple
    with(numtheory):
    A002654 := proc(n)
        local count1, count3, d;
        count1 := 0:
        count3 := 0:
        for d in numtheory[divisors](n) do
            if d mod 4 = 1 then
                count1 := count1+1
            elif d mod 4 = 3 then
                count3 := count3+1
            fi:
        end do:
        count1-count3;
    end proc:
    # second Maple program:
    a:= n-> add(`if`(d::odd, (-1)^((d-1)/2), 0), d=numtheory[divisors](n)):
    seq(a(n), n=1..100);  # Alois P. Heinz, Feb 04 2020
  • Mathematica
    a[n_] := Count[Divisors[n], d_ /; Mod[d, 4] == 1] - Count[Divisors[n], d_ /; Mod[d, 4] == 3]; a/@Range[105] (* Jean-François Alcover, Apr 06 2011, after R. J. Mathar *)
    QP = QPochhammer; CoefficientList[(1/q)*(QP[q^2]^10/(QP[q]*QP[q^4])^4-1)/4 + O[q]^100, q] (* Jean-François Alcover, Nov 24 2015 *)
    f[2, e_] := 1; f[p_, e_] := If[Mod[p, 4] == 1, e + 1, Mod[e + 1, 2]]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 19 2020 *)
    Rest[CoefficientList[Series[EllipticTheta[3, 0, q]^2/4, {q, 0, 100}], q]] (* Vaclav Kotesovec, Mar 10 2023 *)
  • PARI
    direuler(p=2,101,1/(1-X)/(1-kronecker(-4,p)*X))
    
  • PARI
    {a(n) = polcoeff( sum(k=1, n, x^k / (1 + x^(2*k)), x * O(x^n)), n)}
    
  • PARI
    {a(n) = sumdiv( n, d, (d%4==1) - (d%4==3))}
    
  • PARI
    {a(n) = local(A); A = x * O(x^n); polcoeff( eta(x^2 + A)^10 / (eta(x + A) * eta(x^4 + A))^4 / 4, n)} \\ Michael Somos, Jun 03 2005
    
  • PARI
    a(n)=my(f=factor(n>>valuation(n,2))); prod(i=1,#f~, if(f[i,1]%4==1, f[i,2]+1, (f[i,2]+1)%2)) \\ Charles R Greathouse IV, Sep 09 2014
    
  • PARI
    my(B=bnfinit(x^2+1)); vector(100,n,#bnfisintnorm(B,n)) \\ Joerg Arndt, Jun 01 2024
    
  • Python
    from math import prod
    from sympy import factorint
    def A002654(n): return prod(1 if p == 2 else (e+1 if p % 4 == 1 else (e+1) % 2) for p, e in factorint(n).items()) # Chai Wah Wu, May 09 2022

Formula

Dirichlet series: (1-2^(-s))^(-1)*Product (1-p^(-s))^(-2) (p=1 mod 4) * Product (1-p^(-2s))^(-1) (p=3 mod 4) = Dedekind zeta-function of Z[ i ].
Coefficients in expansion of Dirichlet series Product_p (1-(Kronecker(m, p)+1)*p^(-s)+Kronecker(m, p)*p^(-2s))^(-1) for m = -16.
If n=2^k*u*v, where u is product of primes 4m+1, v is product of primes 4m+3, then a(n)=0 unless v is a square, in which case a(n) = number of divisors of u (Jacobi).
Multiplicative with a(p^e) = 1 if p = 2; e+1 if p == 1 (mod 4); (e+1) mod 2 if p == 3 (mod 4). - David W. Wilson, Sep 01 2001
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^4)) where f(u, v, w) = (u - v)^2 - (v - w) * (4*w + 1). - Michael Somos, Jul 19 2004
G.f.: Sum_{n>=1} ((-1)^floor(n/2)*x^((n^2+n)/2)/(1+(-x)^n)). - Vladeta Jovovic, Sep 15 2004
Expansion of (eta(q^2)^10 / (eta(q) * eta(q^4))^4 - 1)/4 in powers of q.
G.f.: Sum_{k>0} x^k / (1 + x^(2*k)) = Sum_{k>0} -(-1)^k * x^(2*k - 1) / (1 - x^(2*k - 1)). - Michael Somos, Aug 17 2005
a(4*n + 3) = a(9*n + 3) = a(9*n + 6) = 0. a(9*n) = a(2*n) = a(n). - Michael Somos, Nov 01 2006
a(4*n + 1) = A008441(n). a(3*n + 1) = A122865(n). a(3*n + 2) = A122856(n). a(12*n + 1) = A002175(n). a(12*n + 5) = 2 * A121444(n). 4 * a(n) = A004018(n) unless n=0.
a(n) = Sum_{k=1..n} A010052(k)*A010052(n-k). a(A022544(n)) = 0; a(A001481(n)) > 0.
- Reinhard Zumkeller, Sep 27 2008
a(n) = A001826(n) - A001842(n). - R. J. Mathar, Mar 23 2011
a(n) = Sum_{d|n} A056594(d-1), n >= 1. See the above comment on A056594(d-1) = h(d) of the Niven-Zuckerman reference. - Wolfdieter Lang, Apr 19 2013
Dirichlet g.f.: zeta(s)*beta(s) = zeta(s)*L(chi_2(4),s). - Ralf Stephan, Mar 27 2015
G.f.: (theta_3(x)^2 - 1)/4, where theta_3() is the Jacobi theta function. - Ilya Gutkovskiy, Apr 17 2018
a(n) = Sum_{ m: m^2|n } A000089(n/m^2). - Andrey Zabolotskiy, May 07 2018
a(n) = A053866(n) + 2 * A025441(n). - Andrey Zabolotskiy, Apr 23 2019
a(n) = Im(Sum_{d|n} i^d). - Ridouane Oudra, Feb 02 2020
a(n) = Sum_{d|n} sin((1/2)*d*Pi). - Ridouane Oudra, Jan 22 2021
Sum_{n>=1} (-1)^n*a(n)/n = Pi*log(2)/4 (Covo, 2010). - Amiram Eldar, Apr 07 2022
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Pi/4 = 0.785398... (A003881). - Amiram Eldar, Oct 11 2022
From Vaclav Kotesovec, Mar 10 2023: (Start)
Sum_{k=1..n} a(k)^2 ~ n * (log(n) + C) / 4, where C = A241011 =
4*gamma - 1 + log(2)/3 - 2*log(Pi) + 8*log(Gamma(3/4)) - 12*Zeta'(2)/Pi^2 = 2.01662154573340811526279685971511542645018417752364748061...
The constant C, published by Ramanujan (1916, formula (22)), 4*gamma - 1 + log(2)/3 - log(Pi) + 4*log(Gamma(3/4)) - 12*Zeta'(2)/Pi^2 = 2.3482276258576... is wrong! (End)

A321836 a(n) = Sum_{d|n, n/d==1 mod 4} d^12 - Sum_{d|n, n/d==3 mod 4} d^12.

Original entry on oeis.org

1, 4096, 531440, 16777216, 244140626, 2176778240, 13841287200, 68719476736, 282429005041, 1000000004096, 3138428376720, 8916083671040, 23298085122482, 56693912371200, 129746094281440, 281474976710656, 582622237229762, 1156829204647936
Offset: 1

Views

Author

N. J. A. Sloane, Nov 24 2018

Keywords

Crossrefs

Cf. A101455.
Cf. A321543 - A321565, A321807 - A321835 for similar sequences.
Glaisher's E'_i (i=0..12): A002654, A050469, A050470, A050471, A050468, A321829, A321830, A321831, A321832, A321833, A321834, A321835, this sequence.

Programs

  • Mathematica
    s[n_,r_] := DivisorSum[n, #^12 &, Mod[n/#,4]==r &]; a[n_] := s[n,1] - s[n,3]; Array[a, 30] (* Amiram Eldar, Nov 26 2018 *)
    s[n_] := If[OddQ[n], (-1)^((n-1)/2), 0]; (* A101455 *)
    f[p_, e_] := (p^(12*e+12) - s[p]^(e+1))/(p^12 - s[p]); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Nov 04 2023 *)
  • PARI
    apply( a(n)=sumdiv(n, d, if(bittest(n\d,0),(2-n\d%4)*d^12)), [1..30]) \\ M. F. Hasler, Nov 26 2018

Formula

G.f.: Sum_{k>=1} k^12*x^k/(1 + x^(2*k)). - Ilya Gutkovskiy, Nov 26 2018
From Amiram Eldar, Nov 04 2023: (Start)
Multiplicative with a(p^e) = (p^(12*e+12) - A101455(p)^(e+1))/(p^12 - A101455(p)).
Sum_{k=1..n} a(k) ~ c * n^13 / 13, where c = beta(13) = 540553*Pi^13/1569592442880 = 0.999999373583..., and beta is the Dirichlet beta function. (End)
a(n) = Sum_{d|n} (n/d)^12*sin(d*Pi/2). - Ridouane Oudra, Sep 27 2024

A091337 a(n) = (2/n), where (k/n) is the Kronecker symbol.

Original entry on oeis.org

0, 1, 0, -1, 0, -1, 0, 1, 0, 1, 0, -1, 0, -1, 0, 1, 0, 1, 0, -1, 0, -1, 0, 1, 0, 1, 0, -1, 0, -1, 0, 1, 0, 1, 0, -1, 0, -1, 0, 1, 0, 1, 0, -1, 0, -1, 0, 1, 0, 1, 0, -1, 0, -1, 0, 1, 0, 1, 0, -1, 0, -1, 0, 1, 0, 1, 0, -1, 0, -1, 0, 1, 0, 1, 0, -1, 0, -1, 0, 1, 0, 1, 0, -1, 0, -1, 0, 1, 0, 1, 0, -1, 0, -1, 0, 1, 0, 1, 0, -1, 0, -1, 0, 1, 0, 1
Offset: 0

Views

Author

Eric W. Weisstein, Dec 30 2003

Keywords

Comments

Sinh(1) in 'reflected factorial' base is 1.01010101010101010101010101010101010101010101... see A073097 for cosh(1). - Robert G. Wilson v, May 04 2005
A non-principal character for the Dirichlet L-series modulo 8, see arXiv:1008.2547 and L-values Sum_{n >= 1} a(n)/n^s in eq (318) by Jolley. - R. J. Mathar, Oct 06 2011 [The other two non-principal characters are A101455 = {(-4/n)} and A188510 = {(-2/n)}. - Jianing Song, Nov 14 2024]
Period 8: repeat [0, 1, 0, -1, 0, -1, 0, 1]. - Wesley Ivan Hurt, Sep 07 2015 [Adapted by Jianing Song, Nov 14 2024 to include a(0) = 0.]
a(n) = (2^(2i+1)/n), where (k/n) is the Kronecker symbol and i >= 0. - A.H.M. Smeets, Jan 23 2018

Examples

			G.f. = x - x^3 - x^5 + x^7 + x^9 - x^11 - x^13 + x^15 + x^17 - x^19 - x^21 + ...
		

References

  • L. B. W. Jolley, Summation of series, Dover (1961).

Crossrefs

Kronecker symbols {(d/n)} where d is a fundamental discriminant with |d| <= 24: A109017 (d=-24), A011586 (d=-23), A289741 (d=-20), A011585 (d=-19), A316569 (d=-15), A011582 (d=-11), A188510 (d=-8), A175629 (d=-7), A101455 (d=-4), A102283 (d=-3), A080891 (d=5), this sequence (d=8), A110161 (d=12), A011583 (d=13), A011584 (d=17), A322829 (d=21), A322796 (d=24).

Programs

  • Magma
    [(n mod 2) * (-1)^((n+1) div 4)  : n in [1..100]]; // Vincenzo Librandi, Oct 31 2014
  • Maple
    A091337:= n -> [0, 1, 0, -1, 0, -1, 0, 1][(n mod 8)+1]: seq(A091337(n), n=1..100); # Wesley Ivan Hurt, Sep 07 2015
  • Mathematica
    KroneckerSymbol[Range[100], 2] (* Alonso del Arte, Oct 30 2014 *)
  • PARI
    {a(n) = (n%2) * (-1)^((n+1)\4)}; /* Michael Somos, Sep 10 2005 */
    
  • PARI
    {a(n) = kronecker( 2, n)}; /* Michael Somos, Sep 10 2005 */
    
  • PARI
    {a(n) = [0, 1, 0, -1, 0, -1, 0, 1][n%8 + 1]}; /* Michael Somos, Jul 17 2009 */
    

Formula

Euler transform of length 8 sequence [0, -1, 0, -1, 0, 0, 0, 1]. - Michael Somos, Jul 17 2009
a(n) is multiplicative with a(2^e) = 0^e, a(p^e) = 1 if p == 1, 7 (mod 8), a(p^e) = (-1)^e if p == 3, 5 (mod 8). - Michael Somos, Jul 17 2009
G.f.: x*(1 - x^2)/(1 + x^4). a(n) = -a(n + 4) = a(-n) for all n in Z. a(2*n) = 0. a(2*n + 1) = A087960(n). - Michael Somos, Apr 10 2011
Transform of Pell numbers A000129 by the Riordan array A102587. - Paul Barry, Jul 14 2005
a(n) = (2/n) = (n/2), Charles R Greathouse IV explained. - Alonso del Arte, Oct 31 2014
a(n) = (1 - (-1)^n)*(-1)^(n/4 - 1/8 - (-1)^n/8 + (-1)^((2*n + 1 - (-1)^n)/4)/4)/2. - Wesley Ivan Hurt, Sep 07 2015
From Jianing Song, Nov 14 2018: (Start)
a(n) = sqrt(2)*sin(Pi*n/2)*sin(Pi*n/4).
E.g.f.: sqrt(2)*cos(x/sqrt(2))*sinh(x/sqrt(2)).
Moebius transform of A035185.
a(n) = A101455(n)*A188510(n). (End)
a(n) = Sum_{i=1..n} (-1)^(i + floor((i-3)/4)). - Wesley Ivan Hurt, Apr 27 2020
Sum_{n>=1} a(n)/n = A196525. Sum_{n>=1} a(n)/n^2 = A328895. Sum_{n>=1} a(n)/n^3 = A329715. Sum_{n>=1} a(n)/n^4 = A346728. - R. J. Mathar, Dec 17 2024

Extensions

a(0) prepended by Jianing Song, Nov 14 2024

A153071 Decimal expansion of L(3, chi4), where L(s, chi4) is the Dirichlet L-function for the non-principal character modulo 4.

Original entry on oeis.org

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

Views

Author

Stuart Clary, Dec 17 2008

Keywords

Examples

			L(3, chi4) = Pi^3/32 = 0.9689461462593693804836348458469186...
		

References

  • Bruce C. Berndt, Ramanujan's Notebooks, Part II, Springer-Verlag, 1989. See page 293, Entry 25 (iii).
  • Leonhard Euler, Introductio in Analysin Infinitorum, First Part, Articles 175, 284 and 287.
  • Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, 2003, Section 1.4.1, p. 20.

Crossrefs

Cf. A003881 (beta(1)=Pi/4), A006752 (beta(2)=Catalan), A175572 (beta(4)), A175571 (beta(5)), A175570 (beta(6)), A258814 (beta(7)), A258815 (beta(8)), A258816 (beta(9)).

Programs

  • Mathematica
    nmax = 1000; First[ RealDigits[Pi^3/32, 10, nmax] ]
  • PARI
    Pi^3/32 \\ Michel Marcus, Aug 15 2018

Formula

chi4(k) = Kronecker(-4, k); chi4(k) is 0, 1, 0, -1 when k reduced modulo 4 is 0, 1, 2, 3, respectively; chi4 is A101455.
Series: L(3, chi4) = Sum_{k>=1} chi4(k) k^{-3} = 1 - 1/3^3 + 1/5^3 - 1/7^3 + 1/9^3 - 1/11^3 + 1/13^3 - 1/15^3 + ...
Series: L(3, chi4) = Sum_{k>=0} tanh((2k+1) Pi/2)/(2k+1)^3. [Ramanujan; see Berndt, page 293]
Closed form: L(3, chi4) = Pi^3/32 = 1/A331095.
Equals Sum_{n>=0} (-1)^n/(2*n+1)^3. - Jean-François Alcover, Mar 29 2013
Equals Product_{k>=3} (1 - tan(Pi/2^k)^4) (Groenman, 1990). - Amiram Eldar, Apr 03 2022
Equals Integral_{x=0..1} arcsinh(x)*arccos(x)/x dx (Kobayashi, 2021). - Amiram Eldar, Jun 23 2023
From Amiram Eldar, Nov 06 2023: (Start)
Equals beta(3), where beta is the Dirichlet beta function.
Equals Product_{p prime >= 3} (1 - (-1)^((p-1)/2)/p^3)^(-1). (End)

Extensions

Offset corrected by R. J. Mathar, Feb 05 2009

A050469 a(n) = Sum_{ d divides n, n/d=1 mod 4} d - Sum_{ d divides n, n/d=3 mod 4} d.

Original entry on oeis.org

1, 2, 2, 4, 6, 4, 6, 8, 7, 12, 10, 8, 14, 12, 12, 16, 18, 14, 18, 24, 12, 20, 22, 16, 31, 28, 20, 24, 30, 24, 30, 32, 20, 36, 36, 28, 38, 36, 28, 48, 42, 24, 42, 40, 42, 44, 46, 32, 43, 62, 36, 56, 54, 40, 60, 48, 36, 60, 58, 48, 62, 60, 42, 64, 84, 40
Offset: 1

Views

Author

N. J. A. Sloane, Dec 23 1999

Keywords

Comments

Multiplicative with a(p^e)=p^e if p=2, (p^(e+1)-1)/(p-1) if p==1 (mod 4), else (p^(e+1)+(-1)^e)/(p+1). - Michael Somos, May 02 2005
Multiplicative because it is the Dirichlet convolution of A000027 = n and A101455 = [1 0 -1 0 1 0 -1 ...], which are both multiplicative. - Christian G. Bower, May 17 2005

Crossrefs

Programs

  • Mathematica
    max = 70; s = Sum[n*x^(n-1)/(1+x^(2*n)), {n, 1, max}] + O[x]^max; CoefficientList[s, x] (* Jean-François Alcover, Dec 02 2015 *)
    f[p_, e_] := Which[p == 2, p^e, Mod[p, 4] == 1, (p^(e + 1) - 1)/(p - 1), Mod[p, 4] == 3, (p^(e + 1) + (-1)^e)/(p + 1)]; a[1] = 1; a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100] (* Amiram Eldar, Nov 06 2022 *)
  • PARI
    a(n)=if(n<1,0,sumdiv(n,d,d*((n/d%4==1)-(n/d%4==3))))
    
  • PARI
    {a(n)=local(A,p,e); if(n<2, n==1, A=factor(n); prod(k=1,matsize(A)[1], if(p=A[k,1], e=A[k,2]; if(p==2, p^e, if(p%4==1, (p^(e+1)-1)/(p-1), (p^(e+1)+(-1)^e)/(p+1)))))) } /* Michael Somos, May 02 2005 */
    
  • PARI
    a(n)=if(n<1,0,polcoeff(sum(k=1,n,k*x^k/(1+x^(2*k)),x*O(x^n)),n))

Formula

G.f.: Sum_{n>=1} n*x^n/(1+x^(2*n)). - Vladeta Jovovic, Oct 16 2002
L.g.f.: Sum_{k>=1} arctan(x^k). - Ilya Gutkovskiy, Dec 16 2019
O.g.f.: Sum_{n >= 1} (-1)^(n+1) * x^(2*n-1)/(1 - x^(2*n-1))^2. - Peter Bala, Jan 04 2021
Sum_{k=1..n} a(k) ~ c * n^2, where c = (1/2) * Product_{primes p == 1 (mod 4)} 1/(1-1/p^2) * Product_{primes p == 3 (mod 4)} 1/(1+1/p^2) = (1/2) * A175647 / A243381 = A006752/2 = 0.4579827970... . - Amiram Eldar, Nov 06 2022, Nov 05 2023
a(n) = Sum_{d|n} (n/d)*sin(d*Pi/2). - Ridouane Oudra, Sep 26 2024

A050470 a(n) = Sum_{d|n, n/d == 1 (mod 4)} d^2 - Sum_{d|n, n/d == 3 (mod 4)} d^2.

Original entry on oeis.org

1, 4, 8, 16, 26, 32, 48, 64, 73, 104, 120, 128, 170, 192, 208, 256, 290, 292, 360, 416, 384, 480, 528, 512, 651, 680, 656, 768, 842, 832, 960, 1024, 960, 1160, 1248, 1168, 1370, 1440, 1360, 1664, 1682, 1536, 1848, 1920, 1898, 2112, 2208, 2048, 2353, 2604
Offset: 1

Views

Author

N. J. A. Sloane, Dec 23 1999

Keywords

Comments

Number 7 of the 74 eta-quotients listed in Table I of Martin (1996).
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Multiplicative because it is the Dirichlet convolution of A000290 = n^2 and A101455 = [1 0 -1 0 1 0 -1 ...], which are both multiplicative. - Christian G. Bower, May 17 2005

Examples

			G.f. = q + 4*q^2 + 8*q^3 + 16*q^4 + 26*q^5 + 32*q^6 + 48*q^7 + 64*q^8 + ...
		

Crossrefs

Programs

  • Haskell
    a050470 n = a050461 n - a050465 n  -- Reinhard Zumkeller, Mar 06 2012
    
  • Magma
    Basis( ModularForms( Gamma1(4), 3), 51) [2]; /* Michael Somos, May 17 2015 */
    
  • Mathematica
    a[ n_] := SeriesCoefficient[ q (QPochhammer[ q^2]^3 (QPochhammer[ q^4] / QPochhammer[ q])^2)^2, {q, 0, n}]; (* Michael Somos, May 17 2015 *)
    a[ n_] := SeriesCoefficient[ (EllipticTheta[ 3, 0, q] EllipticTheta[ 2, 0, q]^2 / 4)^2, {q, 0, n}]; (* Michael Somos, May 17 2015 *)
    a[ n_] := If[ n < 1, 0, Sum[ d^2 Mod[n/d, 2] (-1)^Quotient[n/d, 2], {d, Divisors@n}]]; (* Michael Somos, May 17 2015 *)
    s[n_] := If[OddQ[n], (-1)^((n-1)/2), 0]; (* A101455 *)
    f[p_, e_] := (p^(2*e+2) - s[p]^(e+1))/(p^2 - s[p]); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Nov 04 2023 *)
  • PARI
    {a(n) = if( n<1, 0, sumdiv( n, d, d^2 * (n/d%2) * (-1)^(n/d\2)))};
    
  • PARI
    {a(n) = my(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( eta(x^2 + A)^6 * (eta(x^4 + A) / eta(x + A))^4, n))}; /* Michael Somos, May 17 2015 */
    
  • Python
    from math import prod
    from sympy import factorint
    def A050470(n): return prod((p**(e+1<<1)-(m:=(0,1,0,-1)[p&3]))//(p**2-m) for p, e in factorint(n).items()) # Chai Wah Wu, Jun 21 2024

Formula

G.f.: Sum_{n>=1} n^2*x^n/(1+x^(2*n)). - Vladeta Jovovic, Oct 16 2002
From Michael Somos, Aug 08 2005: (Start)
Euler transform of period 4 sequence [ 4, -2, 4, -6, ...].
Expansion of eta(q^2)^6 * eta(q^4)^4 / eta(q)^4 in powers of q.
G.f.: x Product_{k>0} (1 + x^k)^4 * (1 - x^(2*k))^2 * (1 - x^(4*k))^4.
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^4)) where f(u, v, w) = u*w * (u - 8*v) * (v - 4*w) - v^2 * (v - 8*w)^2. (End)
G.f.: Sum_{k>0} Kronecker(-4, k) * x^k * (1 + x^k) / (1 - x^k)^3. - Michael Somos, Sep 02 2005
Expansion of q * phi(q)^2 * psi(q^2)^4 in powers of q where phi(), psi() are Ramanujan theta functions. - Michael Somos, Aug 15 2007
G.f. is a period 1 Fourier series which satisfies f(-1 / (4 t)) = (1/2) (t/i)^3 g(t) where q = exp(2 Pi i t) and g() is the g.f. for A120030.
a(n) = A050461(n) - A050465(n). - Reinhard Zumkeller, Mar 06 2012
Multiplicative with a(p^e) = ((p^2)^(e+1) - Chi(p)^(e+1))/(p^2 - Chi(p)), Chi = A101455. - Jianing Song, Oct 30 2019
Sum_{k=1..n} a(k) ~ c * n^3 / 3, where c = Pi^3/32 (A153071). - Amiram Eldar, Nov 04 2023
a(n) = Sum_{d|n} (n/d)^2*sin(d*Pi/2). - Ridouane Oudra, Sep 26 2024
Showing 1-10 of 62 results. Next