A006752 -id:A006752 - cp's OEIS frontend

A016754 Odd squares: a(n) = (2n+1)^2. Also centered octagonal numbers.

1, 9, 25, 49, 81, 121, 169, 225, 289, 361, 441, 529, 625, 729, 841, 961, 1089, 1225, 1369, 1521, 1681, 1849, 2025, 2209, 2401, 2601, 2809, 3025, 3249, 3481, 3721, 3969, 4225, 4489, 4761, 5041, 5329, 5625, 5929, 6241, 6561, 6889, 7225, 7569, 7921, 8281, 8649, 9025

Offset: 0

N. J. A. Sloane

The brown rat (rattus norwegicus) breeds very quickly. It can give birth to other rats 7 times a year, starting at the age of three months. The average number of pups is 8. The present sequence gives the total number of rats, when the intervals are 12/7 of a year and a young rat starts having offspring at 24/7 of a year. - Hans Isdahl, Jan 26 2008
Numbers n such that tau(n) is odd where tau(x) denotes the Ramanujan tau function (A000594). - Benoit Cloitre, May 01 2003
If Y is a fixed 2-subset of a (2n+1)-set X then a(n-1) is the number of 3-subsets of X intersecting Y. - Milan Janjic, Oct 21 2007
Binomial transform of [1, 8, 8, 0, 0, 0, ...]; Narayana transform (A001263) of [1, 8, 0, 0, 0, ...]. - Gary W. Adamson, Dec 29 2007
All terms of this sequence are of the form 8k+1. For numbers 8k+1 which aren't squares see A138393. Numbers 8k+1 are squares iff k is a triangular number from A000217. And squares have form 4n(n+1)+1. - Artur Jasinski, Mar 27 2008
Sequence arises from reading the line from 1, in the direction 1, 25, ... and the line from 9, in the direction 9, 49, ..., in the square spiral whose vertices are the squares A000290. - Omar E. Pol, May 24 2008
Equals the triangular numbers convolved with [1, 6, 1, 0, 0, 0, ...]. - Gary W. Adamson & Alexander R. Povolotsky, May 29 2009
First differences: A008590(n) = a(n) - a(n-1) for n>0. - Reinhard Zumkeller, Nov 08 2009
Central terms of the triangle in A176271; cf. A000466, A053755. - Reinhard Zumkeller, Apr 13 2010
Odd numbers with odd abundance. Odd numbers with even abundance are in A088828. Even numbers with odd abundance are in A088827. Even numbers with even abundance are in A088829. - Jaroslav Krizek, May 07 2011
Appear as numerators in the non-simple continued fraction expansion of Pi-3: Pi-3 = K_{k>=1} (1-2*k)^2/6 = 1/(6+9/(6+25/(6+49/(6+...)))), see also the comment in A007509. - Alexander R. Povolotsky, Oct 12 2011
Ulam's spiral (SE spoke). - Robert G. Wilson v, Oct 31 2011
All terms end in 1, 5 or 9. Modulo 100, all terms are among { 1, 9, 21, 25, 29, 41, 49, 61, 69, 81, 89 }. - M. F. Hasler, Mar 19 2012
Right edge of both triangles A214604 and A214661: a(n) = A214604(n+1,n+1) = A214661(n+1,n+1). - Reinhard Zumkeller, Jul 25 2012
Also: Odd numbers which have an odd sum of divisors (= sigma = A000203). - M. F. Hasler, Feb 23 2013
Consider primitive Pythagorean triangles (a^2 + b^2 = c^2, gcd(a, b) = 1) with hypotenuse c (A020882) and respective even leg b (A231100); sequence gives values c-b, sorted with duplicates removed. - K. G. Stier, Nov 04 2013
For n>1 a(n) is twice the area of the irregular quadrilateral created by the points ((n-2)*(n-1),(n-1)*n/2), ((n-1)*n/2,n*(n+1)/2), ((n+1)*(n+2)/2,n*(n+1)/2), and ((n+2)*(n+3)/2,(n+1)*(n+2)/2). - J. M. Bergot, May 27 2014
Number of pairs (x, y) of Z^2, such that max(abs(x), abs(y)) <= n. - Michel Marcus, Nov 28 2014
Except for a(1)=4, the number of active (ON, black) cells in n-th stage of growth of two-dimensional cellular automaton defined by "Rule 737", based on the 5-celled von Neumann neighborhood. - Robert Price, May 23 2016
a(n) is the sum of 2n+1 consecutive numbers, the first of which is n+1. - Ivan N. Ianakiev, Dec 21 2016
a(n) is the number of 2 X 2 matrices with all elements in {0..n} with determinant = 2*permanent. - Indranil Ghosh, Dec 25 2016
Engel expansion of Pi*StruveL_0(1)/2 where StruveL_0(1) is A197037. - Benedict W. J. Irwin, Jun 21 2018
Consider all Pythagorean triples (X,Y,Z=Y+1) ordered by increasing Z; the segments on the hypotenuse {p = a(n)/A001844(n), q = A060300(n)/A001844(n) = A001844(n) - p} and their ratio p/q = a(n)/A060300(n) are irreducible fractions in Q\Z. X values are A005408, Y values are A046092, Z values are A001844. - Ralf Steiner, Feb 25 2020
a(n) is the number of large or small squares that are used to tile primitive squares of type 2 (A344332). - Bernard Schott, Jun 03 2021
Also, positive odd integers with an odd number of odd divisors (for similar sequence with 'even', see A348005). - Bernard Schott, Nov 21 2021
a(n) is the least odd number k = x + y, with 0 < x < y, such that there are n distinct pairs (x,y) for which x*y/k is an integer; for example, a(2) = 25 and the two corresponding pairs are (5,20) and (10,15). The similar sequence with 'even' is A016742 (see Comment of Jan 26 2018). - Bernard Schott, Feb 24 2023
From Peter Bala, Jan 03 2024: (Start)
The sequence terms are the exponents of q in the series expansions of the following infinite products:
1) q*Product_{n >= 1} (1 - q^(16*n))*(1 + q^(8*n)) = q + q^9 + q^25 + q^49 + q^81 + q^121 + q^169 + ....
2) q*Product_{n >= 1} (1 + q^(16*n))*(1 - q^(8*n)) = q - q^9 - q^25 + q^49 + q^81 - q^121 - q^169 + + - - ....
3) q*Product_{n >= 1} (1 - q^(8*n))^3 = q - 3*q^9 + 5*q^25 - 7*q^49 + 9*q^81 - 11*q^121 + 13*q^169 - + ....
4) q*Product_{n >= 1} ( (1 + q^(8*n))*(1 - q^(16*n))/(1 + q^(16*n)) )^3 = q + 3*q^9 - 5*q^25 - 7*q^49 + 9*q^81 + 11*q^121 - 13*q^169 - 15*q^225 + + - - .... (End)

