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.

A370784 a(n) is the denominator of the sum of the reciprocals of the squarefree divisors of the powerful part of n.

This page as a plain text file.
%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