A179890 Values of antiharmonic mean B(q) of the numbers k < q such that gcd(k, q) = 1 is an integer for nonprimes q from A179887, where B(q) = A053818(q) / A023896(q) = A175505(q) / A175506(q).

1, 7, 15, 23, 31, 37, 39, 55, 57, 61, 63, 71, 73, 77, 79, 89, 95, 97, 111, 113, 119, 121, 125, 135, 137, 143, 145, 151, 153, 157, 165, 169, 173, 175, 177, 177, 183, 193, 197, 199, 201, 213, 223, 231, 237, 239, 249, 255, 261, 263, 269, 273, 277, 285, 289, 297, 301, 303, 303

Offset: 1

Author

Jaroslav Krizek, Jul 30 2010, Jul 31 2010

Keywords

• nonn

Examples

			a(6) = 37 because for A179887(6) = 55 holds: B(55) = A053818(55)/A023896(55) = 40700/1100 = 37.
		

Links

• Amiram Eldar, Table of n, a(n) for n = 1..10000

Crossrefs

Cf. A023896, A053818, A175505, A175506.

Cf. A179871, A179872, A179873, A179874, A179875, A179876, A179877, A179878, A179879, A179880, A179882, A179883, A179884, A179885, A179886, A179887, A179890, A179891.

Programs

PARI

list(lim) = print1(1, ", "); for(k = 2, lim, if(!isprime(k), my(f = factor(k), b = 2*k/3 + (1/3) * prod(i = 1, #f~, 1 - f[i, 1])/eulerphi(f)); if(denominator(b) == 1, print1(b, ", ")))); \\ Amiram Eldar, May 26 2025

Extensions

Incorrect formula removed by Amiram Eldar, May 26 2025