cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 29 results. Next

A007814 Exponent of highest power of 2 dividing n, a.k.a. the binary carry sequence, the ruler sequence, or the 2-adic valuation of n.

Original entry on oeis.org

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

Views

Author

John Tromp, Dec 11 1996

Keywords

Comments

This sequence is an exception to my usual rule that when every other term of a sequence is 0 then those 0's should be omitted. In this case we would get A001511. - N. J. A. Sloane
To construct the sequence: start with 0,1, concatenate to get 0,1,0,1. Add + 1 to last term gives 0,1,0,2. Concatenate those 4 terms to get 0,1,0,2,0,1,0,2. Add + 1 to last term etc. - Benoit Cloitre, Mar 06 2003
The sequence is invariant under the following two transformations: increment every element by one (1, 2, 1, 3, 1, 2, 1, 4, ...), put a zero in front and between adjacent elements (0, 1, 0, 2, 0, 1, 0, 3, 0, 1, 0, 2, 0, 1, 0, 4, ...). The intermediate result is A001511. - Ralf Hinze (ralf(AT)informatik.uni-bonn.de), Aug 26 2003
Fixed point of the morphism 0->01, 1->02, 2->03, 3->04, ..., n->0(n+1), ..., starting from a(1) = 0. - Philippe Deléham, Mar 15 2004
Fixed point of the morphism 0->010, 1->2, 2->3, ..., n->(n+1), .... - Joerg Arndt, Apr 29 2014
a(n) is also the number of times to repeat a step on an even number in the hailstone sequence referenced in the Collatz conjecture. - Alex T. Flood (whiteangelsgrace(AT)gmail.com), Sep 22 2006
Let F(n) be the n-th Fermat number (A000215). Then F(a(r-1)) divides F(n)+2^k for r = k mod 2^n and r != 1. - T. D. Noe, Jul 12 2007
The following relation holds: 2^A007814(n)*(2*A025480(n-1)+1) = A001477(n) = n. (See functions hd, tl and cons in [Paul Tarau 2009].)
a(n) is the number of 0's at the end of n when n is written in base 2.
a(n+1) is the number of 1's at the end of n when n is written in base 2. - M. F. Hasler, Aug 25 2012
Shows which bit to flip when creating the binary reflected Gray code (bits are numbered from the right, offset is 0). That is, A003188(n) XOR A003188(n+1) == 2^A007814(n). - Russ Cox, Dec 04 2010
The sequence is squarefree (in the sense of not containing any subsequence of the form XX) [Allouche and Shallit]. Of course it contains individual terms that are squares (such as 4). - Comment expanded by N. J. A. Sloane, Jan 28 2019
a(n) is the number of zero coefficients in the n-th Stern polynomial, A125184. - T. D. Noe, Mar 01 2011
Lemma: For n < m with r = a(n) = a(m) there exists n < k < m with a(k) > r. Proof: We have n=b2^r and m=c2^r with b < c both odd; choose an even i between them; now a(i2^r) > r and n < i2^r < m. QED. Corollary: Every finite run of consecutive integers has a unique maximum 2-adic valuation. - Jason Kimberley, Sep 09 2011
a(n-2) is the 2-adic valuation of A000166(n) for n >= 2. - Joerg Arndt, Sep 06 2014
a(n) = number of 1's in the partition having Heinz number n. We define the Heinz number of a partition p = [p_1, p_2, ..., p_r] as Product_{j=1..r} p_j-th prime (concept used by Alois P. Heinz in A215366 as an "encoding" of a partition). For example, for the partition [1, 1, 2, 4, 10] we get 2*2*3*7*29 = 2436. Example: a(24)=3; indeed, the partition having Heinz number 24 = 2*2*2*3 is [1,1,1,2]. - Emeric Deutsch, Jun 04 2015
a(n+1) is the difference between the two largest parts in the integer partition having viabin number n (0 is assumed to be a part). Example: a(20) = 2. Indeed, we have 19 = 10011_2, leading to the Ferrers board of the partition [3,1,1]. For the definition of viabin number see the comment in A290253. - Emeric Deutsch, Aug 24 2017
Apart from being squarefree, as noted above, the sequence has the property that every consecutive subsequence contains at least one number an odd number of times. - Jon Richfield, Dec 20 2018
a(n+1) is the 2-adic valuation of Sum_{e=0..n} u^e = (1 + u + u^2 + ... + u^n), for any u of the form 4k+1 (A016813). - Antti Karttunen, Aug 15 2020
{a(n)} represents the "first black hat" strategy for the game of countably infinitely many hats, with a probability of success of 1/3; cf. the Numberphile link below. - Frederic Ruget, Jun 14 2021
a(n) is the least nonnegative integer k for which there does not exist i+j=n and a(i)=a(j)=k (cf. A322523). - Rémy Sigrist and Jianing Song, Aug 23 2022