L. Lorentzen and H. Waadeland, Continued Fractions with Applications, North-Holland 1992, p. 586.

Paolo Xausa, Table of n, a(n) for n = 0..10000 (terms 0..1000 from T. D. Noe).
Jeremiah Bartz, Bruce Dearden, and Joel Iiams, Classes of Gap Balancing Numbers, arXiv:1810.07895 [math.NT], 2018.
Bruce C. Berndt and Ken Ono, Ramanujan's unpublished manuscript on the partition and tau functions with proofs and commentary, Séminaire Lotharingien de Combinatoire, B42c (1999), 63 pp.
John Elias, Illustration: 8-fold Square Progression of Ulam's Spiral.
Milan Janjic, Two Enumerative Functions.
Scientific American, Cover of the March 1964 issue.
Amelia Carolina Sparavigna, Groupoids of OEIS A002378 and A016754 Numbers (oblong and odd square numbers), Politecnico di Torino (Italy, 2019).
Leo Tavares, Illustration: Diamond Triangles.
Leo Tavares, Illustration: Diamond Stars.
Eric Weisstein's World of Mathematics, Moore Neighborhood.
R. Yin, J. Mu, and T. Komatsu, The p-Frobenius Number for the Triple of the Generalized Star Numbers, Preprints 2024, 2024072280. See p. 2.
Index entries for sequences related to centered polygonal numbers.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A000290, A000384, A001263, A001539, A001844, A003881, A005408, A006752, A014105, A016742, A016802, A016814, A016826, A016838, A033996, A046092, A060300, A138393, A167661, A167700.
Cf. A000447 (partial sums).
Cf. A005917, A344330, A344332.
Cf. A348005, A379481 [= a(A048673(n)-1)].
Partial sums of A022144.
Positions of odd terms in A341528.
Sequences on the four axes of the square spiral: Starting at 0: A001107, A033991, A007742, A033954; starting at 1: A054552, A054556, A054567, A033951.
Sequences on the four diagonals of the square spiral: Starting at 0: A002939 = 2*A000384, A016742 = 4*A000290, A002943 = 2*A014105, A033996 = 8*A000217; starting at 1: A054554, A053755, A054569, A016754.
Sequences obtained by reading alternate terms on the X and Y axes and the two main diagonals of the square spiral: Starting at 0: A035608, A156859, A002378 = 2*A000217, A137932 = 4*A002620; starting at 1: A317186, A267682, A002061, A080335.
Cf. A014634, A003154.

a016754 n = a016754_list !! n
a016754_list = scanl (+) 1 $ tail a008590_list
-- Reinhard Zumkeller, Apr 02 2012

[n^2: n in [1..100 by 2]]; // Vincenzo Librandi, Jan 03 2017

Range[1, 100, 2]^2 (* Paolo Xausa, Feb 13 2025 *)

A016754(n):=(n+n+1)^2$
makelist(A016754(n),n,0,20); /* Martin Ettl, Nov 12 2012 */

A016754(n)=(n<<1+1)^2 \\ Charles R Greathouse IV, Jun 16 2011, corrected and edited by M. F. Hasler, Apr 11 2023

def A016754(n): return ((n<<1)|1)**2 # Chai Wah Wu, Jul 06 2023

a(n) = 1 + Sum_{i=1..n} 8*i = 1 + 8*A000217(n). - Xavier Acloque, Jan 21 2003; Zak Seidov, May 07 2006; Robert G. Wilson v, Dec 29 2010
O.g.f.: (1+6*x+x^2)/(1-x)^3. - R. J. Mathar, Jan 11 2008
a(n) = 4*n*(n + 1) + 1 = 4*n^2 + 4*n + 1. - Artur Jasinski, Mar 27 2008
a(n) = A061038(2+4n). - Paul Curtz, Oct 26 2008
Sum_{n>=0} 1/a(n) = Pi^2/8 = A111003. - Jaume Oliver Lafont, Mar 07 2009
a(n) = A000290(A005408(n)). - Reinhard Zumkeller, Nov 08 2009
a(n) = a(n-1) + 8*n with n>0, a(0)=1. - Vincenzo Librandi, Aug 01 2010
a(n) = A033951(n) + n. - Reinhard Zumkeller, May 17 2009
a(n) = A033996(n) + 1. - Omar E. Pol, Oct 03 2011
a(n) = (A005408(n))^2. - Zak Seidov, Nov 29 2011
From George F. Johnson, Sep 05 2012: (Start)
a(n+1) = a(n) + 4 + 4*sqrt(a(n)).
a(n-1) = a(n) + 4 - 4*sqrt(a(n)).
a(n+1) = 2*a(n) - a(n-1) + 8.
a(n+1) = 3*a(n) - 3*a(n-1) + a(n-2).
(a(n+1) - a(n-1))/8 = sqrt(a(n)).
a(n+1)*a(n-1) = (a(n)-4)^2.
a(n) = 2*A046092(n) + 1 = 2*A001844(n) - 1 = A046092(n) + A001844(n).
Limit_{n -> oo} a(n)/a(n-1) = 1. (End)
a(n) = binomial(2*n+2,2) + binomial(2*n+1,2). - John Molokach, Jul 12 2013
E.g.f.: (1 + 8*x + 4*x^2)*exp(x). - Ilya Gutkovskiy, May 23 2016
a(n) = A101321(8,n). - R. J. Mathar, Jul 28 2016
Product_{n>=1} A033996(n)/a(n) = Pi/4. - Daniel Suteu, Dec 25 2016
a(n) = A014105(n) + A000384(n+1). - Bruce J. Nicholson, Nov 11 2017
a(n) = A003215(n) + A002378(n). - Klaus Purath, Jun 09 2020
From Amiram Eldar, Jun 20 2020: (Start)
Sum_{n>=0} a(n)/n! = 13*e.
Sum_{n>=0} (-1)^(n+1)*a(n)/n! = 3/e. (End)
Sum_{n>=0} (-1)^n/a(n) = A006752. - Amiram Eldar, Oct 10 2020
From Amiram Eldar, Jan 28 2021: (Start)
Product_{n>=0} (1 + 1/a(n)) = cosh(Pi/2).
Product_{n>=1} (1 - 1/a(n)) = Pi/4 (A003881). (End)
From Leo Tavares, Nov 24 2021: (Start)
a(n) = A014634(n) - A002943(n). See Diamond Triangles illustration.
a(n) = A003154(n+1) - A046092(n). See Diamond Stars illustration. (End)
From Peter Bala, Mar 11 2024: (Start)
Sum_{k = 1..n+1} 1/(k*a(k)*a(k-1)) = 1/(9 - 3/(17 - 60/(33 - 315/(57 - ... - n^2*(4*n^2 - 1)/((2*n + 1)^2 + 2*2^2 ))))).
3/2 - 2*log(2) = Sum_{k >= 1} 1/(k*a(k)*a(k-1)) = 1/(9 - 3/(17 - 60/(33 - 315/(57 - ... - n^2*(4*n^2 - 1)/((2*n + 1)^2 + 2*2^2 - ... ))))).
Row 2 of A142992. (End)
From Peter Bala, Mar 26 2024: (Start)
8*a(n) = (2*n + 1)*(a(n+1) - a(n-1)).
Sum_{n >= 0} (-1)^n/(a(n)*a(n+1)) = 1/2 - Pi/8 = 1/(9 + (1*3)/(8 + (3*5)/(8 + ... + (4*n^2 - 1)/(8 + ... )))). For the continued fraction use Lorentzen and Waadeland, p. 586, equation 4.7.9 with n = 1. Cf. A057813. (End)

