A175505 - OEIS

A175505

Numerator of A053818(n)/A023896(n) = antiharmonic mean of numbers k such that gcd(k,n) = 1, 1 <= k < n.

1, 1, 5, 5, 3, 13, 13, 21, 53, 7, 7, 49, 25, 29, 31, 85, 11, 109, 37, 27, 43, 15, 15, 193, 83, 53, 485, 113, 19, 59, 61, 341, 67, 23, 71, 433, 73, 77, 79, 107, 27, 83, 85, 59, 271, 31, 31, 769, 685, 167, 103, 209, 35, 973, 37, 449, 115, 39, 39, 239, 121, 125, 379, 1365

(list; graph; refs; listen; history; text; internal format)

OFFSET

1,3

COMMENTS

See A175506 - denominators of the antiharmonic means B of numbers k such that gcd(k, n) = 1 for numbers n >= 1 and k < n where B = A053818(n) / A023896(n) = a(n) / A175506(n).

LINKS

Table of n, a(n) for n=1..64.

Wikipedia, Contraharmonic mean.

FORMULA

a(n) = A053818(n) * A175506(n) / A023896(n).

Sum_{k=1..n} a(k)/175506(k) ~ n^2/3. - Amiram Eldar, Dec 07 2023

MAPLE

antiHMean := proc(L)

add(i^2, i=L)/add(i, i=L) ;

end proc:

A175505 := proc(n)

local kset, k ;

kset := [1] ;

for k from 2 to n do

if igcd(k, n) = 1 then

kset := [op(kset), k] ;

end if;

end do:

antiHMean(kset) ;

numer(%) ;

end proc: # R. J. Mathar, Sep 26 2013

MATHEMATICA

f[n_] := 2Plus @@ (Select[ Range@n, GCD[ #, n] == 1 &]^2)/(n*EulerPhi@n); f[1] = 1; Numerator@Array[f, 65] (* Robert G. Wilson v, Jul 01 2010 *)

PROG

(PARI) A175505(n)=numerator((2*n+(-1)^omega(n)*A007947(n)/n)/3) \\ M. F. Hasler, Nov 29 2010

(PARI) a(n) = {my(f = factor(n)); numerator(if(n == 1, 1, 2*n/3 + (1/3) * prod(i = 1, #f~, 1 - f[i, 1])/eulerphi(f))); } \\ Amiram Eldar, Dec 07 2023

CROSSREFS

Cf. A023896, A053818, A175506 (denominators).

Sequence in context: A232609 A225666 A365078 * A158274 A202695 A110986

Adjacent sequences: A175502 A175503 A175504 * A175506 A175507 A175508

KEYWORD

nonn,frac

AUTHOR

Jaroslav Krizek, May 31 2010, Jun 01 2010

EXTENSIONS

More terms from Robert G. Wilson v, Jul 01 2010

STATUS

approved