Examples

			2^3 divides 24, so a(24)=3.
From _Omar E. Pol_, Jun 12 2009: (Start)
Triangle begins:
  0;
  1,0;
  2,0,1,0;
  3,0,1,0,2,0,1,0;
  4,0,1,0,2,0,1,0,3,0,1,0,2,0,1,0;
  5,0,1,0,2,0,1,0,3,0,1,0,2,0,1,0,4,0,1,0,2,0,1,0,3,0,1,0,2,0,1,0;
  6,0,1,0,2,0,1,0,3,0,1,0,2,0,1,0,4,0,1,0,2,0,1,0,3,0,1,0,2,0,1,0,5,0,1,0,2,...
(End)
		

References

  • J.-P. Allouche and J. Shallit, Automatic Sequences, Cambridge Univ. Press, 2003, p. 27.
  • K. Atanassov, On the 37th and the 38th Smarandache Problems, Notes on Number Theory and Discrete Mathematics, Sophia, Bulgaria, Vol. 5 (1999), No. 2, 83-85.
  • Michel Rigo, Formal Languages, Automata and Numeration Systems, 2 vols., Wiley, 2014. Mentions this sequence - see "List of Sequences" in Vol. 2.

Crossrefs

Cf. A011371 (partial sums), A094267 (first differences), A001511 (bisection), A346070 (mod 4).
Bisection of A050605 and |A088705|. Pairwise sums are A050603 and A136480. Difference of A285406 and A281264.
This is Guy Steele's sequence GS(1, 4) (see A135416). Cf. A053398(1,n). Column/row 1 of table A050602.
Cf. A007949 (3-adic), A235127 (4-adic), A112765 (5-adic), A122841 (6-adic), A214411 (7-adic), A244413 (8-adic), A122840 (10-adic).
Cf. A086463 (Dgf at s=2).