Additional description from Terrel Trotter, Jr., Apr 06 2002

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

N. J. A. Sloane, Simon Plouffe

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

			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

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.

Reinhard Zumkeller, Table of n, a(n) for n = 0..1000
Sanjar M. Abrarov and Brendan M. Quine, A Viète-like formula for pi based on infinite sum of the arctangent functions with nested radicals, figshare, 4509014, (2017).
Peter Bala, Arctanh(z) and the Legendre polynomials
Jonathan M. Borwein, Peter B. Borwein, and Karl Dilcher, Pi, Euler numbers and asymptotic expansions, Amer. Math. Monthly, 96 (1989), 681-687.
Antonio Gracia Llorente, Shifting Constants Through Infinite Product Transformations, OSF Preprint, 2024.
Ronald K. Hoeflin, Titan Test.
Richard J. Mathar, Table of Dirichlet L-Series and Prime Zeta Modulo Functions for Small Moduli, arXiv:1008.2547 [math.NT], 2010-2015.
Literate Programs, Pi with Machin's formula (Haskell).
Michael Penn, A surprising appearance of pie!, YouTube video, 2020.
Michael Penn, Transforming normal identities into "crazy" ones, YouTube video, 2022.
Srinivasa Ramanujan, Question 353, J. Ind. Math. Soc.
Eric Weisstein's World of Mathematics, Prime Products.
Wikipedia, Leibniz formula for Pi.
Index entries for sequences from "Goedel, Escher, Bach".
Index entries for transcendental numbers.

Cf. A000796, A001586, A071904, A019669, A197723, A347152.
Cf. A000182, A000364, A024235, A278080, A278195.
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.

-- 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

R:= RealField(100); Pi(R)/4; // G. C. Greubel, Mar 08 2018

Maple
```
evalf(Pi/4) ;
```

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 *)

Pi/4 \\ Charles R Greathouse IV, Jul 07 2014

# 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

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

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

A101455 a(n) = 0 for even n, a(n) = (-1)^((n-1)/2) for odd n. Periodic sequence 1,0,-1,0,...

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

Offset: 0

Gerald McGarvey, Jan 20 2005

Called X(n) (i.e., Chi(n)) in Hardy and Wright (p. 241), who show that X(n*m) = X(n)*X(m) for all n and m (i.e., X(n) is completely multiplicative) since (n*m - 1)/2 - (n - 1)/2 - (m - 1)/2 = (n - 1)*(m - 1)/2 == 0 (mod 2) when n and m are odd.
Same as A056594 but with offset 1.
From R. J. Mathar, Jul 15 2010: (Start)
The sequence is the non-principal Dirichlet character mod 4. (The principal character is A000035.)
Associated Dirichlet L-functions are for example L(1,chi) = Sum_{n>=1} a(n)/n = A003881, or L(2,chi) = Sum_{n>=1} a(n)/n^2 = A006752, or L(3,chi) = Sum_{n>=1} a(n)/n^3 = A153071. (End)
a(n) is a strong elliptic divisibility sequence t_n as given in [Kimberling, p. 16] where x = 0, y = -1, z is arbitrary. - Michael Somos, Nov 27 2019

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

T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1986, page 139, k=4, Chi_2(n).
G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 5th ed., Oxford Univ. Press, 1979, p. 241.

Muniru A Asiru, Table of n, a(n) for n = 0..1000 (a(0) prepended by Jianing Song)
Étienne Fouvry, Claude Levesque, Michel Waldschmidt, Representation of integers by cyclotomic binary forms, arXiv:1712.09019 [math.NT], 2017.
Clark Kimberling, Strong divisibility sequences and some conjectures, Fib. Quart., 17 (1979), 13-17.
Grant Sanderson, Pi hiding in prime regularities, 3Blue1Brown video (2017).
Index to divisibility sequences
Index entries for linear recurrences with constant coefficients, signature (0,-1).

Cf. A002654, A009116, A009454, A056594, A108520, A111335, A126120, A146559.

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), this sequence (d=-4), A102283 (d=-3), A080891 (d=5), A091337 (d=8), A110161 (d=12), A011583 (d=13), A011584 (d=17), A322829 (d=21), A322796 (d=24).

a := [1, 0];; for n in [3..10^2] do a[n] := a[n-2]; od; a; # Muniru A Asiru, Feb 02 2018

m:=75; R:=PowerSeriesRing(Integers(), m); Coefficients(R!(x/(1+x^2))); // G. C. Greubel, Aug 23 2018

a := n -> `if`(n mod 2=0, 0, (-1)^((n-1)/2)):
seq(a(n), n=1..10^3); # Muniru A Asiru, Feb 02 2018

a[ n_] := {1, 0, -1, 0}[[ Mod[ n, 4, 1]]]; (* Michael Somos, Jan 13 2014 *)
LinearRecurrence[{0, -1}, {1, 0}, 75] (* G. C. Greubel, Aug 23 2018 *)

