|
|
A295126
|
|
Denominator of Sum_{d|n} mu(n/d)/d, where mu is the Möbius function A008683.
|
|
2
|
|
|
1, 2, 3, 4, 5, 3, 7, 8, 9, 5, 11, 6, 13, 7, 15, 16, 17, 9, 19, 5, 7, 11, 23, 12, 25, 13, 27, 14, 29, 15, 31, 32, 33, 17, 35, 18, 37, 19, 13, 10, 41, 7, 43, 22, 45, 23, 47, 24, 49, 25, 51, 13, 53, 27, 11, 28, 19, 29, 59, 15, 61, 31, 21, 64, 65, 33, 67, 17, 69, 35
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
a(n) <= n.
n == 0 (mod a(n)).
First occurrence of k: 1, 2, 3, 4, 5, 12, 7, 8, 9, 40, 11, 24, 13, 28, 15, 16, 17, 36, 19, 80, 63, 44, 23, 48, 25, ..., ;
First occurrence of k = a(n)/n: 1, 6, 21, 20, 55, 42, 203, 120, 171, 110, 253, 84, 689, 406, 465, 272, 1751, 342, ..., .
|
|
LINKS
|
|
|
EXAMPLE
|
a(6) = 1 since mu(6)/1 + mu(3)/2 + mu(2)/3 + mu(1)/6 = 1 - 1/2 - 1/3 + 1/6 = 1/3.
|
|
MAPLE
|
f:= n -> denom(add(numtheory:-mobius(n/k)/k, k=numtheory:-divisors(n))):
|
|
MATHEMATICA
|
f[n_] := Block[{d = Divisors@ n}, Plus @@ (MoebiusMu[d]/Reverse@ d)]; Denominator@ Array[f, 70]
|
|
PROG
|
(PARI) a(n) = denominator(sumdiv(n, d, moebius(n/d)/d)); \\ Michel Marcus, Nov 17 2017
|
|
CROSSREFS
|
|
|
KEYWORD
|
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|