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.

A357762 Decimal expansion of -Sum_{k>=1} A106400(k)/k.

This page as a plain text file.
%I A357762 #7 Oct 12 2022 11:39:44
%S A357762 1,1,9,6,2,8,3,2,6,4,3,2,5,2,5,6,4,3,7,2,2,2,2,9,1,6,3,3,2,0,0,8,1,9,
%T A357762 1,8,1,0,1,0,4,2,6,7,4,6,4,0,1,5,9,4,3,8,1,8,9,8,7,2,3,3,3,7,3,0,7,8,
%U A357762 3,7,5,1,6,1,0,9,1,5,8,0,8,7,7,7,9,1,1,9,6,4,5,4,6,2,1,1,0,7,4,8,9,6,3,3,3
%N A357762 Decimal expansion of -Sum_{k>=1} A106400(k)/k.
%C A357762 The asymptotic mean of the excess of the number of odious divisors over the number of evil divisors (A357761, see formula).
%C A357762 The convergence of the partial sums S(m) = -Sum_{k=1..2^m-1} A106400(k)/k is fast: e.g., S(28) is already correct to 100 decimal digits (see also _Jon E. Schoenfield_'s comment in A351404).
%F A357762 Equals -2 * Sum_{k>=1} A106400(2*k-1)/(2*k-1).
%F A357762 Equals lim_{m->oo} (1/m) * Sum_{k=1..m} A357761(k).
%e A357762 1.19628326432525643722229163320081918101042674640159...
%t A357762 sum = 0; m = 1; pow = 2; Do[sum -= (-1)^DigitCount[k, 2, 1]/k; If[k == pow - 1, Print[m, " ", N[sum, 120]]; m++; pow *= 2], {k, 1, 2^30}]
%o A357762 (PARI) default(realprecision, 150);
%o A357762 sm = 0.; m = 1; pow = 2; for(k = 1, 2^30, sm -= (-1)^hammingweight(k)/k; if(k == pow - 1, print(m," ",sm); m++; pow *= 2))
%Y A357762 Cf. A106400, A357761.
%Y A357762 Similar constants: A215016, A351404
%K A357762 nonn,cons
%O A357762 1,3
%A A357762 _Amiram Eldar_, Oct 12 2022