OFFSET
1,2
COMMENTS
Each term appears a finite number of times in the sequence (Hagis and Lord, 1975).
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..290
Peter Hagis, Jr. and Graham Lord, Unitary harmonic numbers, Proc. Amer. Math. Soc., Vol. 51, No. 1 (1975), pp. 1-7.
MATHEMATICA
uh[n_] := n * Times @@ (2/(1 + Power @@@ FactorInteger[n])); uh[1] = 1; PrimeNu[Select[Range[10^6], IntegerQ[uh[#]] &]]
PROG
(PARI) uhmean(n) = {my(f = factor(n)); n*prod(i=1, #f~, 2/(1+f[i, 1]^f[i, 2])); };
lista(kmax) = {my(uh); for(k = 1, kmax, uh = uhmean(k); if(denominator(uh) == 1, print1(omega(k), ", "))); }
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Mar 10 2023
STATUS
approved