{a(n) = if( n%2, (-1)^(n\2))}; /* Michael Somos, Sep 02 2005 */

{a(n) = kronecker( -4, n)}; /* Michael Somos, Mar 30 2012 */

def A101455(n): return (0,1,0,-1)[n&3] # Chai Wah Wu, Jun 21 2024

Multiplicative with a(2^e) = 0, a(p^e) = (-1)^((p^e-1)/2) otherwise. - Mitch Harris May 17 2005
Euler transform of length 4 sequence [0, -1, 0, 1]. - Michael Somos, Sep 02 2005
G.f.: (x - x^3)/(1 - x^4) = x/(1 + x^2). - Michael Somos, Sep 02 2005
G.f. A(x) satisfies: 0 = f(A(x), A(x^2)) where f(u, v) = v - u^2 * (1 + 2*v). - Michael Somos, Aug 04 2011
a(n + 4) = a(n), a(n + 2) = a(-n) = -a(n), a(2*n) = 0, a(2*n + 1) = (-1)^n for all n in Z. - Michael Somos, Aug 04 2011
a(n + 1) = A056594(n). - Michael Somos, Jan 13 2014
REVERT transform is A126120. STIRLING transform of A009454. BINOMIAL transform is A146559. BINOMIAL transform of A009116. BIN1 transform is A108520. MOBIUS transform of A002654. EULER transform is A111335. - Michael Somos, Mar 30 2012
Completely multiplicative with a(p) = 2 - (p mod 4). - Werner Schulte, Feb 01 2018
a(n) = (-(n mod 2))^binomial(n, 2). - Peter Luschny, Sep 08 2018
a(n) = sin(n*Pi/2) = Im(i^n) where i is the imaginary unit. - Jianing Song, Sep 09 2018
From Jianing Song, Nov 14 2018: (Start)
a(n) = ((-4)/n) (or more generally, ((-4^i)/n) for i > 0), where (k/n) is the Kronecker symbol.
E.g.f.: sin(x).
Dirichlet g.f. is the Dirichlet beta function.
a(n) = A091337(n)*A188510(n). (End)

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

A004003 Number of domino tilings (or dimer coverings) of a 2n X 2n square.

Original entry on oeis.org

1, 2, 36, 6728, 12988816, 258584046368, 53060477521960000, 112202208776036178000000, 2444888770250892795802079170816, 548943583215388338077567813208427340288, 1269984011256235834242602753102293934298576249856

Offset: 0

N. J. A. Sloane, Greg Huber

A099390 is the main entry for domino tilings (or dimer tilings) of a rectangle.
The numbers of domino tilings in A006253, A004003, A006125 give the number of perfect matchings in the relevant graphs. There are results of Jockusch and Ciucu that if a planar graph has a rotational symmetry then the number of perfect matchings is a square or twice a square - this applies to these 3 sequences. - Dan Fux (dan.fux(AT)OpenGaia.com or danfux(AT)OpenGaia.com), Apr 12 2001
Christine Bessenrodt points out that Pachter (1997) shows that a(n) is divisible by 2^n (cf. A065072).
a(n) is the number of different ways to cover a 2n X 2n lattice with 2n^2 dominoes. John and Sachs show that a(n) = 2^n*B(n)^2, where B(n) == n+1 (mod 32) when n is even and B(n) == (-1)^((n-1)/2)*n (mod 32) when n is odd. - Yong Kong (ykong(AT)curagen.com), May 07 2000

			The 36 solutions for the 4 X 4 board, from Neven Juric, May 14 2008:
A01 = {(1,2), (3,4), (5,6), (7,8), (9,10), (11,12), (13,14), (15,16)}
A02 = {(1,2), (3,4), (5,6), (7,11), (9,10), (8,12), (13,14), (15,16)}
A03 = {(1,2), (3,4), (5,9), (6,7), (10,11), (8,12), (13,14), (15,16)}
A04 = {(1,2), (3,4), (5,9), (6,10), (7,8), (11,12), (13,14), (15,16)}
A05 = {(1,2), (3,4), (5,9), (6,10), (7,11), (8,12), (13,14), (15,16)}
A06 = {(1,2), (3,4), (5,6), (7,8), (9,10), (13,14), (11,15), (12,16)}
A07 = {(1,2), (3,4), (5,9), (6,10), (7,8), (11,15), (13,14), (12,16)}
A08 = {(1,2), (3,4), (5,6), (7,8), (9,13), (10,14), (11,12), (15,16)}
A09 = {(1,2), (3,4), (5,6), (7,11), (8,12), (9,13), (10,14), (15,16)}
A10 = {(1,2), (3,4), (5,6), (7,8), (9,13), (10,11), (14,15), (12,16)}
A11 = {(1,2), (3,4), (5,6), (7,8), (9,13), (10,14), (11,15), (12,16)}
A12 = {(1,2), (5,6), (3,7), (4,8), (9,10), (11,12), (13,14), (15,16)}
A13 = {(1,2), (3,7), (4,8), (5,9), (6,10), (11,12), (13,14), (15,16)}
A14 = {(1,2), (5,6), (3,7), (4,8), (9,10), (13,14), (11,15), (12,16)}
A15 = {(1,2), (3,7), (4,8), (6,10), (5,9), (11,15), (12,16), (13,14)}
A16 = {(1,2), (3,7), (4,8), (5,6), (9,13), (10,14), (11,12), (15,16)}
A17 = {(1,2), (3,7), (4,8), (5,6), (9,13), (10,11), (14,15), (12,16)}
A18 = {(1,2), (5,6), (3,7), (4,8), (9,13), (10,14), (11,15), (12,16)}
A19 = {(1,5), (2,6), (3,4), (7,8), (9,10), (11,12), (13,14), (15,16)}
A20 = {(1,5), (2,6), (3,4), (7,11), (8,12), (9,10), (13,14), (15,16)}
A21 = {(1,5), (3,4), (2,6), (9,10), (7,8), (11,15), (13,14), (12,16)}
A22 = {(1,5), (2,6), (3,4), (7,8), (9,13), (10,14), (11,12), (15,16)}
A23 = {(1,5), (2,6), (3,4), (7,11), (8,12), (9,13), (10,14), (15,16)}
A24 = {(1,5), (2,6), (3,4), (7,8), (9,13), (10,11), (14,15), (12,16)}
A25 = {(1,5), (2,6), (3,4), (7,8), (9,13), (10,14), (11,15), (12,16)}
A26 = {(1,5), (2,3), (6,7), (4,8), (9,10), (11,12), (13,14), (15,16)}
A27 = {(1,5), (2,6), (3,7), (4,8), (9,10), (11,12), (13,14), (15,16)}
A28 = {(1,5), (2,3), (6,7), (4,8), (9,10), (11,15), (13,14), (12,16)}
A29 = {(1,5), (2,6), (3,7), (4,8), (9,10), (13,14), (11,15), (12,16)}
A30 = {(1,5), (2,3), (6,7), (4,8), (9,13), (10,14), (11,12), (15,16)}
A31 = {(1,5), (2,6), (3,7), (4,8), (9,13), (10,14), (11,12), (15,16)}
A32 = {(1,5), (2,6), (3,7), (4,8), (9,13), (10,14), (11,15), (12,16)}
A33 = {(1,5), (2,6), (3,7), (4,8), (9,13), (10,11), (14,15), (12,16)}
A34 = {(1,5), (2,3), (4,8), (6,10), (7,11), (9,13), (14,15), (12,16)}
A35 = {(1,5), (2,3), (6,7), (4,8), (9,13), (10,14), (11,15), (12,16)}
A36 = {(1,5), (2,3), (6,7), (4,8), (9,13), (10,11), (14,15), (12,16)}

Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 569.
S. R. Finch, Mathematical Constants, Cambridge, 2003, pp. 406-412.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
Darko Veljan, Kombinatorika s teorijom grafova (Croatian) (Combinatorics with Graph Theory) mentions the value 12988816 = 2^4*901^2 for the 8 X 8 case on page 4.

