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.
%I A327340 #25 Feb 16 2025 08:33:58 %S A327340 1,1,8,7,4,8,40,13,64,41,94,59,132,39,4,51,222,43,278,157,346,191,406, %T A327340 227,484,263,562,305,640,178,24,99,280,447,942,169,1052,278,1168,31, %U A327340 1282,689,1422,747,58,819,1686,99,1838,482 %N A327340 Numerator of the rationals r(n) = (1/n^2)*Phi_1(n), with Phi_1(n) = Sum{k=1..n} psi(k), with Dedekind's psi function. %C A327340 The corresponding denominators are given in A327341. %C A327340 Dedekind's psi(k) = k*Product_{p|k}(1 + 1/p), with primes p, and the empty product is set to 1. See psi(k) = A001615(k), k >= 1. In the Walfisz reference psi(k) = phi_1(k). %C A327340 In the Walfisz reference, Satz 2., p. 100, the approximation for Phi_1(x) = (15/(2*Pi^2))*x^2 + O(x*(log(x))^{2/3}) is given (with B instead of the O() notation). For the constant 15/(2*Pi^2) see A323669 . %D A327340 Arnold Walfisz, Weylsche Exponentialsummen in der neueren Zahlentheorie, VEB Deutscher Verlag der Wissenschaften, Berlin, 1963, p. 100, Satz 2. %H A327340 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/DedekindFunction.html">Dedekind Function</a>. %H A327340 Wikipedia, <a href="https://en.wikipedia.org/wiki/Dedekind_psi_function">Dedekind psi function</a>. %F A327340 a(n) = numerator(r(n)), with the rationals r(n) = (1/n^2)*Sum{k=1..n}(k*Product_{p|k}(1 + 1/p)), with distinct prime p divisors of k (with empty product set to 1 for k = 1), for n >= 1. %F A327340 a(n) = numerator(A173290(n)/n^2). - _Amiram Eldar_, Nov 24 2022 %e A327340 The rationals (in lowest terms) begin: 1/1, 1/1, 8/9, 7/8, 4/5, 8/9, 40/49, 13/16, 64/81, 41/50, 94/121, 59/72, 132/169, 39/49, 4/5, 51/64, 222/289, 43/54, 278/361, 157/200, 346/441, 191/242, 406/529, 227/288, 484/625, 263/338, 562/729, 305/392, 640/841, 178/225, 24/31, ... %e A327340 The limit of r(n) for n-> infinity is A323669 = 0.759908877317533285829... %e A327340 r(10^5) is approximatly 0.7599142240 (10 digits). %t A327340 psi[0] = 1; psi[n_] := n * Times @@ (1 + 1/Transpose[FactorInteger[n]][[1]]); a[n_] := Numerator[Sum[psi[k], {k, 1, n}]/n^2]; Array[a, 50] (* _Amiram Eldar_, Sep 03 2019 *) %o A327340 (PARI) dpsi(n) = n * sumdivmult(n, d, issquarefree(d)/d); \\ A001615 %o A327340 a(n) = numerator(sum(k=1, n, dpsi(k))/n^2); \\ _Michel Marcus_, Sep 18 2023 %Y A327340 Cf. A001615, A173290, A323669, A327341 (denominators). %K A327340 nonn,frac,easy %O A327340 1,3 %A A327340 _Wolfdieter Lang_, Sep 03 2019