Programs

  • Haskell
    a007814 n = if m == 0 then 1 + a007814 n' else 0
                where (n', m) = divMod n 2
    -- Reinhard Zumkeller, Jul 05 2013, May 14 2011, Apr 08 2011
    
  • Haskell
    a007814 n | odd n = 0 | otherwise = 1 + a007814 (n `div` 2)
    --  Walt Rorie-Baety, Mar 22 2013
    
  • Magma
    [Valuation(n, 2): n in [1..120]]; // Bruno Berselli, Aug 05 2013
    
  • Maple
    ord := proc(n) local i,j; if n=0 then return 0; fi; i:=0; j:=n; while j mod 2 <> 1 do i:=i+1; j:=j/2; od: i; end proc: seq(ord(n), n=1..111);
    A007814 := n -> padic[ordp](n,2): seq(A007814(n), n=1..111); # Peter Luschny, Nov 26 2010
  • Mathematica
    Table[IntegerExponent[n, 2], {n, 64}] (* Eric W. Weisstein *)
    IntegerExponent[Range[64], 2] (* Eric W. Weisstein, Feb 01 2024 *)
    p=2; Array[ If[ Mod[ #, p ]==0, Select[ FactorInteger[ # ], Function[ q, q[ [ 1 ] ]==p ], 1 ][ [ 1, 2 ] ], 0 ]&, 96 ]
    DigitCount[BitXor[x, x - 1], 2, 1] - 1; a different version based on the same concept: Floor[Log[2, BitXor[x, x - 1]]] (* Jaume Simon Gispert (jaume(AT)nuem.com), Aug 29 2004 *)
    Nest[Join[ #, ReplacePart[ #, Length[ # ] -> Last[ # ] + 1]] &, {0, 1}, 5] (* N. J. Gunther, May 23 2009 *)
    Nest[ Flatten[# /. a_Integer -> {0, a + 1}] &, {0}, 7] (* Robert G. Wilson v, Jan 17 2011 *)
  • PARI
    A007814(n)=valuation(n,2);
    
  • Python
    import math
    def a(n): return int(math.log(n - (n & n - 1), 2)) # Indranil Ghosh, Apr 18 2017
    
  • Python
    def A007814(n): return (~n & n-1).bit_length() # Chai Wah Wu, Jul 01 2022
    
  • R
    sapply(1:100,function(x) sum(gmp::factorize(x)==2)) # Christian N. K. Anderson, Jun 20 2013
    
  • Scheme
    (define (A007814 n) (let loop ((n n) (e 0)) (if (odd? n) e (loop (/ n 2) (+ 1 e))))) ;; Antti Karttunen, Oct 06 2017

Formula

a(n) = A001511(n) - 1.
a(2*n) = A050603(2*n) = A001511(n).
a(n) = A091090(n-1) + A036987(n-1) - 1.
a(n) = 0 if n is odd, otherwise 1 + a(n/2). - Reinhard Zumkeller, Aug 11 2001
Sum_{k=1..n} a(k) = n - A000120(n). - Benoit Cloitre, Oct 19 2002
G.f.: A(x) = Sum_{k>=1} x^(2^k)/(1-x^(2^k)). - Ralf Stephan, Apr 10 2002
G.f. A(x) satisfies A(x) = A(x^2) + x^2/(1-x^2). A(x) = B(x^2) = B(x) - x/(1-x), where B(x) is the g.f. for A001151. - Franklin T. Adams-Watters, Feb 09 2006
Totally additive with a(p) = 1 if p = 2, 0 otherwise.
Dirichlet g.f.: zeta(s)/(2^s-1). - Ralf Stephan, Jun 17 2007
Define 0 <= k <= 2^n - 1; binary: k = b(0) + 2*b(1) + 4*b(2) + ... + 2^(n-1)*b(n-1); where b(x) are 0 or 1 for 0 <= x <= n - 1; define c(x) = 1 - b(x) for 0 <= x <= n - 1; Then: a(k) = c(0) + c(0)*c(1) + c(0)*c(1)*c(2) + ... + c(0)*c(1)...c(n-1); a(k+1) = b(0) + b(0)*b(1) + b(0)*b(1)*b(2) + ... + b(0)*b(1)...b(n-1). - Arie Werksma (werksma(AT)tiscali.nl), May 10 2008
a(n) = floor(A002487(n - 1) / A002487(n)). - Reikku Kulon, Oct 05 2008
Sum_{k=1..n} (-1)^A000120(n-k)*a(k) = (-1)^(A000120(n)-1)*(A000120(n) - A000035(n)). - Vladimir Shevelev, Mar 17 2009
a(A001147(n) + A057077(n-1)) = a(2*n). - Vladimir Shevelev, Mar 21 2009
For n>=1, a(A004760(n+1)) = a(n). - Vladimir Shevelev, Apr 15 2009
2^(a(n)) = A006519(n). - Philippe Deléham, Apr 22 2009
a(n) = A063787(n) - A000120(n). - Gary W. Adamson, Jun 04 2009
a(C(n,k)) = A000120(k) + A000120(n-k) - A000120(n). - Vladimir Shevelev, Jul 19 2009
a(n!) = n - A000120(n). - Vladimir Shevelev, Jul 20 2009
v_{2}(n) = Sum_{r>=1} (r / 2^(r+1)) Sum_{k=0..2^(r+1)-1} e^(2(k*Pi*i(n+2^r))/(2^(r+1))). - A. Neves, Sep 28 2010, corrected Oct 04 2010
a(n) mod 2 = A096268(n-1). - Robert G. Wilson v, Jan 18 2012
a(A005408(n)) = 1; a(A016825(n)) = 3; A017113(a(n)) = 5; A051062(a(n)) = 7; a(n) = (A037227(n)-1)/2. - Reinhard Zumkeller, Jun 30 2012
a((2*n-1)*2^p) = p, p >= 0 and n >= 1. - Johannes W. Meijer, Feb 04 2013
a(n) = A067255(n,1). - Reinhard Zumkeller, Jun 11 2013
a(n) = log_2(n - (n AND n-1)). - Gary Detlefs, Jun 13 2014
a(n) = 1 + A000120(n-1) - A000120(n), where A000120 is the Hamming weight function. - Stanislav Sykora, Jul 14 2014
A053398(n,k) = a(A003986(n-1,k-1)+1); a(n) = A053398(n,1) = A053398(n,n) = A053398(2*n-1,n) = Min_{k=1..n} A053398(n,k). - Reinhard Zumkeller, Aug 04 2014
a((2*x-1)*2^n) = a((2*y-1)*2^n) for positive n, x and y. - Juri-Stepan Gerasimov, Aug 04 2016
a(n) = A285406(n) - A281264(n). - Ralf Steiner, Apr 18 2017
a(n) = A000005(n)/(A000005(2*n) - A000005(n)) - 1. - conjectured by Velin Yanev, Jun 30 2017, proved by Nicholas Stearns, Sep 11 2017
Equivalently to above formula, a(n) = A183063(n) / A001227(n), i.e., a(n) is the number of even divisors of n divided by number of odd divisors of n. - Franklin T. Adams-Watters, Oct 31 2018
a(n)*(n mod 4) = 2*floor(((n+1) mod 4)/3). - Gary Detlefs, Feb 16 2019
Asymptotic mean: lim_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 1. - Amiram Eldar, Jul 11 2020
a(n) = 2*Sum_{j=1..floor(log_2(n))} frac(binomial(n, 2^j)*2^(j-1)/n). - Dario T. de Castro, Jul 08 2022
a(n) = A070939(n) - A070939(A030101(n)). - Andrew T. Porter, Dec 16 2022
a(n) = floor((gcd(n, 2^n)^(n+1) mod (2^(n+1)-1)^2)/(2^(n+1)-1)) (see Lemma 3.4 from Mazzanti's 2002 article). - Lorenzo Sauras Altuzarra, Mar 10 2024
a(n) = 1 - A088705(n). - Chai Wah Wu, Sep 18 2024

Extensions

Formula index adapted to the offset of A025480 by R. J. Mathar, Jul 20 2010
Edited by Ralf Stephan, Feb 08 2014

A033428 a(n) = 3*n^2.

Original entry on oeis.org

0, 3, 12, 27, 48, 75, 108, 147, 192, 243, 300, 363, 432, 507, 588, 675, 768, 867, 972, 1083, 1200, 1323, 1452, 1587, 1728, 1875, 2028, 2187, 2352, 2523, 2700, 2883, 3072, 3267, 3468, 3675, 3888, 4107, 4332, 4563, 4800, 5043, 5292, 5547, 5808, 6075, 6348
Offset: 0

Views

Author

Keywords

Comments

The number of edges of a complete tripartite graph of order 3n, K_n,n,n. - Roberto E. Martinez II, Oct 18 2001
From Floor van Lamoen, Jul 21 2001: (Start)
Write 1,2,3,4,... in a hexagonal spiral around 0; then a(n) is the sequence found by reading the line from 0 in the direction 0,3,.... The spiral begins:
.
33--32--31--30
/ \
34 16--15--14 29
/ / \ \
35 17 5---4 13 28
/ / / \ \ \
36 18 6 0---3--12--27--48-->
/ / / / / / / /
37 19 7 1---2 11 26 47
\ \ \ / / /
38 20 8---9--10 25 46
\ \ / /
39 21--22--23--24 45
\ /
40--41--42--43--44
(End)
Number of edges of the complete bipartite graph of order 4n, K_n,3n. - Roberto E. Martinez II, Jan 07 2002
Also the number of partitions of 6n + 3 into at most 3 parts. - R. K. Guy, Oct 23 2003
Also the number of partitions of 6n into exactly 3 parts. - Colin Barker, Mar 23 2015
Numbers n such that the imaginary quadratic field Q[sqrt(-n)] has six units. - Marc LeBrun, Apr 12 2006
The denominators of Hoehn's sequence (recalled by G. L. Honaker, Jr.) and the numerators of that sequence reversed. The sequence is 1/3, (1+3)/(5+7), (1+3+5)/(7+9+11), (1+3+5+7)/(9+11+13+15), ...; reduced to 1/3, 4/12, 9/27, 16/48, ... . For the reversal, the reduction is 3/1, 12/4, 27/9, 48/16, ... . - Enoch Haga, Oct 05 2007
Right edge of tables in A200737 and A200741: A200737(n, A000292(n)) = A200741(n, A100440(n)) = a(n). - Reinhard Zumkeller, Nov 21 2011
The Wiener index of the crown graph G(n) (n>=3). The crown graph G(n) is the graph with vertex set {x(1), x(2), ..., x(n), y(1), y(2), ..., y(n)} and edge set {(x(i), y(j)): 1<=i, j<=n, i/=j} (= the complete bipartite graph K(n,n) with horizontal edges removed). Example: a(3)=27 because G(3) is the cycle C(6) and 6*1 + 6*2 + 3*3 = 27. The Hosoya-Wiener polynomial of G(n) is n(n-1)(t+t^2)+nt^3. - Emeric Deutsch, Aug 29 2013
From Michel Lagneau, May 04 2015: (Start)
Integer area A of equilateral triangles whose side lengths are in the commutative ring Z[3^(1/4)] = {a + b*3^(1/4) + c*3^(1/2) + d*3^(3/4), a,b,c and d in Z}.
The area of an equilateral triangle of side length k is given by A = k^2*sqrt(3)/4. In the ring Z[3^(1/4)], if k = q*3^(1/4), then A = 3*q^2/4 is an integer if q is even. Example: 27 is in the sequence because the area of the triangle (6*3^(1/4), 6*3^(1/4), 6*3^(1/4)) is 27. (End)
a(n) is 2*sqrt(3) times the area of a 30-60-90 triangle with short side n. Also, 3 times the area of an n X n square. - Wesley Ivan Hurt, Apr 06 2016
Consider the hexagonal tiling of the plane. Extract any four hexagons adjacent by edge. This can be done in three ways. Fold the four hexagons so that all opposite faces occupy parallel planes. For all parallel projections of the resulting object, at least two correspond to area a(n) for side length of n of the original hexagons. - Torlach Rush, Aug 17 2022
The sequence terms are the exponents in the expansion of Product_{n >= 1} (1 - q^(3*n))/(1 + q^(3*n)) = ( Sum_{n in Z} q^(n*(3*n+1)/2) ) / ( Product_{n >= 1} 1 + q^n ) = 1 - 2*q^3 + 2*q^12 - 2*q^27 + 2*q^48 - 2*q^75 + - .... - Peter Bala, Dec 30 2024

Examples

			From _Ilya Gutkovskiy_, Apr 13 2016: (Start)
Illustration of initial terms:
.                                              o
.                                             o o
.                                            o   o
.                          o                o  o  o
.                         o o              o  o o  o
.                        o   o            o  o   o  o
.           o           o  o  o          o  o  o  o  o
.          o o         o  o o  o        o  o  o o  o  o
.         o   o       o  o   o  o      o  o  o   o  o  o
.  o     o  o  o     o  o  o  o  o    o  o  o  o  o  o  o
. o o   o  o o  o   o  o  o o  o  o  o  o  o  o o  o  o  o
. n=1      n=2            n=3                 n=4
(End)
		

Crossrefs

Programs

  • Haskell
    a033428 = (* 3) . (^ 2)
    a033428_list = 0 : 3 : 12 : zipWith (+) a033428_list
       (map (* 3) $ tail $ zipWith (-) (tail a033428_list) a033428_list)
    -- Reinhard Zumkeller, Jul 11 2013
    
  • Magma
    [3*n^2: n in [0..50]]; // Vincenzo Librandi, May 18 2015
    
  • Maple
    seq(3*n^2, n=0..46); # Nathaniel Johnston, Jun 26 2011
  • Mathematica
    3 Range[0, 50]^2
    LinearRecurrence[{3, -3, 1}, {0, 3, 12}, 50] (* Harvey P. Dale, Feb 16 2013 *)
  • Maxima
    makelist(3*n^2,n,0,30); /* Martin Ettl, Nov 12 2012 */
    
  • PARI
    a(n)=3*n^2
    
  • Python
    def a(n): return 3 * (n**2) # Torlach Rush, Aug 25 2022

Formula

a(n) = 3*a(n-1)-3*a(n-2)+a(n-3) for n>2.
G.f.: 3*x*(1+x)/(1-x)^3. - R. J. Mathar, Sep 09 2008
Main diagonal of triangle in A132111: a(n) = A132111(n,n). - Reinhard Zumkeller, Aug 10 2007
A214295(a(n)) = -1. - Reinhard Zumkeller, Jul 12 2012
a(n) = A215631(n,n) for n > 0. - Reinhard Zumkeller, Nov 11 2012
a(n) = A174709(6n+2). - Philippe Deléham, Mar 26 2013
a(n) = a(n-1) + 6*n - 3, with a(0)=0. - Jean-Bernard François, Oct 04 2013
E.g.f.: 3*x*(1 + x)*exp(x). - Ilya Gutkovskiy, Apr 13 2016
a(n) = t(3*n) - 3*t(n), where t(i) = i*(i+k)/2 for any k. Special case (k=1): A000217(3*n) - 3*A000217(n). - Bruno Berselli, Aug 31 2017
a(n) = A000326(n) + A005449(n). - Bruce J. Nicholson, Jan 10 2020
From Amiram Eldar, Jul 03 2020: (Start)
Sum_{n>=1} 1/a(n) = Pi^2/18 (A086463).
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi^2/36. (End)
From Amiram Eldar, Feb 03 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = sqrt(3)*sinh(Pi/sqrt(3))/Pi.
Product_{n>=1} (1 - 1/a(n)) = sqrt(3)*sin(Pi/sqrt(3))/Pi. (End)
a(n) = A003215(n) - A016777(n). - Leo Tavares, Apr 29 2023

Extensions

Better description from N. J. A. Sloane, May 15 1998

A002390 Decimal expansion of natural logarithm of golden ratio.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

The Baxa article proves that every gamma >= this constant is the Lévy constant of a transcendental number. - Michel Marcus, Apr 09 2016
The entropy of the golden mean shift. See Capobianco link. - Michel Marcus, Jan 19 2019
Also the limiting value of the area of the function y = 1/x bounded by the abscissa of consecutive F(n) points (where F(n)=A000045(n) are the Fibonacci numbers and n > 0). - Burak Muslu, May 09 2021

Examples

			0.481211825059603447497758913424368423135184334385660519661...
		

References

  • George Boros and Victor H. Moll, Irresistible integrals, Cambridge University Press (2006), p. 236.
  • W. E. Mansell, Tables of Natural and Common Logarithms. Royal Society Mathematical Tables, Vol. 8, Cambridge Univ. Press, 1964, p. XVIII.
  • B. Muslu, Sayılar ve Bağlantılar 2, Luna, 2021, pages 31-38.
  • 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).

Crossrefs

Programs

Formula

Also equals arcsinh(1/2).
Equals sqrt(5)* A086466 /2. - Seiichi Kirikami, Aug 20 2011
Equals sqrt(5)*(5* A086465 -1)/4. - Jean-François Alcover, Apr 29 2013
Also equals (125*C - 55) / (24*sqrt(5)), where C = Sum_{k>=1} (-1)^(k+1)*1/Cat(k), where Cat(k) = (2k)!/k!/(k+1)! = A000108(k) - k-th Catalan number. See Sep 01 2006 comment at ref. Mathematics in Russian. - Alexander Adamchuk, Dec 27 2013
Equals sqrt(5)/4 * Sum_{n>=0} (-1)^n/((2n+1)*C(2*n,n)) = sqrt(5) *A344041 /4. - Alexander Adamchuk, Dec 27 2013
Equals sqrt((Pi^2/6 - W)/3), where W = Sum_{n>=0} (-1)^n/((2n+1)^2*C(2*n,n)) = A145436, attributed by Alexander Adamchuk to Ramanujan. See Sep 01 2006 comment at ref. Mathematics in Russian. - Alexander Adamchuk, Dec 27 2013
Equals lim_{j->infinity} Sum_{k=F(j)..F(j+1)-1} (1/k), where F = A000045, the Fibonacci sequence. Convergence is slow. For example: Sum_{k=21..33} (1/k) = 0.4910585.... - Richard R. Forberg, Aug 15 2014
Equals Sum_{k>=1} cos(Pi*k/5)/k. - Amiram Eldar, Aug 12 2020
Equals real solution to exp(x)+exp(2*x) = exp(3*x). - Alois P. Heinz, Jul 14 2022
Equals arccoth(sqrt(5)). - Amiram Eldar, Feb 09 2024
Sum_{n >= 1} 1/(n*P(n, sqrt(5))*P(n-1, sqrt(5))), where P(n, x) denotes the n-th Legendre polynomial. The first ten terms of the series gives the approximation log((1 + sqrt(5))/2) = 0.481211825059(39..), correct to 12 decimal places. - Peter Bala, Mar 16 2024
Equals Sum_{n>=0} ((-1)^(n)*binomial(2*n, n))/(2^(4*n + 1)*(2*n + 1)). - Antonio Graciá Llorente, Nov 13 2024

A100044 Decimal expansion of Pi^2/9.

Original entry on oeis.org

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

Views

Author

Eric W. Weisstein, Oct 31 2004

Keywords

Comments

The Dirichlet L-series for the principal character mod 6 (which is A120325 shifted left) evaluated at 2. - R. J. Mathar, Jul 20 2012
Equals the asymptotic mean of the abundancy index of the numbers coprime to 6 (A007310). - Amiram Eldar, May 12 2023

Examples

			1.096622711232150957648276777764...
		

References

  • F. Aubonnet, D. Guinin, and B.Joppin, Précis de Mathématiques, Analyse 2, Classes Préparatoires, Premier Cycle Universitaire, Bréal, 1990, Exercice 908, pages 82 and 91-92.
  • L. B. W. Jolley, Summation of Series, Dover, 1961.

Crossrefs

Programs

  • Mathematica
    RealDigits[Pi^2/9, 10, 110][[1]] (* G. C. Greubel, Feb 17 2017 *)
  • PARI
    default(realprecision, 110); Pi^2/9 \\ G. C. Greubel, Feb 17 2017
    
  • Sage
    numerical_approx(pi^2/9, digits=120) # G. C. Greubel, Jun 02 2021

Formula

Equals 1 + (1/2)*(1/3)*(1/2) + (1/3)*(1*2)/(3*5)*(1/2)^2 + (1/4) *(1*2*3)/(3*5*7)*(1/2)^3 + .... [Jolley eq 277]
Equals 1/1^2 + 1/5^2 + 1/7^2 + 1/11^2 + 1/13^2 + 1/17^2 + .... - R. J. Mathar, Jul 20 2012
Equals 2*Sum_{n>=1} 1/(6*n*(3*n + (-1)^n - 3) - 3*(-1)^n + 5) = 2*Sum_{n>=1} 1/(2*A104777(n)). - Alexander R. Povolotsky, May 18 2014
Equals A019670^2. - Michel Marcus, May 19 2014
Equals 2*A086463 = 2*Sum_{n>=1} 1/A091999(n)^2, equivalent to the formula of 2012 above. - Alexander R. Povolotsky, May 20 2014
Equals 3F2(1,1,1; 3/2,2 ; 1/4), following from Clausen's formula of J. Reine Angew. Math 3 (1828) for squares of 2F1() as noted in A019670. - R. J. Mathar, Oct 16 2015
Equals Product_{n >= 3} prime(n)^2 / (prime(n)^2 - 1), Euler's prime product, excluding first two primes. - Fred Daniel Kline, Jun 09 2016
Equals Integral_{x=0..oo} log(x)/(x^6 - 1) dx. - Amiram Eldar, Aug 12 2020
Equals Sum_{k>=1} A000120(k) * (2*k+1)/(k^2*(k+1)^2) (Allouche and Shallit, 1990). - Amiram Eldar, Jun 01 2021
Equals Integral_{x=0..1} log(1+x+x^2)/x dx (Aubonnet). - Bernard Schott, Feb 04 2022
Equals Sum_{k>=1} A008833(k)/k^4. - Amiram Eldar, Jan 25 2024
Continued fraction expansion: 1/(1 - 1/(13 - 48/(34 - 270/(65 - ... - 2*(2*n-1)*n^3/((5*n^2+6*n+2) - ... ))))). See A130549. - Peter Bala, Feb 16 2024
Equals Sum_{k >= 0} 1/((k + 1)*(2*k + 1)*binomial(2*k, k)). See Catalan, Section 21, equation 30. - Peter Bala, Aug 14 2024

A182448 Decimal expansion of Pi^2/15.

Original entry on oeis.org

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

Views

Author

Mats Granvik, Apr 29 2012

Keywords

Examples

			0.65797362673929...
		

Crossrefs

Programs

  • Mathematica
    RealDigits[N[Sum[1/(n + 0)^2 - 1/(n + 1)^2 + 1/(n + 2)^2 - 1/(n + 3)^2 - 4/(n + 4)^2 - 1/(n + 5)^2 + 1/(n + 6)^2 - 1/(n + 7)^2 + 1/(n + 8)^2 + 4/(n + 9)^2, {n, 1, Infinity, 10}], 90]][[1]]
    RealDigits[N[Sum[LiouvilleLambda[n]/n^2, {n, 1, Infinity}], 90]][[1]]
    RealDigits[Pi^2/15,10,120][[1]] (* Harvey P. Dale, May 28 2017 *)
  • PARI
    Pi^2/15 \\ Michel Marcus, Oct 21 2014

Formula

See Mathematica code.
Equals Gamma(4)*zeta(4)/Pi^2 = zeta(4)/zeta(2) = A013662/A013661 = Product_{p prime} (p^2/(p^2+1)). - Stanislav Sykora, Oct 21 2014
Equals (1/10) * Sum_{n >= 0} (-1)^n*( 1/(n + 1/3)^2 - 1/(n + 2/3)^2 ). - Peter Bala, Oct 31 2019
Equals Sum_{k>=1} A008836(k)/k^2. - Amiram Eldar, Jun 23 2020
Equals (1/10) * Sum_{k>=1} (5*t(k-1) + 3*t(k))/k^2, where t(k) = A010060(k) (Tóth, 2022). - Amiram Eldar, Feb 04 2024
Equals 3/5 + (1/5) * Sum_{n>=1} 1/(n^2*(n+1)^2). - Davide Rotondo, May 28 2025
Equals 1/A082020 = A164102/30 = A195055/5. - Hugo Pfoertner, May 28 2025

Extensions

Offset corrected and more terms added by Rick L. Shepherd, Jan 08 2014

A197863 Smallest powerful number that is a multiple of n.

Original entry on oeis.org

1, 4, 9, 4, 25, 36, 49, 8, 9, 100, 121, 36, 169, 196, 225, 16, 289, 36, 361, 100, 441, 484, 529, 72, 25, 676, 27, 196, 841, 900, 961, 32, 1089, 1156, 1225, 36, 1369, 1444, 1521, 200, 1681, 1764, 1849, 484, 225, 2116, 2209, 144, 49, 100, 2601, 676, 2809, 108, 3025
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • Haskell
    a197863 n = product $
       zipWith (^) (a027748_row n) (map (max 2) $ a124010_row n)
    -- Reinhard Zumkeller, Jan 06 2012
  • Mathematica
    With[{pwrnos=Join[{1},Select[Range[5000],Min[Transpose[ FactorInteger[#]] [[2]]]>1&]]},Flatten[Table[Select[pwrnos,Divisible[#,n]&,1],{n,60}]]] (* Harvey P. Dale, Aug 14 2012 *)
    f[p_, e_] := p^Max[e, 2]; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Jul 09 2022 *)
  • PARI
    a(n)=local(fm=factor(n));prod(k=1,matsize(fm)[1],fm[k,1]^max(fm[k,2],2))
    

Formula

Multiplicative with a(p^e) = p^max(e,2).
Sum_{n>=1} 1/a(n) = Product_{p prime} (1 + (2*p-1)/(p^2*(p-1))) = 2.71098009471568319328... . - Amiram Eldar, Jul 29 2022
Sum_{k=1..n} a(k) ~ c * n^3, where c = (Pi^2/18) * Product_{p prime} (1 - 2/p^2 + 2/p^4 - 1/p^5) = 0.2165355664... . - Amiram Eldar, Nov 19 2022
a(n) = n * A055231(n). - Amiram Eldar, Sep 01 2023

A353908 Decimal expansion of Pi^2/36.

Original entry on oeis.org

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

Views

Author

Omar E. Pol, May 10 2022

Keywords

Comments

Ratio between the volume of the stepped pyramid with an infinite number of levels described in A245092 and that of the circumscribed cube (see the first formula).
See also Vaclav Kotesovec's formula (2016) in A175254.
Volume shared by a sphere inscribed in a cube of volume Pi and one of the six pyramids inscribed in the cube. - Omar E. Pol, Sep 01 2024

Examples

			0.2741556778080377394120691944410041982031583168677997396225930382283345784...
		

Crossrefs

Programs

  • Maple
    evalf(Pi^2/36, 121);  # Alois P. Heinz, May 11 2022
  • Mathematica
    RealDigits[Pi^2/36, 10, 100][[1]] (* Amiram Eldar, May 11 2022 *)
  • PARI
    Pi^2/36
    
  • PARI
    zeta(2)/6

Formula

Equals lim_{n->oo} A175254(n)/n^3.
Equals A002388/36.
Equals A102753/18.
Equals A195055/12.
Equals A091476/9.
Equals A013661/6.
Equals A100044/4.
Equals A072691/3.
Equals A086463/2.
Equals A086729*2.
Equals A019673^2.
Equals Re(dilog((1+sqrt(3)*i)/2)). - Mohammed Yaseen, Jul 03 2024

A051553 Decimal expansion of sin (Pi degrees). (Of course sin (Pi radians) = 0.)

Original entry on oeis.org

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

Views

Author

Hsu, Po-Wei (Benny) (arsene_lupin(AT)intekom.co.za)

Keywords

Examples

			0.054803665148789530887748713539833254... [corrected by _Georg Fischer_, Jul 29 2021]
		

Crossrefs

Cf. A086463 (Pi^2/18).

Programs

A086464 Decimal expansion of 17/36*zeta(4).

Original entry on oeis.org

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

Views

Author

Eric W. Weisstein, Jul 21 2003

Keywords

Examples

			0.51109708258581525710477952336666262075474350507273...
		

Crossrefs

Programs

  • Mathematica
    RealDigits[17*Zeta[4]/36, 10, 120][[1]] (* Amiram Eldar, May 25 2023 *)
  • PARI
    zeta(4)*17/36 \\ Michel Marcus, Jul 31 2015

Formula

Equals Sum_{n>=1} 1/(n^4 * binomial(2*n,n)).

A091259 Numerator of sigma_3(n)/sigma(n).

Original entry on oeis.org

1, 3, 7, 73, 21, 21, 43, 39, 757, 63, 111, 73, 157, 129, 147, 151, 273, 2271, 343, 219, 301, 333, 507, 273, 15751, 471, 511, 3139, 813, 441, 931, 4161, 777, 819, 903, 55261, 1333, 1029, 1099, 819, 1641, 903, 1807, 8103, 15897, 1521, 2163, 1057, 39331, 47253
Offset: 1

Views

Author

Labos Elemer, Feb 12 2004

Keywords

Crossrefs

Cf. A032766.

Programs

  • Magma
    [Numerator(DivisorSigma(3,n)/DivisorSigma(1,n)): n in [1..50]]; // Vincenzo Librandi, Jan 26 2018
  • Maple
    seq(numer(numtheory:-sigma[3](n)/numtheory:-sigma(n)),n=1..100); # Robert Israel, Jan 25 2018
  • Mathematica
    Array[Numerator[DivisorSigma[3,#]/DivisorSigma[1,#]]&,50] (* Harvey P. Dale, Feb 29 2016 *)
  • PARI
    a(n) = numerator(sigma(n, 3)/sigma(n)); \\ Michel Marcus, Jan 26 2018
    

Formula

a(p) = A002061(p), for prime p. - Robert Israel, Jan 25 2018
Sum_{k=1..n} a(k)/A091258(k) ~ c * n^3, where c = (Pi^2/18)*zeta(3)^2 * Product_{p prime} (1 - 2/p^2 - 1/p^3 + 5/p^5 - 3/p^6) = 0.2382648075... . - Amiram Eldar, Nov 21 2022
Conjecture: a(n) (mod 3) = A353816(n). - Michel Marcus, Aug 11 2024
Showing 1-10 of 29 results. Next