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 A370784 #8 Mar 02 2024 03:50:45 %S A370784 1,1,1,2,1,1,1,2,3,1,1,2,1,1,1,2,1,3,1,2,1,1,1,2,5,1,3,2,1,1,1,2,1,1, %T A370784 1,1,1,1,1,2,1,1,1,2,3,1,1,2,7,5,1,2,1,3,1,2,1,1,1,2,1,1,3,2,1,1,1,2, %U A370784 1,1,1,1,1,1,5,2,1,1,1,2,3,1,1,2,1,1,1 %N A370784 a(n) is the denominator of the sum of the reciprocals of the squarefree divisors of the powerful part of n. %H A370784 Amiram Eldar, <a href="/A370784/b370784.txt">Table of n, a(n) for n = 1..10000</a> %H A370784 Rafael Jakimczuk, <a href="http://dx.doi.org/10.13140/RG.2.2.12174.13124">Arithmetical Functions over the Powerful Part of an Integer</a>, ResearchGate, 2024. %F A370784 a(n) = A332881(A057521(n)). %F A370784 a(n) = 1 if n is squarefree (A005117). %t A370784 a[n_] := Denominator[Times @@ (1 + 1/Select[FactorInteger[n], Last[#] > 1 &][[;; , 1]])]; Array[a, 100] %o A370784 (PARI) a(n) = {my(f = factor(n)); denominator(prod(i = 1, #f~, if(f[i,2] == 1, 1, 1 + 1/f[i,1])));} %Y A370784 Cf. A005117, A057521, A295295, A332881, A370783 (numerators). %K A370784 nonn,easy,frac %O A370784 1,4 %A A370784 _Amiram Eldar_, Mar 02 2024