login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

Decimal expansion of the sum of the alternating series of reciprocals of composite numbers with distinct prime factors.
3

%I #34 Feb 26 2020 08:14:03

%S 1,1,6,5,9,1,0,2,4,8,4

%N Decimal expansion of the sum of the alternating series of reciprocals of composite numbers with distinct prime factors.

%C The denominators in the terms used to generate the sequence exclude the powers of primes (A000961).

%C From _Jon E. Schoenfield_, Apr 19 2018: (Start)

%C The k-th term in the sum is t(k) = (-1)^(k+1)/A024619(k). The partial sums S(j) = Sum_{k=1..j} t(k) converge rather slowly; e.g., since A024619(921265) = 10^6, each partial sum is still differing from its predecessor by an absolute difference of about 0.000001 after 921265 terms.

%C However, smoothing the sequence of partial sums by counting only half of the last term in each partial sum, i.e., S'(j) = (Sum_{k=1..j-1} t(k)) + (1/2)*t(j), yields a sequence that converges much more quickly. Faster still is the convergence of the sequence of partial sums S'(j) when only those values of j for which 1/t(j) is odd are considered, similar to the approach described at A275712. Here, however, the relative smoothness of the plot of S'(j), considering only those values of j for which 1/t(j) is odd, is interrupted by a significant step upward or downward each time the first such 1/t(j) is encountered that exceeds a new power of 2. (The absolute value of each such discontinuity is on the order of 1/(t(j))^2, and its sign depends on the parity of the number of terms in A024619 that are less than that power of 2.) In spite of these discontinuities, however, apparent convergence beyond the digits currently listed in the Data section is easy to obtain; see the table in the Example section.

%C (End)

%F Sum_{k>=1} (-1)^(k+1)/A024619(k) = 1/6 - 1/10 + 1/12 - 1/14 + 1/15 - 1/18 + ...

%e 0.11659102484...

%e From _Jon E. Schoenfield_, Apr 19 2018: (Start)

%e For b = 3..23, the table below lists the minimum and maximum values of the smoothed partial sum S'(j) (see Comments) among values of j such that 1/t(j) is an odd number in the interval [2^b, 2^(b+1)]. (The table begins at b=3 because the first odd denominator in the alternating sum is A024619(5) = 15, which lies in the interval [2^3, 2^4].)

%e .

%e b min S'(j) (1/t(j) odd) max S'(j) (1/t(j) odd)

%e == ======================= =======================

%e 3 0.111904761904761904... 0.111904761904761904...

%e 4 0.115873015873015873... 0.115873015873015873...

%e 5 0.116751763887986798... 0.116780441780441780...

%e 6 0.116515193298484059... 0.116521345165903716...

%e 7 0.116576739754580331... 0.116579123025777238...

%e 8 0.116593064475078642... 0.116593574635271689...

%e 9 0.116589740175207247... 0.116589805144898644...

%e 10 0.116590699983114246... 0.116590715183253815...

%e 11 0.116590952081202610... 0.116590954393007789...

%e 12 0.116591012248587507... 0.116591013407441019...

%e 13 0.116591027701733147... 0.116591027880912807...

%e 14 0.116591024064904443... 0.116591024080256316...

%e 15 0.116591025004372396... 0.116591025007229633...

%e 16 0.116591024772646678... 0.116591024773252190...

%e 17 0.116591024831104914... 0.116591024831253967...

%e 18 0.116591024845789528... 0.116591024845809974...

%e 19 0.116591024842154904... 0.116591024842161473...

%e 20 0.116591024841250054... 0.116591024841250922...

%e 21 0.116591024841022776... 0.116591024841022992...

%e 22 0.116591024840966103... 0.116591024840966133...

%e 23 0.116591024840980341... 0.116591024840980345...

%e 24 0.116591024840976791... 0.116591024840976792...

%e (End)

%o (Sage)

%o composites = (i for i in NN if i>3 and not i.is_prime() and not i.is_prime_power())

%o s = RLF(0); s

%o RealField(110)(s)

%o for i in range(0,50000000): s += (-1)^i / next(composites)

%o print(s) # _Terry D. Grant_, Aug 01 2016

%o (PARI) approx(lim)=my(s=0.,r=1,t); forcomposite(n=6,lim, if(ispower(n,,&t) && isprime(t), next); s += r/n; r=-r); t=floor(lim)+1; if(isprime(t), t++); t=s+r/t; [min(t,s), max(t,s)] \\ Gives lower and upper bounds

%o \\ _Charles R Greathouse IV_, Jul 31 2016

%Y Cf. A024619, A269229.

%K nonn,cons,more

%O 0,3

%A _Terry D. Grant_, Jul 16 2016

%E 2 more digits from _Charles R Greathouse IV_, Aug 02 2016