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.

A378681 a(n) = denominator(Sum_{k=1..n} 1/P_2(k)), where P_2(k) = A087040(k) is the second largest prime dividing the k-th composite number.

This page as a plain text file.
%I A378681 #7 Dec 06 2024 11:29:53
%S A378681 2,1,2,6,3,6,3,3,6,2,1,3,6,3,15,30,30,15,5,10,30,15,15,15,30,10,5,15,
%T A378681 30,30,15,30,210,210,70,35,105,105,210,210,105,35,70,210,105,21,21,42,
%U A378681 14,70,210,105,105,210,210,210,105,35,70,210,210,105,105,210
%N A378681 a(n) = denominator(Sum_{k=1..n} 1/P_2(k)), where P_2(k) = A087040(k) is the second largest prime dividing the k-th composite number.
%C A378681 See A378680 for more details.
%H A378681 Amiram Eldar, <a href="/A378681/b378681.txt">Table of n, a(n) for n = 1..1000</a>
%t A378681 p2[c_] := Module[{f = FactorInteger[c]}, If[f[[-1, 2]] > 1, f[[-1, 1]], f[[-2, 1]]]]; Denominator@ Accumulate[Table[1/p2[c], {c, Select[Range[50], CompositeQ]}]]
%o A378681 (PARI) lista(nmax) = {my(s = 0); forcomposite(n = 1, nmax, f = factor(n); s += if(f[#f~, 2] > 1, 1/f[#f~, 1], 1/f[#f~ - 1, 1]); print1(denominator(s), ", "));}
%Y A378681 Cf. A087039, A087040, A378680 (numerators).
%K A378681 nonn,easy,frac
%O A378681 1,1
%A A378681 _Amiram Eldar_, Dec 03 2024