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.

A254196 a(n) is the numerator of Product_{i=1..n} (1/(1-1/prime(i))) - 1.

This page as a plain text file.
%I A254196 #53 Feb 16 2025 08:33:24
%S A254196 1,2,11,27,61,809,13945,268027,565447,2358365,73551683,2734683311,
%T A254196 112599773191,4860900544813,9968041656757,40762420985117,
%U A254196 83151858555707,5085105491885327,341472595155548909,24295409051193284539
%N A254196 a(n) is the numerator of Product_{i=1..n} (1/(1-1/prime(i))) - 1.
%C A254196 The denominators are A038110(n+1).
%C A254196 a(n)/A038110(n+1) = Sum_{k >=2} 1/k where k is a positive integer whose prime factors are among the first n primes. In particular, for n=1,2,3,4,5, a(n)/A038110(n+1) is the sum of the reciprocals of the terms (excepting the first, 1) in A000079, A003586, A051037, A002473, A051038.
%C A254196 Appears to be a duplicate of A161527. - _Michel Marcus_, Aug 05 2019
%H A254196 Robert Israel, <a href="/A254196/b254196.txt">Table of n, a(n) for n = 1..422</a>
%H A254196 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/SmoothNumber.html">Smooth Number</a>
%F A254196 a(n) = A038111(n+1)/prime(n+1)-A038110(n+1). - _Robert Israel_, Jan 28 2015, corrected Jul 07 2019.
%e A254196 a(1)=1 because 1/2 + 1/4 + 1/8 + 1/16 + ... = 1/1.
%e A254196 a(2)=2 because 1/2 + 1/3 + 1/4 + 1/6 + 1/8 + 1/9 + 1/12 + ... = 2/1.
%e A254196 a(3)=11 because 1/2 + 1/3 + 1/4 + 1/5 + 1/6 + 1/8 + 1/9 + 1/10 + 1/12 + 1/15 + ... = 11/4.
%e A254196 a(4)=27 because Sum_{n>=2} 1/A002473(n) = 27/8.
%e A254196 a(5)=61 because Sum_{n>=2} 1/A051038(n) = 61/16.
%p A254196 seq(numer(mul(1/(1-1/ithprime(i)),i=1..n)-1),n=1..20); # _Robert Israel_, Jan 28 2015
%t A254196 Numerator[Table[Product[1/(1 - 1/p), {p, Prime[Range[n]]}] - 1, {n,1,20}]]
%t A254196 b[0] := 0; b[n_] := b[n - 1] + (1 - b[n - 1]) / Prime[n]
%t A254196 Numerator@ Table[b[n], {n, 1, 20}] (* _Fred Daniel Kline_, Jun 27 2017 *)
%o A254196 (PARI) a(n) = numerator(prod(i=1, n, (1/(1-1/prime(i)))) - 1); \\ _Michel Marcus_, Jun 29 2017
%Y A254196 Cf. A038110, A000079, A003586, A051037, A002473, A051038, A161527.
%K A254196 nonn
%O A254196 1,2
%A A254196 _Geoffrey Critzer_, Jan 26 2015