Alois P. Heinz, Table of n, a(n) for n = 0..44 (first 26 terms from T. D. Noe)
M. Aanjaneya and S. P. Pal, Faultfree tromino tilings of rectangles, arXiv:math/0610925 [math.CO], 2006.
N. Allegra, Exact solution of the 2d dimer model: Corner free energy, correlation functions and combinatorics, arXiv:1410.4131 [cond-mat.stat-mech], 2014.
M. Ciucu, Enumeration of perfect matchings in graphs with reflective symmetry, Journal of Combinatorial Theory, Series A, Volume 77, Issue 1, January 1997, Pages 67-97.
H. Cohn, 2-adic behavior of numbers of domino tilings, arXiv:math/0008222 [math.CO], 2000.
Steven R. Finch, The Dimer Problem [From Steven Finch, Apr 20 2019]
M. E. Fisher, Statistical mechanics of dimers on a plane lattice, Physical Review, 124 (1961), 1664-1672.
P. Flajolet and R. Sedgewick, Analytic Combinatorics, 2009; see page 363
Laura Florescu, Daniela Morar, David Perkinson, Nicholas Salter, and Tianyuan Xu, Sandpiles and Dominos, Electronic Journal of Combinatorics, Volume 22(1), 2015.
Laura Florescu, Daniela Morar, David Perkinson, Nicholas Salter, and Tianyuan Xu, Illustration for a(2) = 36 [Fig. 9 from "Sandpiles and Dominos"]
W. Jockusch, Perfect matchings and perfect squares J. Combin. Theory Ser. A 67 (1994), no. 1, 100-115.
Peter E. John and Horst Sachs, On a strange observation in the theory of the dimer problem, Discrete Math. 216 (2000), no. 1-3, 211-219.
Adrien Kassel, Le modèle de dimères, Images des Mathématiques, CNRS, 2016. [In French]
P. W. Kasteleyn, The Statistics of Dimers on a Lattice, Physica, 27 (1961), 1209-1225.
P. W. Kasteleyn, Dimer statistics and phase transitions, J. Mathematical Phys. 4 1963 287-293. MR0153427 (27 #3394).
Viet-Ha Nguyen, Kévin Perrot, and Mathieu Vallet, NP-completeness of the game Kingdomino(TM), Theoretical Computer Science (2020) Vol. 822, 23-35. See also arXiv:1909.02849, [cs.CC], 2019.
Lior Pachter, Combinatorial approaches and conjectures ..., Elec. J. Comb. 4 (1997) #R29.
James Propp, Some 2-Adic Conjectures Concerning Polyomino Tilings of Aztec Diamonds, Integers (2023) Vol. 23, Art. A30.
Jaime Rangel-Mondragon, Polyominoes and Related Families, The Mathematica Journal, 9:3 (2005), 609-640.
N. J. A. Sloane, Illustration for a(2) = 36 [Slide from an old talk I gave]
R. P. Stanley, A combinatorial miscellany
Eric Weisstein's World of Mathematics, Domino Tiling
Eric Weisstein, Illustration for a(2) = 36, from Domino Tilings web page (see previous link) [Included with permission]
Index entries for sequences related to dominoes

Cf. A028420, A006253, A006125, A065072, A124997, A239273, A256043.

Main diagonal of array A099390 or A187596.

f := n->2^(2*n^2)*product(product(cos(i*Pi/(2*n+1))^2+cos(j*Pi/(2*n+1))^2,j=1..n),i=1..n); for k from 0 to 12 do round(evalf(f(k),300)) od;

a[n_] := Round[ N[ Product[ 2*Cos[(2i*Pi)/(2n+1)] + 2*Cos[(2j*Pi)/(2n+1)] + 4,  {i, 1, n}, {j, 1, n}], 300] ]; Table[a[n], {n, 0, 12}] (* Jean-François Alcover, Jan 04 2012, after Maple *)
Table[Sqrt[Resultant[ChebyshevU[2*n, x/2], ChebyshevU[2*n, I*x/2], x]], {n, 0, 12}] (* Vaclav Kotesovec, Dec 30 2020 *)

{a(n) = sqrtint(polresultant(polchebyshev(2*n, 2, x/2), polchebyshev(2*n, 2, I*x/2)))} \\ Seiichi Manyama, Apr 13 2020

from math import isqrt
from sympy.abc import x
from sympy import I, resultant, chebyshevu
def A004003(n): return isqrt(resultant(chebyshevu(n<<1,x/2),chebyshevu(n<<1,I*x/2))) if n else 1 # Chai Wah Wu, Nov 07 2023

a(n) = A099390(2n,2n).
a(n) = Product_{j=1..n} Product_{k=1..n} (4*cos(j*Pi/(2*n+1))^2 + 4*cos(k*Pi/(2*n+1))^2). - N. J. A. Sloane, Mar 16 2015
a(n) = 2^n * A065072(n)^2. - Alois P. Heinz, Nov 22 2018
a(n)^2 = Resultant(U(2*n,x/2), U(2*n,i*x/2)), where U(n,x) is a Chebyshev polynomial of the second kind and i = sqrt(-1). - Seiichi Manyama, Apr 13 2020
a(n) ~ 2 * (sqrt(2)-1)^(2*n+1) * exp(G*(2*n+1)^2/Pi), where G is Catalan's constant A006752. - Vaclav Kotesovec, Dec 30 2020

Corrected and extended by David Radcliffe

A003983 Array read by antidiagonals with T(n,k) = min(n,k).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 2, 3, 2, 1, 1, 2, 3, 3, 2, 1, 1, 2, 3, 4, 3, 2, 1, 1, 2, 3, 4, 4, 3, 2, 1, 1, 2, 3, 4, 5, 4, 3, 2, 1, 1, 2, 3, 4, 5, 5, 4, 3, 2, 1, 1, 2, 3, 4, 5, 6, 5, 4, 3, 2, 1, 1, 2, 3, 4, 5, 6, 6, 5, 4, 3, 2, 1, 1, 2, 3, 4, 5, 6, 7, 6, 5, 4, 3, 2, 1, 1, 2, 3, 4, 5, 6, 7, 7, 6, 5, 4, 3, 2, 1

Offset: 1

Marc LeBrun

Also, "correlation triangle" for the constant sequence 1. - Paul Barry, Jan 16 2006
Antidiagonal sums are in A002620.
As a triangle, row sums are A002620. T(2n,n)=n+1. Diagonal sums are A001399. Construction: Take antidiagonal triangle of MM^T where M is the sequence array for the constant sequence 1 (lower triangular matrix with all 1's). - Paul Barry, Jan 16 2006
From Franklin T. Adams-Watters, Sep 25 2011: (Start)
As a triangle, count up to ceiling(n/2) and back down again (repeating the central term when n is even).
When the first two instances of each number are removed from the sequence, the original sequence is recovered.
(End)

			Triangle version begins
  1;
  1, 1;
  1, 2, 1;
  1, 2, 2, 1;
  1, 2, 3, 2, 1;
  1, 2, 3, 3, 2, 1;
  1, 2, 3, 4, 3, 2, 1;
  1, 2, 3, 4, 4, 3, 2, 1;
  1, 2, 3, 4, 5, 4, 3, 2, 1;
  ...

Reinhard Zumkeller, Rows n=1..100 of triangle, flattened

Cf. A002620, A001399, A087062, A115236, A115237, A124258, A006752, A120268, A173945, A173947, A173948, A173949, A173953, A173954, A173955, A173973, A173982-A173986, A004197.

a003983 n k = a003983_tabl !! (n-1) !! (k-1)
a003983_tabl = map a003983_row [1..]
a003983_row n = hs ++ drop m (reverse hs)
   where hs = [1..n' + m]
         (n',m) = divMod n 2
-- Reinhard Zumkeller, Aug 14 2011

a(n) = min(floor(1/2 + sqrt(2*n)) - (2*n + round(sqrt(2*n)) - round(sqrt(2*n))^2)/2+1, (2*n + round(sqrt(2*n)) - round(sqrt(2*n))^2)/2) # Leonid Bedratyuk, Dec 13 2009

Flatten[Table[Min[n-k+1, k], {n, 1, 14}, {k, 1, n}]] (* Jean-François Alcover, Feb 23 2012 *)

T(n,k) = min(n,k) \\ Charles R Greathouse IV, Feb 06 2017

from math import isqrt
def A003983(n):
    a = (m:=isqrt(k:=n<<1))+(k>m*(m+1))
    x = n-(a*(a-1)>>1)
    return min(x,a-x+1) # Chai Wah Wu, Jun 14 2025

Number triangle T(n, k) = Sum_{j=0..n} [j<=k][j<=n-k]. - Paul Barry, Jan 16 2006
G.f.: 1/((1-x)*(1-x*y)*(1-x^2*y)). - Christian G. Bower, Jan 17 2006
a(n) = min(floor( 1/2 + sqrt(2*n)) - (2*n + round(sqrt(2*n)) - round(sqrt(2*n))^2)/2+1, (2*n + round(sqrt(2*n)) - round(sqrt(2*n))^2)/2). - Leonid Bedratyuk, Dec 13 2009

More terms from Larry Reeves (larryr(AT)acm.org), Nov 08 2000

Entry revised by N. J. A. Sloane, Dec 05 2006

A118323 (Greedy) Egyptian fraction expansion of Catalan constant.

Original entry on oeis.org

2, 3, 13, 176, 36543, 1394774578, 12493702893882521837, 265316559833226727589598741150947701321

Offset: 1

Eric W. Weisstein, Apr 23 2006

			Catalan constant = 1/2 + 1/3 + 1/13 + 1/176 + 1/36543 + ...

Amiram Eldar, Table of n, a(n) for n = 1..12
Mohammad K. Azarian, Sylvester's Sequence and the Infinite Egyptian Fraction Decomposition of 1, Problem 958, College Mathematics Journal, Vol. 42, No. 4, September 2011, p. 330.
Mohammad K. Azarian, Sylvester's Sequence and the Infinite Egyptian Fraction Decomposition of 1, Solution College Mathematics Journal, Vol. 43, No. 4, September 2012, pp. 340-342.
Eric Weisstein's World of Mathematics, Egyptian Fraction.
Index entries for sequences related to Egyptian fractions.

Cf. A006752, A104338, A014538, A153069, A153070, A054543. - Stuart Clary, Dec 17 2008

a = {}; k = N[Catalan, 1000]; Do[s = Ceiling[1/k]; AppendTo[a, s]; k = k - 1/s, {n, 1, 10}]; a (* Artur Jasinski, Sep 22 2008 *)

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

Stuart Clary, Dec 17 2008

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

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.

J. T. Groenman, Problem 1511, Crux Mathematicorum, Vol. 16, No. 2 (1990), p. 43; Solution to Problem 1511, by Beatriz Margolis, ibid., Vol. 17, No. 3 (1991), pp. 92-93.
Qing-Hu Hou and Zhi-Wei Sun, A q-analogue of the identity Sum_{k>=0}(-1)^k/(2k+1)^3 = Pi^3/32, arXiv:1808.04717 [math.CO], 2018.
Masato Kobayashi, Integral representations for zeta(3) with the inverse sine function, arXiv:2108.01247 [math.NT], 2021.
Richard J. Mathar, Table of Dirichlet L-series and prime zeta modulo functions for small moduli, arXiv:1008.2547 [math.NT], 2010-2015, section 2.2 entry L(m=4,r=2,s=3).
Michael Penn, An infinite tangent product, YouTube video, 2020.
Eric Weisstein's World of Mathematics, Dirichlet Beta Function.
Wikipedia, Dirichlet beta function.
Index entries for transcendental numbers.

Cf. A101455, A153072, A153073, A153074.
Cf. A233091, A251809, A331095.
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)).

nmax = 1000; First[ RealDigits[Pi^3/32, 10, nmax] ]

PARI
```
Pi^3/32 \\ Michel Marcus, Aug 15 2018
```

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)

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

A007341 Number of spanning trees in n X n grid.

Original entry on oeis.org

1, 4, 192, 100352, 557568000, 32565539635200, 19872369301840986112, 126231322912498539682594816, 8326627661691818545121844900397056, 5694319004079097795957215725765328371712000, 40325021721404118513276859513497679249183623593590784, 2954540993952788006228764987084443226815814190099484786032640000

Offset: 1

N. J. A. Sloane

Kreweras calls this the complexity of the n X n grid.
a(n) is the number of perfect mazes made from a grid of n X n cells. - Leroy Quet, Sep 08 2007
Also number of domino tilings of the (2n-1) X (2n-1) square with upper left corner removed.  For n=2 the 4 domino tilings of the 3 X 3 square with upper left corner removed are:
. ._. . ._. . ._. . ._.
.|__| .|__| .| | | .|___|
| |_| | | | | | ||| |_| |
||__| |||_| ||__| |_|_| - Alois P. Heinz, Apr 15 2011
Indeed, more is true. Let L denote the (2*n - 1) X (2*n - 1) square lattice graph with vertices (i,j), 1 <= i,j <= 2*n-1. Call a vertex (i,j) odd if both coordinates i and j are odd. Then there is a bijection between the set of spanning trees on the square n X n grid and the set of domino tilings of L with an odd boundary point removed. See Tzeng and Wu, 2002. This is a divisibility sequence, i.e., if n divides m then a(n) divides a(m). - Peter Bala, Apr 29 2014
Also, a(n) is the order of the sandpile group of the (n-1)X(n-1) grid graph. This is because the n X n grid is dual to (n-1)X(n-1) grid + sink vertex, and the latter is related to the sandpiles by the burning bijection. See Járai, Sec. 4.1, or Redig, Sec. 2.2. In M. F. Hasler's comment below, index n refers to the size of the grid underlying the sandpile. - Andrey Zabolotskiy, Mar 27 2018
From M. F. Hasler, Mar 07 2018: (Start)
The sandpile addition (+) of two n X n matrices is defined as the ordinary addition, followed by the topple-process in which each element larger than 3 is decreased by 4 and each of its von Neumann neighbors is increased by 1.
For any n, there is a neutral element e_n such that the set S(n) = { A in M_n({0..3}) | A (+) e_n = A } of matrices invariant under sandpile addition of e_n, forms a group, i.e., each element A in S(n) has an inverse A' in S(n) such that A (+) A' = e_n. (For n > 1, e_n cannot be the zero matrix O_n, because for this choice S(n) would include, e.g., the all 1's matrix 1_n which cannot have an inverse X such that 1_n (+) X = O_n. The element e_n is the unique nonzero matrix such that e_n (+) e_n = e_n.)
The present sequence lists the size of the abelian group (S(n), (+), e_n). See the example section for the e_n. The elements of S(2) are listed as A300006 and their inverses are listed as A300007. (End)

			From _M. F. Hasler_, Mar 07 2018: (Start)
For n = 1, there exists only one 0 X 0 matrix, e_0 = []; it is the neutral element of the singleton group S(0) = {[]}.
For n = 2, the sandpile addition is isomorphic to addition in Z/4Z, the neutral element is e_1 = [0] and we get the group S(1) isomorphic to (Z/4Z, +).
For n = 3, one finds that e_2 = [2,2;2,2] is the neutral element of the sandpile addition restricted to S(2), having 192 elements, listed in A300006.
For n = 4, one finds that e_3 = [2,1,2;1,0,1;2,1,2] is the neutral element of the sandpile addition restricted to S(3), having 100352 elements.
For n = 5, the neutral element is e_4 = [2,3,3,2; 3,2,2,3; 3,2,2,3; 2,3,3,2]. (End)

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

Alois P. Heinz, Table of n, a(n) for n = 1..45
Anakin Dey, Sam Ruggerio, and Melkior Ornik, Optimizing a Model-Agnostic Measure of Graph Counterdeceptiveness via Reattachment, arXiv:2311.15093 [math.OC], 2023. See p. 10.
Noah Doman, The Identity of the Abelian Sandpile Group, Bachelor Thesis, University of Groningen (Netherlands 2020).
Laura Florescu, Daniela Morar, David Perkinson, Nick Salter and Tianyuan Xu, Sandpiles and Dominos, Electronic Journal of Combinatorics, Volume 22, Issue 1 (2015), Paper #P1.66
Luis David Garcia-Puente and Brady Haran, Sandpiles, Numberphile video, on YouTube.com, Jan. 13, 2017
Antal A. Járai, Sandpile models, arXiv:1401.0354 [math.PR], 2014.
Germain Kreweras, Complexite et circuits Euleriens dans les sommes tensorielles de graphes, J. Combin. Theory, B 24 (1978), 202-212.
Lionel Levine and James Propp, What is... a sandpile?, Notices of the AMS, Volume 57 (2010), Number 8, 976-979.
F. Redig, Mathematical aspects of the abelian sandpile model (2005)
W.-J. Tzeng, F. Y. Wu, Spanning Trees on Hypercubic Lattices and Non-orientable Surfaces. arXiv:cond-mat/0001408v1 [cond-mat.stat-mech], Jan 2000.
W.-J. Tzeng and F. Y. Wu, Dimers on a simple-quartic net with a vacancy, arXiv:cond-mat/0203149v2 [cond-mat.stat-mech], Mar 2002.
Eric Weisstein's World of Mathematics, Grid Graph
Eric Weisstein's World of Mathematics, Spanning Tree
David B. Wilson, Local statistics of the abelian sandpile model (2014)
F. Y. Wu, Number of spanning trees on a lattice, J. Phys. A: Math. Gen., 10 (1977) no. 6, L113-L115.
Index to divisibility sequences

Main diagonal of A116469.
Cf. A300006, A300007, A300008, A300009; A256043, A256045.
Cf. A080690 (number of acyclic orientations), A080691 (number of spanning forests), A349718 (number of spanning trees, reduced for symmetry).

a:= n-> round(evalf(2^(n^2-1) /n^2 *mul(mul(`if`(j<>0 or k<>0, 2 -cos(Pi*j/n) -cos(Pi*k/n), 1), k=0..n-1), j=0..n-1), 15 +n*(n+1)/2)): seq(a(n), n=1..20);  # Alois P. Heinz, Apr 15 2011
# uses expression as a resultant
seq(resultant(simplify(ChebyshevU(n-1, x/2)), simplify(ChebyshevU(n-1, (4-x)/2)), x), n = 1 .. 24); # Peter Bala, Apr 29 2014

Table[2^((n-1)^2) Product[(2 - Cos[Pi i/n] - Cos[Pi j/n]), {i, 1, n-1}, {j, 1, n-1}], {n, 12}] // Round
Table[Resultant[ChebyshevU[n-1, x/2], ChebyshevU[n-1, (4-x)/2], x], {n, 1, 12}] (* Vaclav Kotesovec, Apr 15 2020 *)

{a(n) = polresultant( polchebyshev(n-1, 2, x/2), polchebyshev(n-1, 2, (4-x)/2) )}; /* Michael Somos, Aug 12 2017 */

a(n) = 2^(n^2-1) / n^2 * product_{n1=0..n-1, n2=0..n-1, n1 and n2 not both 0} (2 - cos(Pi*n1/n) - cos(Pi*n2/n) ). - Sharon Sela (sharonsela(AT)hotmail.com), Jun 04 2002
Equivalently, a(n) = Resultant( U(n-1,x/2), U(n-1,(4-x)/2) ), where U(n,x) is a Chebyshev polynomial of the second kind. - Peter Bala, Apr 29 2014
From Vaclav Kotesovec, Dec 30 2020: (Start)
a(n) ~ 2^(1/4) * Gamma(1/4) * exp(4*G*n^2/Pi) / (Pi^(3/4)*sqrt(n)*(1+sqrt(2))^(2*n)), where G is Catalan's constant A006752.
a(n) = n * 2^(n-1) * A007726(n)^2. (End)

More terms and better description from Roberto E. Martinez II, Jan 07 2002

A016814 a(n) = (4*n + 1)^2.

Original entry on oeis.org

1, 25, 81, 169, 289, 441, 625, 841, 1089, 1369, 1681, 2025, 2401, 2809, 3249, 3721, 4225, 4761, 5329, 5929, 6561, 7225, 7921, 8649, 9409, 10201, 11025, 11881, 12769, 13689, 14641, 15625, 16641, 17689, 18769, 19881, 21025, 22201, 23409, 24649, 25921, 27225, 28561, 29929

Offset: 0

N. J. A. Sloane

A bisection of A016754. Sequence arises from reading the line from 1, in the direction 1, 25, ..., in the square spiral whose vertices are the squares A000290. - Omar E. Pol, May 24 2008

G. C. Greubel, Table of n, a(n) for n = 0..10000 (terms 0..200 from Ivan Panchenko).
Leo Tavares, Illustration: Square trapeziums
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Sequences of the form (m*n+1)^2: A000012 (m=0), A000290 (m=1), A016754 (m=2), A016778 (m-3), this sequence (m=4), A016862 (m=5), A016922 (m=6), A016994 (m=7), A017078 (m=8), A017174 (m=9), A017282 (m=10), A017402 (m=11), A017534 (m=12), A134934 (m=14).
Cf. A001539, A006752, A016742, A016802, A016813.
Cf. A016826, A016838, A017126, A068467, A087197.
Cf. A272399, A014105.

[(4*n+1)^2: n in [0..40]]; // G. C. Greubel, Dec 28 2022

A016814:=n->(4*n+1)^2; seq(A016814(k), k=0..100); # Wesley Ivan Hurt, Nov 02 2013

(4*Range[0,40] +1)^2 (* or *) LinearRecurrence[{3,-3,1}, {1,25,81}, 40] (* Harvey P. Dale, Nov 20 2012 *)
Accumulate[32Range[0, 47] - 8] + 9 (* Alonso del Arte, Aug 19 2017 *)

a(n)=(4*n+1)^2 \\ Charles R Greathouse IV, Oct 07 2015

[(4*n+1)^2 for n in range(41)] # G. C. Greubel, Dec 28 2022

a(n) = a(n-1) + 32*n - 8, n > 0. - Vincenzo Librandi, Dec 15 2010
From George F. Johnson, Sep 28 2012: (Start)
G.f.: (1 + 22*x + 9*x^2)/(1 - x)^3.
a(n+1) = a(n) + 16 + 8*sqrt(a(n)).
a(n+1) = 2*a(n) - a(n-1) + 32 = 3*a(n) - 3*a(n-1) + a(n-2).
a(n-1)*a(n+1) = (a(n) - 16)^2 ; a(n+1) - a(n-1) = 16*sqrt(a(n)).
a(n) = A016754(2*n) = (A016813(n))^2. (End)
Sum_{n>=0} 1/a(n) = G/2 + Pi^2/16, where G is the Catalan constant (A006752). - Amiram Eldar, Jun 28 2020
Product_{n>=1} (1 - 1/a(n)) = 2*Gamma(5/4)^2/sqrt(Pi) = 2 * A068467^2 * A087197. - Amiram Eldar, Feb 01 2021
From G. C. Greubel, Dec 28 2022: (Start)
a(2*n) = A017078(n).
a(2*n+1) = A017126(n).
E.g.f.: (1 + 24*x + 16*x^2)*exp(x). (End)
a(n) = A272399(n+1) - A014105(n). - Leo Tavares, Dec 24 2023

cp's OEIS Frontend

A378025 Decimal expansion of 1/2 - log(2)/4 - G/Pi, where G = A006752.

Views

Author

Keywords

Examples

References

Crossrefs

Programs

Mathematica

Formula

A016754 Odd squares: a(n) = (2n+1)^2. Also centered octagonal numbers.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Mathematica

Maxima

PARI

Python

Formula

Extensions

A003881 Decimal expansion of Pi/4.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

SageMath

Formula

Extensions

A101455 a(n) = 0 for even n, a(n) = (-1)^((n-1)/2) for odd n. Periodic sequence 1,0,-1,0,...

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

PARI

Python

Formula

Extensions

A004003 Number of domino tilings (or dimer coverings) of a 2n X 2n square.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula