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.

Previous Showing 31-40 of 1164 results. Next

A168356 A000796(n-2) - A000796(n).

Original entry on oeis.org

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

Views

Author

Mark Lessel (mark(AT)lessel.us), Nov 23 2009

Keywords

Comments

Difference between digit n of pi and digit n+2.

Crossrefs

Programs

  • Mathematica
    Last[#]-First[#]&/@Partition[RealDigits[\[Pi],10,150][[1]], 3,1]  (* Harvey P. Dale, Mar 18 2011 *)

A180661 Decimal expansion of the constant whose continued fraction representation is [Pi^0; Pi^1, Pi^2, Pi^3, Pi^4, ...] where Pi is the ratio of a circle's circumference to its diameter (A000796) and the exponents cycle through all nonnegative integers.

Original entry on oeis.org

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

Views

Author

Jonathan D. B. Hodgson, Sep 15 2010

Keywords

Examples

			1.30839602...
		

Crossrefs

Extensions

More digits from R. J. Mathar, Sep 19 2010

A232191 A walk based on the digits of Pi (A000796).

Original entry on oeis.org

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

Views

Author

Philippe Deléham, Nov 20 2013

Keywords

Comments

Pi = 3.14159265358979323846....
Between 3 and 1 we place 2.
Between 1 and 4 we place 2 and 3.
Between 4 and 1 we place 3 and 2.
Between 1 and 5 we place 2, 3, and 4.
Between 5 and 9 we place 6, 7, and 8.
Between 9 and 2 we place 8, 7, 6, 5, 4 and 3.
Between 2 and 6 we place 3, 4 and 5.
Between 6 and 5 we place nothing.
Between 5 and 3 we place 4, and so on.
This gives:
3, 2, 1, 2, 3, 4, 3, 2, 1, 2, 3, 4, 5, 6, 7, 8, 9, 8, 7, 6, 5, 4, 3, 2, 3, 4, 5, 6, 5, 4, 3, 4, 5, ...
This could be called a walk (or promenade) on the digits of Pi.

Crossrefs

A237818 Decimal expansion of b^b^b^..., where b equals Pi-3 (A000796).

Original entry on oeis.org

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

Views

Author

Robert G. Wilson v, Mar 01 2014

Keywords

Examples

			=0.4307960410656323271949005495484538406836564717220423497579856062612692009443…
		

Crossrefs

Cf. A221566.

Programs

  • Mathematica
    b = N[Pi - 3, 128]; f[n_] := Nest[b^# &, b, n]; RealDigits[f[1500], 10, 111][[1]] (* Or *)
    RealDigits[ -LambertW[-Log[Pi - 3]]/Log[Pi - 3], 10, 105][[1]]

A002378 Oblong (or promic, pronic, or heteromecic) numbers: a(n) = n*(n+1).

Original entry on oeis.org

0, 2, 6, 12, 20, 30, 42, 56, 72, 90, 110, 132, 156, 182, 210, 240, 272, 306, 342, 380, 420, 462, 506, 552, 600, 650, 702, 756, 812, 870, 930, 992, 1056, 1122, 1190, 1260, 1332, 1406, 1482, 1560, 1640, 1722, 1806, 1892, 1980, 2070, 2162, 2256, 2352, 2450, 2550
Offset: 0

Views

Author

Keywords

Comments

4*a(n) + 1 are the odd squares A016754(n).
The word "pronic" (used by Dickson) is incorrect. - Michael Somos
According to the 2nd edition of Webster, the correct word is "promic". - R. K. Guy
a(n) is the number of minimal vectors in the root lattice A_n (see Conway and Sloane, p. 109).
Let M_n denote the n X n matrix M_n(i, j) = (i + j); then the characteristic polynomial of M_n is x^(n-2) * (x^2 - a(n)*x - A002415(n)). - Benoit Cloitre, Nov 09 2002
The greatest LCM of all pairs (j, k) for j < k <= n for n > 1. - Robert G. Wilson v, Jun 19 2004
First differences are a(n+1) - a(n) = 2*n + 2 = 2, 4, 6, ... (while first differences of the squares are (n+1)^2 - n^2 = 2*n + 1 = 1, 3, 5, ...). - Alexandre Wajnberg, Dec 29 2005
25 appended to these numbers corresponds to squares of numbers ending in 5 (i.e., to squares of A017329). - Lekraj Beedassy, Mar 24 2006
A rapid (mental) multiplication/factorization technique -- a generalization of Lekraj Beedassy's comment: For all bases b >= 2 and positive integers n, c, d, k with c + d = b^k, we have (n*b^k + c)*(n*b^k + d) = a(n)*b^(2*k) + c*d. Thus the last 2*k base-b digits of the product are exactly those of c*d -- including leading 0(s) as necessary -- with the preceding base-b digit(s) the same as a(n)'s. Examples: In decimal, 113*117 = 13221 (as n = 11, b = 10 = 3 + 7, k = 1, 3*7 = 21, and a(11) = 132); in octal, 61*67 = 5207 (52 is a(6) in octal). In particular, for even b = 2*m (m > 0) and c = d = m, such a product is a square of this type. Decimal factoring: 5609 is immediately seen to be 71*79. Likewise, 120099 = 301*399 (k = 2 here) and 99990000001996 = 9999002*9999998 (k = 3). - Rick L. Shepherd, Jul 24 2021
Number of circular binary words of length n + 1 having exactly one occurrence of 01. Example: a(2) = 6 because we have 001, 010, 011, 100, 101 and 110. Column 1 of A119462. - Emeric Deutsch, May 21 2006
The sequence of iterated square roots sqrt(N + sqrt(N + ...)) has for N = 1, 2, ... the limit (1 + sqrt(1 + 4*N))/2. For N = a(n) this limit is n + 1, n = 1, 2, .... For all other numbers N, N >= 1, this limit is not a natural number. Examples: n = 1, a(1) = 2: sqrt(2 + sqrt(2 + ...)) = 1 + 1 = 2; n = 2, a(2) = 6: sqrt(6 + sqrt(6 + ...)) = 1 + 2 = 3. - Wolfdieter Lang, May 05 2006
Nonsquare integers m divisible by ceiling(sqrt(m)), except for m = 0. - Max Alekseyev, Nov 27 2006
The number of off-diagonal elements of an (n + 1) X (n + 1) matrix. - Artur Jasinski, Jan 11 2007
a(n) is equal to the number of functions f:{1, 2} -> {1, 2, ..., n + 1} such that for a fixed x in {1, 2} and a fixed y in {1, 2, ..., n + 1} we have f(x) <> y. - Aleksandar M. Janjic and Milan Janjic, Mar 13 2007
Numbers m >= 0 such that round(sqrt(m+1)) - round(sqrt(m)) = 1. - Hieronymus Fischer, Aug 06 2007
Numbers m >= 0 such that ceiling(2*sqrt(m+1)) - 1 = 1 + floor(2*sqrt(m)). - Hieronymus Fischer, Aug 06 2007
Numbers m >= 0 such that fract(sqrt(m+1)) > 1/2 and fract(sqrt(m)) < 1/2 where fract(x) is the fractional part (fract(x) = x - floor(x), x >= 0). - Hieronymus Fischer, Aug 06 2007
X values of solutions to the equation 4*X^3 + X^2 = Y^2. To find Y values: b(n) = n(n+1)(2n+1). - Mohamed Bouhamida, Nov 06 2007
Nonvanishing diagonal of A132792, the infinitesimal Lah matrix, so "generalized factorials" composed of a(n) are given by the elements of the Lah matrix, unsigned A111596, e.g., a(1)*a(2)*a(3) / 3! = -A111596(4,1) = 24. - Tom Copeland, Nov 20 2007
If Y is a 2-subset of an n-set X then, for n >= 2, a(n-2) is the number of 2-subsets and 3-subsets of X having exactly one element in common with Y. - Milan Janjic, Dec 28 2007
a(n) coincides with the vertex of a parabola of even width in the Redheffer matrix, directed toward zero. An integer p is prime if and only if for all integer k, the parabola y = kx - x^2 has no integer solution with 1 < x < k when y = p; a(n) corresponds to odd k. - Reikku Kulon, Nov 30 2008
The third differences of certain values of the hypergeometric function 3F2 lead to the squares of the oblong numbers i.e., 3F2([1, n + 1, n + 1], [n + 2, n + 2], z = 1) - 3*3F2([1, n + 2, n + 2], [n + 3, n + 3], z = 1) + 3*3F2([1, n + 3, n + 3], [n + 4, n + 4], z = 1) - 3F2([1, n + 4, n + 4], [n + 5, n + 5], z = 1) = (1/((n+2)*(n+3)))^2 for n = -1, 0, 1, 2, ... . See also A162990. - Johannes W. Meijer, Jul 21 2009
Generalized factorials, [a.(n!)] = a(n)*a(n-1)*...*a(0) = A010790(n), with a(0) = 1 are related to A001263. - Tom Copeland, Sep 21 2011
For n > 1, a(n) is the number of functions f:{1, 2} -> {1, ..., n + 2} where f(1) > 1 and f(2) > 2. Note that there are n + 1 possible values for f(1) and n possible values for f(2). For example, a(3) = 12 since there are 12 functions f from {1, 2} to {1, 2, 3, 4, 5} with f(1) > 1 and f(2) > 2. - Dennis P. Walsh, Dec 24 2011
a(n) gives the number of (n + 1) X (n + 1) symmetric (0, 1)-matrices containing two ones (see [Cameron]). - L. Edson Jeffery, Feb 18 2012
a(n) is the number of positions of a domino in a rectangled triangular board with both legs equal to n + 1. - César Eliud Lozada, Sep 26 2012
a(n) is the number of ordered pairs (x, y) in [n+2] X [n+2] with |x-y| > 1. - Dennis P. Walsh, Nov 27 2012
a(n) is the number of injective functions from {1, 2} into {1, 2, ..., n + 1}. - Dennis P. Walsh, Nov 27 2012
a(n) is the sum of the positive differences of the partition parts of 2n + 2 into exactly two parts (see example). - Wesley Ivan Hurt, Jun 02 2013
a(n)/a(n-1) is asymptotic to e^(2/n). - Richard R. Forberg, Jun 22 2013
Number of positive roots in the root system of type D_{n + 1} (for n > 2). - Tom Edgar, Nov 05 2013
Number of roots in the root system of type A_n (for n > 0). - Tom Edgar, Nov 05 2013
From Felix P. Muga II, Mar 18 2014: (Start)
a(m), for m >= 1, are the only positive integer values t for which the Binet-de Moivre formula for the recurrence b(n) = b(n-1) + t*b(n-2) with b(0) = 0 and b(1) = 1 has a root of a square. PROOF (as suggested by Wolfdieter Lang, Mar 26 2014): The sqrt(1 + 4t) appearing in the zeros r1 and r2 of the characteristic equation is (a positive) integer for positive integer t precisely if 4t + 1 = (2m + 1)^2, that is t = a(m), m >= 1. Thus, the characteristic roots are integers: r1 = m + 1 and r2 = -m.
Let m > 1 be an integer. If b(n) = b(n-1) + a(m)*b(n-2), n >= 2, b(0) = 0, b(1) = 1, then lim_{n->oo} b(n+1)/b(n) = m + 1. (End)
Cf. A130534 for relations to colored forests, disposition of flags on flagpoles, and colorings of the vertices (chromatic polynomial) of the complete graphs (here simply K_2). - Tom Copeland, Apr 05 2014
The set of integers k for which k + sqrt(k + sqrt(k + sqrt(k + sqrt(k + ...) ... is an integer. - Leslie Koller, Apr 11 2014
a(n-1) is the largest number k such that (n*k)/(n+k) is an integer. - Derek Orr, May 22 2014
Number of ways to place a domino and a singleton on a strip of length n - 2. - Ralf Stephan, Jun 09 2014
With offset 1, this appears to give the maximal number of crossings between n nonconcentric circles of equal radius. - Felix Fröhlich, Jul 14 2014
For n > 1, the harmonic mean of the n values a(1) to a(n) is n + 1. The lowest infinite sequence of increasing positive integers whose cumulative harmonic mean is integral. - Ian Duff, Feb 01 2015
a(n) is the maximum number of queens of one color that can coexist without attacking one queen of the opponent's color on an (n+2) X (n+2) chessboard. The lone queen can be placed in any position on the perimeter of the board. - Bob Selcoe, Feb 07 2015
With a(0) = 1, a(n-1) is the smallest positive number not in the sequence such that Sum_{i = 1..n} 1/a(i-1) has a denominator equal to n. - Derek Orr, Jun 17 2015
The positive members of this sequence are a proper subsequence of the so-called 1-happy couple products A007969. See the W. Lang link there, eq. (4), with Y_0 = 1, with a table at the end. - Wolfdieter Lang, Sep 19 2015
For n > 0, a(n) is the reciprocal of the area bounded above by y = x^(n-1) and below by y = x^n for x in the interval [0, 1]. Summing all such areas visually demonstrates the formula below giving Sum_{n >= 1} 1/a(n) = 1. - Rick L. Shepherd, Oct 26 2015
It appears that, except for a(0) = 0, this is the set of positive integers n such that x*floor(x) = n has no solution. (For example, to get 3, take x = -3/2.) - Melvin Peralta, Apr 14 2016
If two independent real random variables, x and y, are distributed according to the same exponential distribution: pdf(x) = lambda * exp(-lambda * x), lambda > 0, then the probability that n - 1 <= x/y < n is given by 1/a(n). - Andres Cicuttin, Dec 03 2016
a(n) is equal to the sum of all possible differences between n different pairs of consecutive odd numbers (see example). - Miquel Cerda, Dec 04 2016
a(n+1) is the dimension of the space of vector fields in the plane with polynomial coefficients up to order n. - Martin Licht, Dec 04 2016
It appears that a(n) + 3 is the area of the largest possible pond in a square (A268311). - Craig Knecht, May 04 2017
Also the number of 3-cycles in the (n+3)-triangular honeycomb acute knight graph. - Eric W. Weisstein, Jul 27 2017
Also the Wiener index of the (n+2)-wheel graph. - Eric W. Weisstein, Sep 08 2017
The left edge of a Floyd's triangle that consists of even numbers: 0; 2, 4; 6, 8, 10; 12, 14, 16, 18; 20, 22, 24, 26, 28; ... giving 0, 2, 6, 12, 20, ... The right edge generates A028552. - Waldemar Puszkarz, Feb 02 2018
a(n+1) is the order of rowmotion on a poset obtained by adjoining a unique minimal (or maximal) element to a disjoint union of at least two chains of n elements. - Nick Mayers, Jun 01 2018
From Juhani Heino, Feb 05 2019: (Start)
For n > 0, 1/a(n) = n/(n+1) - (n-1)/n.
For example, 1/6 = 2/3 - 1/2; 1/12 = 3/4 - 2/3.
Corollary of this:
Take 1/2 pill.
Next day, take 1/6 pill. 1/2 + 1/6 = 2/3, so your daily average is 1/3.
Next day, take 1/12 pill. 2/3 + 1/12 = 3/4, so your daily average is 1/4.
And so on. (End)
From Bernard Schott, May 22 2020: (Start)
For an oblong number m >= 6 there exists a Euclidean division m = d*q + r with q < r < d which are in geometric progression, in this order, with a common integer ratio b. For b >= 2 and q >= 1, the Euclidean division is m = qb*(qb+1) = qb^2 * q + qb where (q, qb, qb^2) are in geometric progression.
Some examples with distinct ratios and quotients:
6 | 4 30 | 25 42 | 18
----- ----- -----
2 | 1 , 5 | 1 , 6 | 2 ,
and also:
42 | 12 420 | 100
----- -----
6 | 3 , 20 | 4 .
Some oblong numbers also satisfy a Euclidean division m = d*q + r with q < r < d that are in geometric progression in this order but with a common noninteger ratio b > 1 (see A335064). (End)
For n >= 1, the continued fraction expansion of sqrt(a(n)) is [n; {2, 2n}]. For n=1, this collapses to [1; {2}]. - Magus K. Chu, Sep 09 2022
a(n-2) is the maximum irregularity over all trees with n vertices. The extremal graphs are stars. (The irregularity of a graph is the sum of the differences between the degrees over all edges of the graph.) - Allan Bickle, May 29 2023
For n > 0, number of diagonals in a regular 2*(n+1)-gon that are not parallel to any edge (cf. A367204). - Paolo Xausa, Mar 30 2024
a(n-1) is the maximum Zagreb index over all trees with n vertices. The extremal graphs are stars. (The Zagreb index of a graph is the sum of the squares of the degrees over all vertices of the graph.) - Allan Bickle, Apr 11 2024
For n >= 1, a(n) is the determinant of the distance matrix of a cycle graph on 2*n + 1 vertices (if the length of the cycle is even such a determinant is zero). - Miquel A. Fiol, Aug 20 2024
For n > 1, the continued fraction expansion of sqrt(16*a(n)) is [2n+1; {1, 2n-1, 1, 8n+2}]. - Magus K. Chu, Nov 20 2024
For n>=2, a(n) is the number of faces on a n+1-zone rhombic zonohedron. Each pair of a collection of great circles on a sphere intersects at two points, so there are 2*binomial(n+1,2) intersections. The dual of the implied polyhedron is a rhombic zonohedron, its faces corresponding to the intersections. - Shel Kaphan, Aug 12 2025

Examples

			a(3) = 12, since 2(3)+2 = 8 has 4 partitions with exactly two parts: (7,1), (6,2), (5,3), (4,4). Taking the positive differences of the parts in each partition and adding, we get: 6 + 4 + 2 + 0 = 12. - _Wesley Ivan Hurt_, Jun 02 2013
G.f. = 2*x + 6*x^2 + 12*x^3 + 20*x^4 + 30*x^5 + 42*x^6 + 56*x^7 + ... - _Michael Somos_, May 22 2014
From _Miquel Cerda_, Dec 04 2016: (Start)
a(1) = 2, since 45-43 = 2;
a(2) = 6, since 47-45 = 2 and 47-43 = 4, then 2+4 = 6;
a(3) = 12, since 49-47 = 2, 49-45 = 4, and 49-43 = 6, then 2+4+6 = 12. (End)
		

References

  • W. W. Berman and D. E. Smith, A Brief History of Mathematics, 1910, Open Court, page 67.
  • J. H. Conway and R. K. Guy, The Book of Numbers, 1996, p. 34.
  • J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag.
  • L. E. Dickson, History of the Theory of Numbers, Vol. 1: Divisibility and Primality. New York: Chelsea, p. 357, 1952.
  • L. E. Dickson, History of the Theory of Numbers, Vol. 2: Diophantine Analysis. New York: Chelsea, pp. 6, 232-233, 350 and 407, 1952.
  • H. Eves, An Introduction to the History of Mathematics, revised, Holt, Rinehart and Winston, 1964, page 72.
  • Nicomachus of Gerasa, Introduction to Arithmetic, translation by Martin Luther D'Ooge, Ann Arbor, University of Michigan Press, 1938, p. 254.
  • Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §8.6 Figurate Numbers, p. 291.
  • Granino A. Korn and Theresa M. Korn, Mathematical Handbook for Scientists and Engineers, McGraw-Hill Book Company, New York (1968), pp. 980-981.
  • C. S. Ogilvy and J. T. Anderson, Excursions in Number Theory, Oxford University Press, 1966, pp. 61-62.
  • Alfred S. Posamentier, Math Charmers, Tantalizing Tidbits for the Mind, Prometheus Books, NY, 2003, pages 54-55.
  • 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).
  • F. J. Swetz, From Five Fingers to Infinity, Open Court, 1994, p. 219.
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 2-6.

Crossrefs

Partial sums of A005843 (even numbers). Twice triangular numbers (A000217).
1/beta(n, 2) in A061928.
A036689 and A036690 are subsequences. Cf. numbers of the form n*(n*k-k+4)/2 listed in A226488. - Bruno Berselli, Jun 10 2013
Row n=2 of A185651.
Cf. A007745, A169810, A213541, A005369 (characteristic function).
Cf. A281026. - Bruno Berselli, Jan 16 2017
Cf. A045943 (4-cycles in triangular honeycomb acute knight graph), A028896 (5-cycles), A152773 (6-cycles).
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.
A335064 is a subsequence.
Second column of A003506.
Cf. A002378, A046092, A028896 (irregularities of maximal k-degenerate graphs).
Cf. A347213 (Dgf at s=4).
Cf. A002378, A152811, A371912 (Zagreb indices of maximal k-degenerate graphs).

Programs

Formula

G.f.: 2*x/(1-x)^3. - Simon Plouffe in his 1992 dissertation.
a(n) = a(n-1) + 2*n, a(0) = 0.
Sum_{n >= 1} a(n) = n*(n+1)*(n+2)/3 (cf. A007290, partial sums).
Sum_{n >= 1} 1/a(n) = 1. (Cf. Tijdeman)
Sum_{n >= 1} (-1)^(n+1)/a(n) = log(4) - 1 = A016627 - 1 [Jolley eq (235)].
1 = 1/2 + Sum_{n >= 1} 1/(2*a(n)) = 1/2 + 1/4 + 1/12 + 1/24 + 1/40 + 1/60 + ... with partial sums: 1/2, 3/4, 5/6, 7/8, 9/10, 11/12, 13/14, ... - Gary W. Adamson, Jun 16 2003
a(n)*a(n+1) = a(n*(n+2)); e.g., a(3)*a(4) = 12*20 = 240 = a(3*5). - Charlie Marion, Dec 29 2003
Sum_{k = 1..n} 1/a(k) = n/(n+1). - Robert G. Wilson v, Feb 04 2005
a(n) = A046092(n)/2. - Zerinvary Lajos, Jan 08 2006
Log 2 = Sum_{n >= 0} 1/a(2n+1) = 1/2 + 1/12 + 1/30 + 1/56 + 1/90 + ... = (1 - 1/2) + (1/3 - 1/4) + (1/5 - 1/6) + (1/7 - 1/8) + ... = Sum_{n >= 0} (-1)^n/(n+1) = A002162. - Gary W. Adamson, Jun 22 2003
a(n) = A110660(2*n). - N. J. A. Sloane, Sep 21 2005
a(n-1) = n^2 - n = A000290(n) - A000027(n) for n >= 1. a(n) is the inverse (frequency distribution) sequence of A000194(n). - Mohammad K. Azarian, Jul 26 2007
(2, 6, 12, 20, 30, ...) = binomial transform of (2, 4, 2). - Gary W. Adamson, Nov 28 2007
a(n) = 2*Sum_{i=0..n} i = 2*A000217(n). - Artur Jasinski, Jan 09 2007, and Omar E. Pol, May 14 2008
a(n) = A006503(n) - A000292(n). - Reinhard Zumkeller, Sep 24 2008
a(n) = A061037(4*n) = (n+1/2)^2 - 1/4 = ((2n+1)^2 - 1)/4 = (A005408(n)^2 - 1)/4. - Paul Curtz, Oct 03 2008 and Klaus Purath, Jan 13 2022
a(0) = 0, a(n) = a(n-1) + 1 + floor(x), where x is the minimal positive solution to fract(sqrt(a(n-1) + 1 + x)) = 1/2. - Hieronymus Fischer, Dec 31 2008
E.g.f.: (x+2)*x*exp(x). - Geoffrey Critzer, Feb 06 2009
Product_{i >= 2} (1-1/a(i)) = -2*sin(Pi*A001622)/Pi = -2*sin(A094886)/A000796 = 2*A146481. - R. J. Mathar, Mar 12 2009, Mar 15 2009
E.g.f.: ((-x+1)*log(-x+1)+x)/x^2 also Integral_{x = 0..1} ((-x+1)*log(-x+1) + x)/x^2 = zeta(2) - 1. - Stephen Crowley, Jul 11 2009
a(A007018(n)) = A007018(n+1), i.e., A007018(n+1) = A007018(n)-th oblong numbers. - Jaroslav Krizek, Sep 13 2009
a(n) = floor((n + 1/2)^2). a(n) = A035608(n) + A004526(n+1). - Reinhard Zumkeller, Jan 27 2010
a(n) = 2*(2*A006578(n) - A035608(n)). - Reinhard Zumkeller, Feb 07 2010
a(n-1) = floor(n^5/(n^3 + n^2 + 1)). - Gary Detlefs, Feb 11 2010
For n > 1: a(n) = A173333(n+1, n-1). - Reinhard Zumkeller, Feb 19 2010
a(n) = A004202(A000217(n)). - Reinhard Zumkeller, Feb 12 2011
a(n) = A188652(2*n+1) + 1. - Reinhard Zumkeller, Apr 13 2011
For n > 0 a(n) = 1/(Integral_{x=0..Pi/2} 2*(sin(x))^(2*n-1)*(cos(x))^3). - Francesco Daddi, Aug 02 2011
a(n) = A002061(n+1) - 1. - Omar E. Pol, Oct 03 2011
a(0) = 0, a(n) = A005408(A034856(n)) - A005408(n-1). - Ivan N. Ianakiev, Dec 06 2012
a(n) = A005408(A000096(n)) - A005408(n). - Ivan N. Ianakiev, Dec 07 2012
a(n) = A001318(n) + A085787(n). - Omar E. Pol, Jan 11 2013
Sum_{n >= 1} 1/(a(n))^(2s) = Sum_{t = 1..2*s} binomial(4*s - t - 1, 2*s - 1) * ( (1 + (-1)^t)*zeta(t) - 1). See Arxiv:1301.6293. - R. J. Mathar, Feb 03 2013
a(n)^2 + a(n+1)^2 = 2 * a((n+1)^2), for n > 0. - Ivan N. Ianakiev, Apr 08 2013
a(n) = floor(n^2 * e^(1/n)) and a(n-1) = floor(n^2 / e^(1/n)). - Richard R. Forberg, Jun 22 2013
a(n) = 2*C(n+1, 2), for n >= 0. - Felix P. Muga II, Mar 11 2014
A005369(a(n)) = 1. - Reinhard Zumkeller, Jul 05 2014
Binomial transform of [0, 2, 2, 0, 0, 0, ...]. - Alois P. Heinz, Mar 10 2015
a(2n) = A002943(n) for n >= 0, a(2n-1) = A002939(n) for n >= 1. - M. F. Hasler, Oct 11 2015
For n > 0, a(n) = 1/(Integral_{x=0..1} (x^(n-1) - x^n) dx). - Rick L. Shepherd, Oct 26 2015
a(n) = A005902(n) - A007588(n). - Peter M. Chema, Jan 09 2016
For n > 0, a(n) = lim_{m -> oo} (1/m)*1/(Sum_{i=m*n..m*(n+1)} 1/i^2), with error of ~1/m. - Richard R. Forberg, Jul 27 2016
From Ilya Gutkovskiy, Jul 28 2016: (Start)
Dirichlet g.f.: zeta(s-2) + zeta(s-1).
Convolution of nonnegative integers (A001477) and constant sequence (A007395).
Sum_{n >= 0} a(n)/n! = 3*exp(1). (End)
From Charlie Marion, Mar 06 2020: (Start)
a(n)*a(n+2k-1) + (n+k)^2 = ((2n+1)*k + n^2)^2.
a(n)*a(n+2k) + k^2 = ((2n+1)*k + a(n))^2. (End)
Product_{n>=1} (1 + 1/a(n)) = cosh(sqrt(3)*Pi/2)/Pi. - Amiram Eldar, Jan 20 2021
A generalization of the Dec 29 2003 formula, a(n)*a(n+1) = a(n*(n+2)), follows. a(n)*a(n+k) = a(n*(n+k+1)) + (k-1)*n*(n+k+1). - Charlie Marion, Jan 02 2023
a(n) = A016742(n) - A049450(n). - Leo Tavares, Mar 15 2025

Extensions

Additional comments from Michael Somos
Comment and cross-reference added by Christopher Hunt Gribble, Oct 13 2009

A023900 Dirichlet inverse of Euler totient function (A000010).

Original entry on oeis.org

1, -1, -2, -1, -4, 2, -6, -1, -2, 4, -10, 2, -12, 6, 8, -1, -16, 2, -18, 4, 12, 10, -22, 2, -4, 12, -2, 6, -28, -8, -30, -1, 20, 16, 24, 2, -36, 18, 24, 4, -40, -12, -42, 10, 8, 22, -46, 2, -6, 4, 32, 12, -52, 2, 40, 6, 36, 28, -58, -8, -60, 30, 12, -1, 48, -20, -66, 16, 44, -24, -70, 2, -72, 36, 8, 18, 60, -24, -78, 4, -2
Offset: 1

Views

Author

Keywords

Comments

Also called reciprocity balance of n.
Apart from different signs, same as Sum_{d divides n} core(d)*mu(n/d), where core(d) (A007913) is the squarefree part of d. - Benoit Cloitre, Apr 06 2002
Main diagonal of A191898. - Mats Granvik, Jun 19 2011

Examples

			x - x^2 - 2*x^3 - x^4 - 4*x^5 + 2*x^6 - 6*x^7 - x^8 - 2*x^9 + 4*x^10 - ...
		

References

  • T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 37.
  • D. M. Burton, Elementary Number Theory, Allyn and Bacon Inc. Boston, MA, 1976, p. 125.

Crossrefs

Moebius transform is A055615.
Cf. A027748, A173557 (gives the absolute values), A295876.
Cf. A253905 (Dgf at s=3).

Programs

  • Haskell
    a023900 1 = 1
    a023900 n = product $ map (1 -) $ a027748_row n
    -- Reinhard Zumkeller, Jun 01 2015
    
  • Maple
    A023900 := n -> mul(1-i,i=numtheory[factorset](n)); # Peter Luschny, Oct 26 2010
  • Mathematica
    a[ n_] := If[ n < 1, 0, Sum[ d MoebiusMu @ d, { d, Divisors[n]}]] (* Michael Somos, Jul 18 2011 *)
    Array[ Function[ n, 1/Plus @@ Map[ #*MoebiusMu[ # ]/EulerPhi[ # ]&, Divisors[ n ] ] ], 90 ]
    nmax = 81; Drop[ CoefficientList[ Series[ Sum[ MoebiusMu[k] k x^k/(1 - x^k), {k, 1, nmax} ], {x, 0, nmax} ], x ], 1 ] (* Stuart Clary, Apr 15 2006 *)
    t[n_, 1] = 1; t[1, k_] = 1; t[n_, k_] :=  t[n, k] = If[n < k, If[n > 1 && k > 1, Sum[-t[k - i, n], {i, 1, n - 1}], 0], If[n > 1 && k > 1, Sum[-t[n - i, k], {i, 1, k - 1}], 0]]; Table[t[n, n], {n, 36}] (* Mats Granvik, Robert G. Wilson v, Jun 25 2011 *)
    Table[DivisorSum[m, # MoebiusMu[#] &], {m, 90}] (* Jan Mangaldan, Mar 15 2013 *)
    f[p_, e_] := (1 - p); a[1] = 1; a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100] (* Amiram Eldar, Oct 14 2020 *)
  • PARI
    {a(n) = direuler( p=2, n, (1 - p*X) / (1 - X))[n]}
    
  • PARI
    {a(n) = if( n<1, 0, sumdiv( n, d, d * moebius(d)))} /* Michael Somos, Jul 18 2011 */
    
  • PARI
    a(n)=sumdivmult(n,d, d*moebius(d)) \\ Charles R Greathouse IV, Sep 09 2014
    
  • Python
    from sympy import divisors, mobius
    def a(n): return sum([d*mobius(d) for d in divisors(n)]) # Indranil Ghosh, Apr 29 2017
    
  • Python
    from math import prod
    from sympy import primefactors
    def A023900(n): return prod(1-p for p in primefactors(n)) # Chai Wah Wu, Sep 08 2023
    
  • Scheme
    ;; With memoization-macro definec.
    (definec (A023900 n) (if (= 1 n) 1 (* (- 1 (A020639 n)) (A023900 (A028234 n))))) ;; Antti Karttunen, Nov 28 2017

Formula

a(n) = Sum_{ d divides n } d*mu(d) = Product_{p|n} (1-p).
a(n) = 1 / (Sum_{ d divides n } mu(d)*d/phi(d)).
Dirichlet g.f.: zeta(s)/zeta(s-1). - Michael Somos, Jun 04 2000
a(n+1) = det(n+1)/det(n) where det(n) is the determinant of the n X n matrix M_(i, j) = i/gcd(i, j) = lcm(i, j)/j. - Benoit Cloitre, Aug 19 2003
a(n) = phi(n)*moebius(A007947(n))*A007947(n)/n. Logarithmic g.f.: Sum_{n >= 1} a(n)*x^n/n = log(F(x)) where F(x) is the g.f. of A117209 and satisfies: 1/(1-x) = Product_{n >= 1} F(x^n). - Paul D. Hanna, Mar 03 2006
G.f.: A(x) = Sum_{k >= 1} mu(k) k x^k/(1 - x^k) where mu(k) is the Moebius (Mobius) function, A008683. - Stuart Clary, Apr 15 2006
G.f.: A(x) is x times the logarithmic derivative of A117209(x). - Stuart Clary, Apr 15 2006
Row sums of triangle A134842. - Gary W. Adamson, Nov 12 2007
G.f.: x/(1-x) = Sum_{n >= 1} a(n)*x^n/(1-x^n)^2. - Paul D. Hanna, Aug 16 2008
a(n) = phi(rad(n)) *(-1)^omega(n) = A000010(A007947(n)) *(-1)^A001221(n). - Enrique Pérez Herrero, Aug 24 2010
a(n) = Product_{i = 2..n} (1-i)^( (pi(i)-pi(i-1)) * floor( (cos(n*Pi/i))^2 ) ), where pi = A000720, Pi = A000796. - Wesley Ivan Hurt, May 24 2013
a(n) = -limit of zeta(s)*(Sum_{d divides n} moebius(d)/exp(d)^(s-1)) as s->1 for n>1. - Mats Granvik, Jul 31 2013
a(n) = Sum_{d divides n} mu(d)*rad(d), where rad is A007947. - Enrique Pérez Herrero, May 29 2014
Conjecture for n>1: Let n = 2^(A007814(n))*m = 2^(ruler(n))*odd_part(n), where m = A000265(n), then a(n) = (-1)^(m=n)*(0+Sum_{i=1..m and gcd(i,m)=1} (4*min(i,m-i)-m)) = (-1)^(m1} (4*min(i,m-i)-m)). - I. V. Serov, May 02 2017
a(n) = (-1)^A001221(n) * A173557(n). - R. J. Mathar, Nov 02 2017
a(1) = 1; for n > 1, a(n) = (1-A020639(n)) * a(A028234(n)), because multiplicative with a(p^e) = (1-p). - Antti Karttunen, Nov 28 2017
a(n) = 1 - Sum_{d|n, d > 1} d*a(n/d). - Ilya Gutkovskiy, Apr 26 2019
From Richard L. Ollerton, May 07 2021: (Start)
For n>1, Sum_{k=1..n} a(gcd(n,k)) = 0.
For n>1, Sum_{k=1..n} a(n/gcd(n,k))*phi(gcd(n,k))/phi(n/gcd(n,k)) = 0. (End)
a(n) = rad(n)*(-1)^omega(n)*phi(n)/n = A062953(n)*A000010(n)/n. - Amrit Awasthi, Jan 30 2022
a(n) = mu(n)*phi(n) = A008683(n)*A000010(n) whenever n is squarefree. - Amrit Awasthi, Feb 03 2022
From Peter Bala, Jan 24 2024: (Start)
a(n) = Sum_{d divides n} core(d)*mu(d). Cf. Comment by Benoit Cloitre, dated Apr 06 2002.
a(n) = Sum_{d|n, e|n} n/gcd(d, e) * mu(n/d) * mu(n/e) (the sum is a multiplicative function of n by Tóth, and takes the value 1 - p for n = p^e, a prime power). (End)
From Peter Bala, Feb 01 2024: (Start)
G.f. Sum_{n >= 1} (2*n-1)*moebius(2*n-1)*x^(2*n-1)/(1 + x^(2n-1)).
a(n) = (-1)^(n+1) * Sum_{d divides n, d odd} d*moebius(d). (End)

A019669 Decimal expansion of Pi/2.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

With offset 2, decimal expansion of 5*Pi. - Omar E. Pol, Oct 03 2013
Decimal expansion of the number of radians in a quadrant. - John W. Nicholson, Oct 07 2013
Not the same as A085679. First differing term occurs at 10^-49, as list -49, or 51st counting term (a(-49)= 5 and A085679(-49) = 4). - John W. Nicholson, Oct 07 2013
5*Pi is also the surface area of a sphere whose diameter equals the square root of 5. More generally x*Pi is also the surface area of a sphere whose diameter equals the square root of x. - Omar E. Pol, Dec 22 2013
Pi/2 is also the radius of a sphere whose surface area equals the volume of the circumscribed cube. - Omar E. Pol, Dec 27 2013

Examples

			Pi/2 = 1.570796326794896619231321691639751442098584699...
5*Pi = 15.70796326794896619231321691639751442098584699...
		

References

  • Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, 2003, Sections 1.4.1 and 1.4.2, pp. 20-21.

Crossrefs

Cf. A053300 (continued fraction), A060294 (2/Pi).
Cf. A000796, A019692, A122952, A019694 (Pi through 4*Pi), A106854.

Programs

Formula

Pi/2 = log(i)/i, where i = sqrt(-1). - Eric Desbiaux, Jun 27 2009
Pi/2 = Product_{n>=1} (n/(n+1))^((-1)^n) = 2 * 2/3 * 4/3 * 4/5 * 6/5 * 6/7 * 8/7 * 8/9 * 10/9 * ... (Wallis formula). - William Keith and Alonso del Arte, Jun 24 2012
Equals Sum_{k>1} 2^k/binomial(2*k,k). - Bruno Berselli, Sep 11 2015
The previous result is the particular case n = 1 of the more general identity: Pi/2 = 4^(n-1) * n!/(2*n)! * Sum_{k >= 2} 2^(k+1)*(k + n - 1)!*(k + 2*n - 2)!/(2*k + 2*n - 2)! valid for n = 0,1,2,... . - Peter Bala, Oct 26 2016
Pi/2 = Product_{n>=1} (4*n^2)/(4*n^2-1). - Fred Daniel Kline, Oct 29 2016
Pi/2 = lim_{n->oo} F(2^(n+3))/2, with one half of the area of a regular 2^(n+3)-gon, for n >= 0, inscribed in the unit circle, written as iterated square roots of 2 as F(2^(n+3))/2 = 2^n*sqrt(2 + sq2(n)), with sq2(n) = sqrt(2 + sq2(n-1)), n >= 1, with input sq2(0) = 0 (2 appears n times in sq2(n)). Viète's infinite product formula works with the partial product F(2^(n+2))/2 = Product_{j=1..n} (2/sq2(j)), n >= 1, which corresponds to the above given formula. - Wolfdieter Lang, Jul 06 2018
Pi/2 = Integral_{x = 0..oo} sin(x)^2/x^2 dx = 1/2 + Sum_{n >= 1} sin(n)^2/n^2, by the Abel-Plana formula. - Peter Bala, Nov 05 2019
From Amiram Eldar, Aug 15 2020: (Start)
Equals Sum_{k>=0} k!/(2*k + 1)!!.
Equals Sum_{k>=0} (-1)^k/(k + 1/2).
Equals Integral_{x=0..oo} 1/(x^2 + 1) dx.
Equals Integral_{x=0..oo} sin(x)/x dx.
Equals Integral_{x=0..oo} exp(x/2)/(exp(x) + 1) dx.
Equals Product_{p prime > 2} p/(p + (-1)^((p-1)/2)). (End)
Pi/2 = Integral_{x = 0..oo} 1/(1 - x^2 + x^4) dx = (1 + 2/3 + 1/5) - (1/7 + 2/9 + 1/11) + (1/13 + 2/15 + 1/17) - .... - Peter Bala, Jul 22 2022
Equals arcsin(9/10) + sqrt(19)*Sum_{k >= 1} A106854(k-1)/(k*10^k) (see Bailey and Crandall, 2001). - Paolo Xausa, Jul 15 2024
Equals 2F1(1/2,1/2 ; 3/2; 1). - R. J. Mathar, Aug 20 2024
Pi/2 = [1;1,1/2,1/3,...,1/n,...] by Wallis's approximation. - Thomas Ordowski, Oct 19 2024
From Stefano Spezia, Oct 21 2024: (Start)
Equals Sum_{k>=0} 2^k/((2*k + 1)*binomial(2*k,k)) (see Finch).
Equals Limit_{n->oo} 2^(4*n)/((2*n + 1)*binomial(2*n,n)^2) (see Finch). (End)
Equals Integral_{x=-oo..oo} sech((2*x^3 + x^2 - 5*x)/(x^2 - 1)) dx. - Kritsada Moomuang, May 29 2025

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

A049541 Decimal expansion of 1/Pi.

Original entry on oeis.org

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

Views

Author

N. J. A. Sloane, Dec 11 1999

Keywords

Comments

The ratio of the volume of a regular octahedron to the volume of the circumscribed sphere (which has circumradius a*sqrt(2)/2 = a*A010503, where a is the octahedron's edge length; see MathWorld link). For similar ratios for other Platonic solids, see A165922, A165952, A165953 and A165954. - Rick L. Shepherd, Oct 01 2009
Corresponds to a gauge point marked "M" on slide rule calculating devices in the 20th century. The Pickworth reference notes its use in calculating the area of the curved surface of a cylinder. - Peter Munn, Aug 14 2020

Examples

			0.3183098861837906715377675267450287240689192914809128974953...
		

References

  • J.-P. Delahaye, Pi - die Story (German translation), Birkhäuser, 1999 Baasel, p. 245. French original: Le fascinant nombre Pi, Pour la Science, Paris, 1997.
  • C. N. Pickworth, The Slide Rule, 24th Ed., Pitman, London, 1945, p. 53, Gauge Points.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987, p. 27.

Crossrefs

Programs

Formula

Equals (1/(12-16*A002162))*Sum_{n>=0} A002894(n)*H(n)/(A001025(n) * A016754(n-1)), where H(n) denotes the n-th harmonic number. - John M. Campbell, Aug 28 2016
1/Pi = Sum_{m>=0} binomial(2*m, m)^3 * (42*m+5)/(2^(12*m+4)), Ramanujan, from the J.-P. Delahaye reference. - Wolfdieter Lang, Sep 18 2018; corrected by Bernard Schott, Mar 26 2020
1/Pi = 12*Sum_{n >= 0} (-1)^n*((6*n)!/(n!^3*(3*n)!))*(13591409 + 545140134*n)/640320^(3*n + 3/2) [Chudnovsky]. - Sanjar Abrarov, Mar 31 2020
1/Pi = (sqrt(8)/9801) * Sum_{n >= 0} ((4*n)!/((n!)^4)) * (26390*n + 1103)/(396^(4*n)) [Ramanujan, 1914]. - Bernard Schott, Mar 26 2020
Equal Sum_{k>=2} tan(Pi/2^k)/2^k. - Amiram Eldar, Aug 05 2020
Floor((3/8)*Sum_{n>=1} sigma[3](n)*n/exp(Pi*n/(10^((1/5)*k+(1/5))))) mod 10, will give the k-th digit of 1/Pi. - Simon Plouffe, Dec 19 2023

A002161 Decimal expansion of square root of Pi.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

Also Gamma(1/2). - Franklin T. Adams-Watters, Apr 07 2006
The integral of the Gaussian function exp(-x^2) over the real line. - Richard Chapling (r.chappers(AT)gmail.com), Jun 05 2008
Also equals the average distance between two points in two dimensions where coordinates are independent normally distributed random variables with mean 0 and variance 1. - Jean-François Alcover, Oct 31 2014, after Steven Finch
Also diameter of a sphere whose surface area equals Pi^2. More generally, the square root of x is also the diameter of a sphere whose surface area equals x*Pi. - Omar E. Pol, Nov 11 2018
Convergents of continued fractions: 7/4, 16/9, 23/13, 39/22, 257/145, 296/167, 8545/4821, ... - R. J. Mathar, Jan 29 2025

Examples

			1.7724538509055160272981674833411451827975494561223871282138...
		

References

  • George Boros and Victor H. Moll, Irresistible integrals, Cambridge University Press (2006), p. 190.
  • Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, 2003, Section 1.5.4, p. 33.
  • W. E. Mansell, Tables of Natural and Common Logarithms. Royal Society Mathematical Tables, Vol. 8, Cambridge Univ. Press, 1964, p. XVIII.
  • 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).
  • Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 43, page 413.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 40.

Crossrefs

Cf. decimal expansions of Gamma(1/k): A073005 (k=3), A068466 (k=4), A175380 (k=5), A175379 (k=6), A220086 (k=7), A203142 (k=8).

Programs

  • Magma
    R:= RealField(100); Sqrt(Pi(R));  // G. C. Greubel, Mar 10 2018
  • Maple
    evalf(sqrt(Pi),120); # Muniru A Asiru, Nov 11 2018
  • Mathematica
    RealDigits[N[Sqrt[Pi], 120]][[1]] (* Richard Chapling (r.chappers(AT)gmail.com), Jun 05 2008 *)
  • PARI
    default(realprecision, 20080); x=sqrt(Pi); for (n=1, 20000, d=floor(x); x=(x-d)*10; write("b002161.txt", n, " ", d)); \\ Harry J. Smith, May 01 2009
    

Formula

Equals (1/2) * Sum_{n>=0} ((-1)^n * (4*n+1) * (1/8)^(n+1) * (2^(n+1))^3 * Gamma(n+1/2)^3 / Gamma(n+1)^3). - Alexander R. Povolotsky, Mar 25 2013
Equals Integral_{x=0..1} 1/sqrt(-log(x)) dx. - Jean-François Alcover, Apr 29 2013
Equals Sum_{k>=0} (k+1/2)!/(k+2)!. - Amiram Eldar, Jun 19 2023
Equals Integral_{x=0..oo} exp(-x)/sqrt(x) dx. - Michal Paulovic, Sep 24 2023
Equals Integral_{x=0..oo} 4/(exp(x^2)*(2*x^2 + 1)^2) dx. - Kritsada Moomuang, Jun 05 2025

Extensions

More terms from Franklin T. Adams-Watters, Apr 07 2006
Previous Showing 31-40 of 1164 results. Next