Search: a013679 -id:a013679
|
|
A002117
|
|
Apéry's number or Apéry's constant zeta(3). Decimal expansion of zeta(3) = Sum_{m >= 1} 1/m^3.
(Formerly M0020)
|
|
+10
431
|
|
|
1, 2, 0, 2, 0, 5, 6, 9, 0, 3, 1, 5, 9, 5, 9, 4, 2, 8, 5, 3, 9, 9, 7, 3, 8, 1, 6, 1, 5, 1, 1, 4, 4, 9, 9, 9, 0, 7, 6, 4, 9, 8, 6, 2, 9, 2, 3, 4, 0, 4, 9, 8, 8, 8, 1, 7, 9, 2, 2, 7, 1, 5, 5, 5, 3, 4, 1, 8, 3, 8, 2, 0, 5, 7, 8, 6, 3, 1, 3, 0, 9, 0, 1, 8, 6, 4, 5, 5, 8, 7, 3, 6, 0, 9, 3, 3, 5, 2, 5, 8, 1, 4, 6, 1, 9, 9, 1, 5
(list;
constant;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Sometimes called Apéry's constant.
"A natural question is whether Zeta(3) is a rational multiple of Pi^3. This is not known, though in 1978 R. Apéry succeeded in proving that Zeta(3) is irrational. In Chapter 8 we pointed out that the probability that two random integers are relatively prime is 6/Pi^2, which is 1/Zeta(2). This generalizes to: The probability that k random integers are relatively prime is 1/Zeta(k) ... ." [Stan Wagon]
In 2001 Tanguy Rivoal showed that there are infinitely many odd (positive) integers at which zeta is irrational, including at least one value j in the range 5 <= j <= 21 (refined the same year by Zudilin to 5 <= j <= 11), at which zeta(j) is irrational. See the Rivoal link for further information and references.
The reciprocal of this constant is the probability that three integers chosen randomly using uniform distribution are relatively prime. - Joseph Biberstine (jrbibers(AT)indiana.edu), Apr 13 2005
Also the value of zeta(1,2), the double zeta-function of arguments 1 and 2. - R. J. Mathar, Oct 10 2011
Also the length of minimal spanning tree for large complete graph with uniform random edge lengths between 0 and 1, cf. link to John Baez's comment. - M. F. Hasler, Sep 26 2017
This number is the average value of sigma_2(n)/n^2 where sigma_2(n) is the sum of the squares of the divisors of n. - Dimitri Papadopoulos, Jan 07 2022
|
|
REFERENCES
|
S. R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, pp. 40-53.
A. Fletcher, J. C. P. Miller, L. Rosenhead, and L. J. Comrie, An Index of Mathematical Tables. Vols. 1 and 2, 2nd ed., Blackwell, Oxford and Addison-Wesley, Reading, MA, 1962, Vol. 1, p. 84.
R. William Gosper, Strip Mining in the Abandoned Orefields of Nineteenth Century Mathematics, Computers in Mathematics (Stanford CA, 1986); Lecture Notes in Pure and Appl. Math., Dekker, New York, 125 (1990), 261-284; MR 91h:11154.
Xavier Gourdon, Analyse, Les Maths en tête, Ellipses, 1994, Exemple 3, page 224.
Richard K. Guy, Unsolved Problems in Number Theory, 3rd Edition, Springer, 2004, Section F17, Series associated with the zeta-function, p. 391.
G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, Oxford University Press; 6 edition (2008), pp. 47, 268-269.
Paul Levrie, The Ubiquitous Apéry Number, Math. Intelligencer, Vol. 45, No. 2, 2023, pp. 118-119.
A. A. Markoff, Mémoire sur la transformation de séries peu convergentes en séries très convergentes, Mém. de l'Acad. Imp. Sci. de St. Pétersbourg, XXXVII, 1890.
Paul J. Nahin, In Pursuit of Zeta-3, Princeton University Press, 2021.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
Stan Wagon, Mathematica In Action, W. H. Freeman and Company, NY, 1991, page 354.
David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987, p. 33.
A. M. Yaglom and I. M. Yaglom, Challenging Mathematical Problems with Elementary Solutions, Dover (1987), Ex. 92-93.
|
|
LINKS
|
|
|
FORMULA
|
Lima gives an approximation to zeta(3) as (236*log(2)^3)/197 - 283/394*Pi*log(2)^2 + 11/394*Pi^2*log(2) + 209/394*log(sqrt(2) + 1)^3 - 5/197 + (93*Catalan*Pi)/197. - Jonathan Vos Post, Oct 14 2009 [Corrected by Wouter Meeussen, Apr 04 2010]
zeta(3) = 5/2*Integral_(x=0..2*log((1+sqrt(5))/2), x^2/(exp(x)-1)) + 10/3*(log((1+sqrt(5))/2))^3. - Seiichi Kirikami, Aug 12 2011
zeta(3) = -4/3*Integral_{x=0..1} log(x)/x*log(1+x) = Integral_{x=0..1} log(x)/x*log(1-x) = -4/7*Integral_{x=0..1} log(x)/x*log((1+x)/(1-x)) = 4*Integral_{x=0..1} 1/x*log(1+x)^2 = 1/2*Integral_{x=0..1} 1/x*log(1-x)^2 = -16/7*Integral_{x=0..Pi/2} x*log(2*cos(x)) = -4/Pi*Integral_{x=0..Pi/2} x^2*log(2*cos(x)). - Jean-François Alcover, Apr 02 2013, after R. J. Mathar
zeta(3) = (16/7)*Sum_{k even} (k^3 + k^5)/(k^2 - 1)^4.
zeta(3) - 1 = Sum_{k >= 1} 1/(k^3 + 4*k^7) = 1/(5 - 1^6/(21 - 2^6/(55 - 3^6/(119 - ... - (n - 1)^6/((2*n - 1)*(n^2 - n + 5) - ...))))) (continued fraction).
More generally, there is a sequence of polynomials P(n,x) (of degree 2*n) such that
zeta(3) - Sum_{k = 1..n} 1/k^3 = Sum_{k >= 1} 1/( k^3*P(n,k-1)*P(n,k) ) = 1/((2*n^2 + 2*n + 1) - 1^6/(3*(2*n^2 + 2*n + 3) - 2^6/(5*(2*n^2 + 2*n + 7) - 3^6/(7*(2*n^2 + 2*n + 13) - ...)))) (continued fraction). See A143003 and A143007 for details.
Series acceleration formulas:
zeta(3) = (5/2)*Sum_{n >= 1} (-1)^(n+1)/( n^3*binomial(2*n,n) )
= (5/2)*Sum_{n >= 1} P(n)/( (2*n(2*n - 1))^3*binomial(4*n,2*n) )
= (5/2)*Sum_{n >= 1} (-1)^(n+1)*Q(n)/( (3*n(3*n - 1)*(3*n - 2))^3*binomial(6*n,3*n) ), where P(n) = 24*n^3 + 4*n^2 - 6*n + 1 and Q(n) = 9477*n^6 - 11421*n^5 + 5265*n^4 - 1701*n^3 + 558*n^2 - 108*n + 8 (Bala, section 7). (End)
zeta(3) = Sum_{n >= 1} (A010052(n)/n^(3/2)) = Sum_{n >= 1} ( (floor(sqrt(n)) - floor(sqrt(n-1)))/n^(3/2) ). - Mikael Aaltonen, Feb 22 2015
zeta(3) = 4*(2*log(2) - 1 - 2*Sum_{k>=2} zeta(2*k+1)/2^(2*k+1)). - Jorge Coveiro, Jun 21 2020
zeta(3) = (4*zeta'''(1/2)*(zeta(1/2))^2-12*zeta(1/2)*zeta'(1/2)*zeta''(1/2)+8*(zeta'(1/2))^3-Pi^3*(zeta(1/2))^3)/(28*(zeta(1/2))^3). - Artur Jasinski, Jun 27 2020
zeta(3) = (5/4)*Li_3(1/f^2) + Pi^2*log(f)/6 - 5*log(f)^3/6,
zeta(3) = (8/7)*Li_3(1/2) + (2/21)*Pi^2 log(2) - (4/21) log(2)^3, where f is golden ratio (A001622) and Li_3 is the polylogarithm function, formulas published by John Landen in 1780, p. 118. (End)
zeta(3) = (1/2)*Integral_{x=0..oo} x^2/(e^x-1) dx (Gourdon). - Bernard Schott, Apr 28 2021
zeta(3) = 1 + Sum_{n >= 1} 1/(n^3*(4*n^4 + 1)) = 25/24 + (2!)^4*Sum_{n >= 1} 1/(n^3*(4*n^4 + 1)*(4*n^4 + 2^4)) = 28333/27000 + (3!)^4*Sum_{n >= 1} 1/(n^3*(4*n^4 + 1)*(4*n^4 + 2^4)*(4*n^4 + 3^4)). In general, for k >= 1, we have zeta(3) = r(k) + (k!)^4*Sum_{n >= 1} 1/(n^3*(4*n^4 + 1)*...*(4*n^4 + k^4)), where r(k) is rational.
zeta(3) = (6/7) + (64/7)*Sum_{n >= 1} n/(4*n^2 - 1)^3.
More generally, for k >= 0, it appears that zeta(3) = a(k) + b(k)*Sum_{n >= 1} n/( (4*n^2 - 1)*(4*n^2 - 9)*...*(4*n^2 - (2*k+1)^2) )^3, where a(k) and b(k) are rational.
zeta(3) = (10/7) - (128/7)*Sum_{n >= 1} n/(4*n^2 - 1)^4.
More generally, for k >= 0, it appears that zeta(3) = c(k) + d(k)*Sum_{n >= 1} n/( (4*n^2 - 1)*(4*n^2 - 9)*...*(4*n^2 - (2*k+1)^2) )^4, where c(k) and d(k) are rational. [added Nov 27 2023: for the values of a(k), b(k), c(k) and d(k) see the Bala 2023 link, Sections 8 and 9.]
zeta(3) = 2/3 + (2^13)/(3*7)*Sum_{n >= 1} n^3/(4*n^2 - 1)^6. (End)
zeta(3) = -Psi(2)(1/2)/14 (the second derivative of digamma function evaluated at 1/2). - Artur Jasinski, Mar 18 2022
zeta(3) = -(8*Pi^2/9) * Sum_{k>=0} zeta(2*k)/((2*k+1)*(2*k+3)*4^k) = (2*Pi^2/9) * (log(2) + 2 * Sum_{k>=0} zeta(2*k)/((2*k+3)*4^k)) (Scheufens, 2011, Glasser Math. Comp. 22 1968). - Amiram Eldar, May 28 2022
zeta(3) = Sum_{k>=1} (30*k-11) / (4*(2k-1)*k^3*(binomial(2k,k))^2) (Gosper, 1986 and Richard K. Guy reference). - Bernard Schott, Jul 20 2022
zeta(3) = (4/3)*Integral_{x >= 1} x*log(x)*(1 + log(x))*log(1 + 1/x^x) dx = (2/3)*Integral_{x >= 1} x^2*log(x)^2*(1 + log(x))/(1 + x^x) dx. - Peter Bala, Nov 27 2023
zeta_3(n) = 1/180*(-360*n^3*f(-3, n/4) + Pi^3*(n^4 + 20*n^2 + 16))/(n*(n^2 + 4)), where f(-3, n) = Sum_{k>=1} 1/(k^3*(exp(Pi*k/n) - 1)). Will give at least 1 digit of precision/term, example: zeta_3(5) = 1.202056944732.... - Simon Plouffe, Dec 21 2023
zeat(3) = 1 + (1/2)*Sum_{n >= 1} (2*n + 1)/(n^3*(n + 1)^3) = 5/4 - (1/4)*Sum_{n >= 1} (2*n + 1)/(n^4*(n + 1)^4) = 147/120 + (2/15)*Sum_{n >= 1} (2*n + 1)/(n^5*(n + 1)^5) - (64/15)*Sum_{n >= 1} (n + 1)/(n^5*(n + 2)^5) = 19/16 + (128/21)*Sum_{n >= 1} (n + 1)/(n^6*(n + 2)^6) - (1/21)*Sum_{n >= 1} (2*n + 1)/(n^6*(n + 1)^6). - Peter Bala, Apr 15 2024
|
|
EXAMPLE
|
1.2020569031595942853997...
|
|
MAPLE
|
# Calculates an approximation with n exact decimal places (small deviation
# in the last digits are possible). Goes back to ideas of A. A. Markoff 1890.
zeta3 := proc(n) local s, w, v, k; s := 0; w := -1; v := 4;
for k from 2 by 2 to 7*n/2 do
w := -w*v/k;
v := v + 8;
s := s + 1/(w*k^3);
od; 20*s; evalf(%, n) end:
|
|
MATHEMATICA
|
RealDigits[ N[ Zeta[3], 100] ] [ [1] ]
(* Second program (historical interest): *)
d[n_] := 34*n^3 + 51*n^2 + 27*n + 5; 6/Fold[Function[d[#2-1] - #2^6/#1], 5, Reverse[Range[100]]] // N[#, 108]& // RealDigits // First
|
|
PROG
|
(PARI) default(realprecision, 20080); x=zeta(3); for (n=1, 20000, d=floor(x); x=(x-d)*10; write("b002117.txt", n, " ", d)); \\ Harry J. Smith, Apr 19 2009
(Maxima) fpprec : 100$ ev(bfloat(zeta(3)))$ bfloat(%); /* Martin Ettl, Oct 21 2012 */
(Python)
from mpmath import mp, apery
mp.dps=109
print([int(z) for z in list(str(apery).replace('.', ''))[:-1]]) # Indranil Ghosh, Jul 08 2017
(Magma) L:=RiemannZeta(: Precision:=100); Evaluate(L, 3); // G. C. Greubel, Aug 21 2018
|
|
CROSSREFS
|
Cf. A013631, A013679, A013661, A013663, A013667, A013669, A013671, A013675, A013677, A059956 (6/Pi^2), A084225; A084226.
Cf. sums of inverses: A152623 (tetrahedral numbers), A175577 (octahedral numbers), A295421 (dodecahedral numbers), A175578 (icosahedral numbers).
|
|
KEYWORD
|
|
|
AUTHOR
|
|
|
EXTENSIONS
|
Quotation from Stan Wagon corrected by N. J. A. Sloane on Dec 24 2005. Thanks to Jose Brox for noticing this error.
|
|
STATUS
|
approved
|
|
|
|
|
A013661
|
|
Decimal expansion of Pi^2/6 = zeta(2) = Sum_{m>=1} 1/m^2.
|
|
+10
374
|
|
|
1, 6, 4, 4, 9, 3, 4, 0, 6, 6, 8, 4, 8, 2, 2, 6, 4, 3, 6, 4, 7, 2, 4, 1, 5, 1, 6, 6, 6, 4, 6, 0, 2, 5, 1, 8, 9, 2, 1, 8, 9, 4, 9, 9, 0, 1, 2, 0, 6, 7, 9, 8, 4, 3, 7, 7, 3, 5, 5, 5, 8, 2, 2, 9, 3, 7, 0, 0, 0, 7, 4, 7, 0, 4, 0, 3, 2, 0, 0, 8, 7, 3, 8, 3, 3, 6, 2, 8, 9, 0, 0, 6, 1, 9, 7, 5, 8, 7, 0
(list;
constant;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
"In 1736 he [Leonard Euler, 1707-1783] discovered the limit to the infinite series, Sum 1/n^2. He did it by doing some rather ingenious mathematics using trigonometric functions that proved the series summed to exactly Pi^2/6. How can this be? ... This demonstrates one of the most startling characteristics of mathematics - the interconnectedness of, seemingly, unrelated ideas." - Clawson [See Hardy and Wright, Theorems 332 and 333. - N. J. A. Sloane, Jan 20 2017]
Also Integral_{x>=0} x/(exp(x)-1) dx. [Abramowitz-Stegun, 23.2.7., for s=2, p. 807]
Pi^2/6 is also the length of the circumference of a circle whose diameter equals the ratio of volume of an ellipsoid to the circumscribed cuboid. Pi^2/6 is also the length of the circumference of a circle whose diameter equals the ratio of surface area of a sphere to the circumscribed cube. - Omar E. Pol, Oct 07 2011
Volume of a sphere inscribed in a cube of volume Pi. More generally, Pi^x/6 is the volume of an ellipsoid inscribed in a cuboid of volume Pi^(x-1). - Omar E. Pol, Feb 17 2016
Surface area of a sphere inscribed in a cube of surface area Pi. More generally, Pi^x/6 is the surface area of a sphere inscribed in a cube of surface area Pi^(x-1). - Omar E. Pol, Feb 19 2016
zeta(2)+1 is a weighted average of the integers, n > 2, using zeta(n)-1 as the weights for each n. We have: Sum_{n >= 2} (zeta(n)-1) = 1 and Sum_{n >= 2} n*(zeta(n)-1) = zeta(2)+1. - Richard R. Forberg, Jul 14 2016
zeta(2) is the expected value of sigma(n)/n. - Charlie Neder, Oct 22 2018
Graham shows that a rational number x can be expressed as a finite sum of reciprocals of distinct squares if and only if x is in [0, Pi^2/6-1) U [1, Pi^2/6). See section 4 for other results and Theorem 5 for the underlying principle. - Charles R Greathouse IV, Aug 04 2020
|
|
REFERENCES
|
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 811.
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.
Calvin C. Clawson, Mathematical Mysteries, The Beauty and Magic of Numbers, Perseus Books, 1996, p. 97.
W. Dunham, Euler: The Master of Us All, The Mathematical Association of America, Washington, D.C., 1999, p. xxii.
Hardy and Wright, 'An Introduction to the Theory of Numbers'. See Theorems 332 and 333.
A. A. Markoff, Mémoire sur la transformation de séries peu convergentes en séries très convergentes, Mém. de l'Acad. Imp. Sci. de St. Pétersbourg, XXXVII, 1890.
G. F. Simmons, Calculus Gems, Section B.15, B.24, pp. 270-271, 323-325, McGraw Hill, 1992.
Arnold Walfisz, Weylsche Exponentialsummen in der neueren Zahlentheorie, Deutscher Verlag der Wissenschaften, Berlin, 1963, p. 99, Satz 1.
A. Weil, Number theory: an approach through history; from Hammurapi to Legendre, Birkhäuser, Boston, 1984; see p. 261.
David Wells, The Penguin Dictionary of Curious and Interesting Numbers, Revised Edition, Penguin Books, London, England, 1997, page 23.
|
|
LINKS
|
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
|
|
FORMULA
|
Limit_{n->oo} (1/n)*(Sum_{k=1..n} frac((n/k)^(1/2))) = zeta(2) and in general we have lim_{n->oo} (1/n)*(Sum_{k=1..n} frac((n/k)^(1/m))) = zeta(m), m >= 2. - Yalcin Aktar, Jul 14 2005
Equals Integral_{x=0..1} (log(x)/(x-1)) dx or Integral_{x>=1} (log(x/(x-1))/x) dx. - Jean-François Alcover, May 30 2013
For s >= 2 (including Complex), zeta(s) = Product_{n >= 1} prime(n)^s/(prime(n)^s - 1). - Fred Daniel Kline, Apr 10 2014
zeta(2) = Sum_{n>=1} ((floor(sqrt(n)) - floor(sqrt(n-1)))/n). - Mikael Aaltonen, Jan 10 2015
zeta(2) = Sum_{n>=1} (((sqrt(5)-1)/2/sqrt(5))^n/n^2) + Sum_{n>=1} (((sqrt(5)+1)/2/sqrt(5))^n/ n^2) + log((sqrt(5)-1)/2/sqrt(5))log((sqrt(5)+1)/2/sqrt(5)). - Seiichi Kirikami, Oct 14 2015
The above formula can also be written zeta(2) = dilog(x) + dilog(y) + log(x)*log(y) where x = (1-1/sqrt(5))/2 and y=(1+1/sqrt(5))/2. - Peter Luschny, Oct 16 2015
zeta(2) = Integral_{x>=0} 1/(1 + e^x^(1/2)) dx, because (1 - 1/2^(s-1))*Gamma[1 + s]*Zeta[s] = Integral_{x>=0} 1/(1 + e^x^(1/s)) dx. After Jean-François Alcover in A002162. - Mats Granvik, Sep 12 2016
zeta(2) = Product_{n >= 1} (144*n^4)/(144*n^4 - 40*n^2 + 1). - Fred Daniel Kline, Oct 29 2016
Equals ((m+1)/m) * Integral_{x=0..1} log(Sum _{k=0..m} x^k )/x dx, m > 0 (Aubonnet reference). - Bernard Schott, Feb 11 2022
|
|
EXAMPLE
|
1.6449340668482264364724151666460251892189499012067984377355582293700074704032...
|
|
MAPLE
|
# Calculates an approximation with n exact decimal places (small deviation
# in the last digits are possible). Goes back to ideas of A. A. Markoff 1890.
zeta2 := proc(n) local q, s, w, v, k; q := 0; s := 0; w := 1; v := 4;
for k from 2 by 2 to 7*n/2 do
w := w*v/k;
q := q + v;
v := v + 8;
s := s + 1/(w*q);
od; 12*s; evalf[n](%) end:
|
|
MATHEMATICA
|
RealDigits[N[Pi^2/6, 100]][[1]]
|
|
PROG
|
(PARI) default(realprecision, 200); Pi^2/6
(PARI) default(realprecision, 200); dilog(1)
(PARI) default(realprecision, 200); zeta(2)
(PARI) A013661(n)={localprec(n+2); Pi^2/.6\10^n%10} \\ Corrected and improved by M. F. Hasler, Apr 20 2021
(PARI) default(realprecision, 20080); x=Pi^2/6; for (n=1, 20000, d=floor(x); x=(x-d)*10; write("b013661.txt", n, " ", d)); \\ Harry J. Smith, Apr 29 2009
(Maxima) fpprec : 100$ ev(bfloat(zeta(2)))$ bfloat(%); /* Martin Ettl, Oct 21 2012 */
(Magma) pi:=Pi(RealField(110)); Reverse(Intseq(Floor(10^105*pi^2/6))); // Vincenzo Librandi, Oct 13 2015
(Python) # Use some guard digits when computing.
# BBP formula (3 / 16) P(2, 64, 6, (16, -24, -8, -6, 1, 0)).
from decimal import Decimal as dec, getcontext
def BBPzeta2(n: int) -> dec:
getcontext().prec = n
s = dec(0); f = dec(1); g = dec(64)
for k in range(int(n * 0.5536546824812272) + 1):
sixk = dec(6 * k)
s += f * ( dec(16) / (sixk + 1) ** 2 - dec(24) / (sixk + 2) ** 2
- dec(8) / (sixk + 3) ** 2 - dec(6) / (sixk + 4) ** 2
+ dec(1) / (sixk + 5) ** 2 )
f /= g
return (s * dec(3)) / dec(16)
|
|
CROSSREFS
|
|
|
KEYWORD
|
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|
|
A046951
|
|
a(n) is the number of squares dividing n.
|
|
+10
130
|
|
|
1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 2, 1, 1, 1, 3, 1, 2, 1, 2, 1, 1, 1, 2, 2, 1, 2, 2, 1, 1, 1, 3, 1, 1, 1, 4, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 3, 2, 2, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1, 1, 2, 4, 1, 1, 1, 2, 1, 1, 1, 4, 1, 1, 2, 2, 1, 1, 1, 3, 3, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 3, 1, 2, 2, 4, 1, 1, 1, 2, 1, 1, 1, 4, 1, 1, 1, 3, 1, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,4
|
|
COMMENTS
|
Rediscovered by the HR automatic theory formation program.
a(n) depends only on prime signature of n (cf. A025487, A046523). So a(24) = a(375) since 24 = 2^3*3 and 375 = 3*5^3 both have prime signature (3, 1).
Number of 2-generated Abelian groups of order n, if n > 1. - Álvar Ibeas, Dec 22 2014 [In other words, number of order-n abelian groups with rank <= 2. Proof: let b(n) be such number. A finite abelian group is the inner direct product of all Sylow-p subgroups, so {b(n)} is multiplicative. Obviously b(p^e) = floor(e/2)+1 (corresponding to the groups C_(p^r) X C_(p^(e-r)) for 0 <= r <= floor(e/2)), hence b(n) = a(n) for all n. - Jianing Song, Nov 05 2022]
Number of ways of writing n = r*s such that r|s. - Eric M. Schmidt, Jan 08 2015
The number of divisors of the square root of the largest square dividing n. - Amiram Eldar, Jul 07 2020
|
|
LINKS
|
|
|
FORMULA
|
a(n) = |{(i, j) : i*j = n AND i|j}| = |{(i, j) : i*j^2 = n}|. Also tau(A000188(n)), where tau = A000005.
Inverse Moebius transform of characteristic function of squares (A010052). Dirichlet g.f.: zeta(s)*zeta(2s).
a(n) = Sum_{k = 1..n} ( floor(n/k^2) - floor((n-1)/k^2) ). - Peter Bala, Feb 17 2014
(End)
G.f.: Sum_{k>0}(theta(q^k)-1)/2, where theta(q)=1+2q+2q^4+2q^9+2q^16+... - Mamuka Jibladze, Dec 04 2016
(End)
L.g.f.: -log(Product_{k>=1} (1 - x^(k^2))^(1/k^2)) = Sum_{n>=1} a(n)*x^n/n. - Ilya Gutkovskiy, Jul 30 2018
a(n) = Sum_{d divides n} mu(core(d)^2), where core(n) = A007913(n). - Peter Bala, Jan 24 2024
|
|
EXAMPLE
|
a(16) = 3 because the squares 1, 4, and 16 divide 16.
G.f. = x + x^2 + x^3 + 2*x^4 + x^5 + x^6 + x^7 + 2*x^8 + 2*x^9 + x^10 + ...
|
|
MAPLE
|
local a, s;
a := 1 ;
for p in ifactors(n)[2] do
a := a*(1+floor(op(2, p)/2)) ;
end do:
a ;
|
|
MATHEMATICA
|
a[n_] := Length[ Select[ Divisors[n], IntegerQ[Sqrt[#]]& ] ]; Table[a[n], {n, 1, 105}] (* Jean-François Alcover, Jun 26 2012 *)
Table[Length[Intersection[Divisors[n], Range[10]^2]], {n, 100}] (* Alonso del Arte, Dec 10 2012 *)
a[ n_] := If[ n < 1, 0, Sum[ Mod[ DivisorSigma[ 0, d], 2], {d, Divisors @ n}]]; (* Michael Somos, Jun 13 2014 *)
a[ n_] := If[ n < 2, Boole[ n == 1], Times @@ (Quotient[ #[[2]], 2] + 1 & /@ FactorInteger @ n)]; (* Michael Somos, Jun 13 2014 *)
a[ n_] := If[ n < 0, 0, SeriesCoefficient[ Sum[ x^k^2 / (1 - x^k^2), {k, Sqrt @ n}], {x, 0, n}]]; (* Michael Somos, Jun 13 2014 *)
f[p_, e_] := 1 + Floor[e/2]; a[1] = 1; a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100] (* Amiram Eldar, Sep 15 2020 *)
|
|
PROG
|
(PARI) a(n)=my(f=factor(n)); for(i=1, #f[, 1], f[i, 2]\=2); numdiv(factorback(f)) \\ Charles R Greathouse IV, Dec 11 2012
(PARI) a(n) = direuler(p=2, n, 1/((1-X^2)*(1-X)))[n]; \\ Michel Marcus, Mar 08 2015
(Haskell)
a046951 = sum . map a010052 . a027750_row
(Scheme)
(define (A008619 n) (+ 1 (/ (- n (modulo n 2)) 2)))
(Magma) [#[d: d in Divisors(n)|IsSquare(d)]:n in [1..120]]; // Marius A. Burtea, Jan 21 2020
(Python)
from math import prod
from sympy import factorint
def A046951(n): return prod((e>>1)+1 for e in factorint(n).values()) # Chai Wah Wu, Aug 04 2024
|
|
CROSSREFS
|
Cf. A000005, A000188, A004101, A005117 (positions of 1's), A008619, A008833, A013936 (partial sums), A038538, A046952, A052304, A056595, A159631, A007814, A010052, A027750, A239930, A007862, A046523, A064989, A065704, A130279, A156552, A278161.
Differs from A096309 for the first time at n=32, where a(32) = 3, while A096309(32) = 2 (and also A185102(32) = 2).
Sum of the k-th powers of the square divisors of n for k=0..10: this sequence (k=0), A035316 (k=1), A351307 (k=2), A351308 (k=3), A351309 (k=4), A351310 (k=5), A351311 (k=6), A351313 (k=7), A351314 (k=8), A351315 (k=9), A351315 (k=10).
Sequences of the form n^k * Sum_{d^2|n} 1/d^k for k = 0..10: this sequence (k=0), A340774 (k=1), A351600 (k=2), A351601 (k=3), A351602 (k=4), A351603 (k=5), A351604 (k=6), A351605 (k=7), A351606 (k=8), A351607 (k=9), A351608 (k=10).
|
|
KEYWORD
|
nice,nonn,mult
|
|
AUTHOR
|
Simon Colton (simonco(AT)cs.york.ac.uk)
|
|
EXTENSIONS
|
Data section filled up to 125 terms and wrong claim deleted from Crossrefs section by Antti Karttunen, Nov 14 2016
|
|
STATUS
|
approved
|
|
|
|
|
A013631
|
|
Continued fraction for zeta(3).
|
|
+10
31
|
|
|
1, 4, 1, 18, 1, 1, 1, 4, 1, 9, 9, 2, 1, 1, 1, 2, 7, 1, 1, 7, 11, 1, 1, 1, 3, 1, 6, 1, 30, 1, 4, 1, 1, 4, 1, 3, 1, 2, 7, 1, 3, 1, 2, 2, 1, 16, 1, 1, 3, 3, 1, 2, 2, 1, 6, 1, 1, 1, 6, 1, 1, 4, 428, 5, 1, 1, 3, 1, 1, 11, 2, 4, 4, 5, 4, 1, 5, 14, 1, 3, 1, 2, 19, 1, 2, 5, 1, 7, 1, 1, 1, 1, 1, 57, 3, 2, 14, 2
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
|
|
EXAMPLE
|
zeta(3) = 1.2020569031595942... = 1 + 1/(4 + 1/(1 + 1/(18 + 1/(1 + ...)))). - Harry J. Smith, Apr 20 2009
|
|
MATHEMATICA
|
ContinuedFraction[ Zeta[3], 100]
|
|
PROG
|
(PARI) { allocatemem(932245000); default(realprecision, 21000); x=contfrac(zeta(3)); for (n=1, 20000, write("b013631.txt", n-1, " ", x[n])); } \\ Harry J. Smith, Apr 20 2009
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,cofr,nice
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|
|
A013680
|
|
Continued fraction for zeta(4).
|
|
+10
23
|
|
|
1, 12, 6, 1, 3, 1, 4, 183, 1, 1, 2, 1, 3, 1, 1, 5, 4, 2, 7, 23, 1, 1, 1, 1, 3, 2, 4, 2, 2, 22, 1, 13, 5, 1, 4, 2, 1, 3, 1, 1, 1, 6, 11, 40, 1, 7, 5, 2, 4, 1, 2, 3, 14, 9, 1, 33, 78, 1, 12, 4, 1, 2, 551, 1, 1, 1, 1, 1, 1, 2, 1, 9, 2, 7, 3, 1, 3, 2, 15, 1, 1, 2, 2
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
|
|
EXAMPLE
|
zeta(4) = 1 + 1/(12 + 1/(6 + 1/(1 + 1/(3 + ...)))). - Harry J. Smith, Apr 29 2009
|
|
MATHEMATICA
|
|
|
PROG
|
(PARI) { allocatemem(932245000); default(realprecision, 21000); x=contfrac(Pi^4/90); for (n=1, 20000, write("b013680.txt", n-1, " ", x[n])); } \\ Harry J. Smith, Apr 29 2009
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,cofr
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|
|
A222171
|
|
Decimal expansion of Pi^2/24.
|
|
+10
20
|
|
|
4, 1, 1, 2, 3, 3, 5, 1, 6, 7, 1, 2, 0, 5, 6, 6, 0, 9, 1, 1, 8, 1, 0, 3, 7, 9, 1, 6, 6, 1, 5, 0, 6, 2, 9, 7, 3, 0, 4, 7, 3, 7, 4, 7, 5, 3, 0, 1, 6, 9, 9, 6, 0, 9, 4, 3, 3, 8, 8, 9, 5, 5, 7, 3, 4, 2, 5, 0, 1, 8, 6, 7, 6, 0, 0, 8, 0, 0, 2, 1, 8, 4, 5, 8, 4, 0, 7, 2, 2, 5, 1, 5, 4, 9, 3, 9, 6, 7, 6, 3
(list;
constant;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,1
|
|
REFERENCES
|
George Boros and Victor H. Moll, Irresistible integrals, Cambridge University Press, 2006, p. 242.
Ovidiu Furdui, Limits, Series, and Fractional Part Integrals: Problems in Mathematical Analysis, New York: Springer, 2013. See Problem 3.45, p. 158 and 199-200.
|
|
LINKS
|
|
|
FORMULA
|
Equals Integral_{x=0..Pi/2} log(sec(x))/tan(x) dx.
Equals (1/10) * Sum_{k>=1} d(k^2)/k^2, where d(k) is the number of divisors of k (A000005). - Amiram Eldar, Jun 27 2020
Equals Sum_{n >= 0} 1/((2*n+1)*(6*n+3)). - Peter Bala, Feb 02 2022
Equals Sum_{n>=0} ((-1)^n * (Sum_{k>=n+1} (-1)^k/k)^2) (Furdui, 2013). - Amiram Eldar, Mar 26 2022
|
|
EXAMPLE
|
0.411233516712056609118103791661506297304737475301699609433889557342501867600...
|
|
MATHEMATICA
|
RealDigits[Pi^2/24, 10, 100] // First
|
|
PROG
|
(Magma) pi:=Pi(RealField(110)); Reverse(Intseq(Floor(10^100*(pi)^2/24))); // Vincenzo Librandi, Sep 25 2015
|
|
CROSSREFS
|
|
|
KEYWORD
|
|
|
AUTHOR
|
|
|
EXTENSIONS
|
Leading 0 term removed (to make offset correct) by Rick L. Shepherd, Jan 01 2014
|
|
STATUS
|
approved
|
|
|
|
|
A013696
|
|
Continued fraction for zeta(20).
|
|
+10
18
|
|
|
1, 1048259, 1, 2, 1, 18, 3, 1, 9, 7, 1, 1, 2, 1, 13, 3, 1, 1, 1, 2, 4, 2, 10, 2, 1, 1, 2, 8, 1, 1, 1, 3, 1, 3, 9, 2, 1, 2, 1, 1, 4, 2, 2, 56, 2, 2, 1, 1, 1, 6, 5, 2, 15, 1, 5, 2, 2, 1, 5, 1, 1, 39, 1, 6, 2, 6, 1, 1, 1, 3, 24, 11, 1, 1, 4
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
|
|
MATHEMATICA
|
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,cofr
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|
|
A013681
|
|
Continued fraction for zeta(5).
|
|
+10
5
|
|
|
1, 27, 12, 1, 1, 15, 1, 5, 1, 2, 19, 1, 1, 32, 1, 13, 1, 1, 1, 3, 1, 3, 2, 16, 1, 12, 4, 1, 5, 1, 1, 1, 1, 1, 2, 2, 6, 1, 8, 8, 6, 2, 3, 2, 2, 1, 30, 1, 17, 116, 1, 7, 1, 1, 1, 1, 1, 1, 2, 2, 12, 1, 4, 1, 1, 94, 1, 1, 3, 3, 6, 6, 1, 1, 2, 1, 17, 1, 1, 7, 1, 1, 11
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
|
|
MATHEMATICA
|
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,cofr
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|
|
A013682
|
|
Continued fraction for zeta(6).
|
|
+10
2
|
|
|
1, 57, 1, 1, 1, 15, 1, 6, 3, 61, 1, 5, 3, 1, 6, 1, 3, 3, 6, 1, 10, 1, 3, 2, 1, 4, 1, 1, 5, 1, 61, 1, 3, 1, 2, 1, 3, 2, 1, 3, 1, 2, 2, 28, 1, 2, 18, 53, 2, 1, 17, 11, 3, 4, 3, 5, 2, 1, 27, 9, 8, 3, 3, 3, 9, 5, 1, 3, 29, 1, 4, 1, 2, 40, 4, 8, 1, 3, 1, 2, 2, 1, 4
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
|
|
MATHEMATICA
|
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,cofr
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|
|
A013683
|
|
Continued fraction for zeta(7).
|
|
+10
2
|
|
|
1, 119, 1, 3, 2, 1, 2, 1, 39, 2, 3, 12, 3, 1, 1, 1, 2, 6, 5, 1, 5, 1, 2, 1, 23, 2, 1, 5, 34, 2, 1, 1, 3, 47, 2, 1, 8, 16, 1, 4, 1, 2, 1, 1, 1, 10, 72, 1, 1, 1, 1, 1, 2, 3, 13, 1, 2, 1, 5, 1, 27, 2, 9283, 1, 36, 1, 1, 1, 1, 3, 3, 23, 27, 5, 2, 4, 1, 3, 16, 1, 4
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
|
|
MATHEMATICA
|
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,cofr
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
Search completed in 0.018 seconds
|