%I #21 Nov 03 2023 17:01:58
%S 1,2,6,6,30,15,105,105,105,210,2310,2310,30030,15015,5005,5005,85085,
%T 85085,1616615,1616615,4849845,9699690,223092870,223092870,223092870,
%U 111546435,111546435,111546435,3234846615,2156564410,66853496710
%N Denominator of Sum_{k=1..n} mu(k)/k.
%H Amiram Eldar, <a href="/A070889/b070889.txt">Table of n, a(n) for n = 1..2370</a>
%e a(6) = 15 because 1 - 1/2 - 1/3 - 1/5 + 1/6 = 4/30 = 2/15.
%t Table[ Denominator[ Sum[ MoebiusMu[k]/k, {k, 1, n}]], {n, 1, 37}]
%o (PARI) t = 0; v = []; for( n = 1, 30, t = t + moebius( n) / n; v = concat( v, denominator( t))); v
%o (Python)
%o from functools import lru_cache
%o from sympy import harmonic
%o @lru_cache(maxsize=None)
%o def f(n):
%o if n <= 1:
%o return 1
%o c, j = 1, 2
%o k1 = n//j
%o while k1 > 1:
%o j2 = n//k1 + 1
%o c += (harmonic(j-1)-harmonic(j2-1))*f(k1)
%o j, k1 = j2, n//j2
%o return c+harmonic(j-1)-harmonic(n)
%o def A070889(n): return f(n).denominator # _Chai Wah Wu_, Nov 03 2023
%Y Cf. A008683, A068337, A070888 (numerators).
%K frac,nonn
%O 1,2
%A _Donald S. McDonald_, May 17 2002
%E Edited by _Robert G. Wilson v_, Jun 10 2002
|