OFFSET
0,2
COMMENTS
a(n) is the least k such that A386591(k) = n.
LINKS
David A. Corneth, Table of n, a(n) for n = 0..752 (first 251 terms from Robert Israel)
EXAMPLE
a(3) = 16 because 16 has exactly 3 divisors that are not balanced numbers, namely 4, 8 and 16, and no smaller number works.
MAPLE
g:= proc(n) option remember; numtheory:-sigma(n) mod numtheory:-phi(n) <> 0 end proc:
f:= n -> nops(select(g, numtheory:-divisors(n))):
N:= 60: # for a(0) to a(N)
V:= Array(0..N, -1): count:= 0:
for i from 1 while count < N+1 do
v:= f(i);
if V[v] = -1 then count:= count+1; V[v]:= i; fi
od:
convert(V, list);
PROG
(PARI) a(n) = my(k=1); while (sumdiv(k, d, sigma(d)%eulerphi(d) != 0) != n, k++); k; \\ Michel Marcus, Aug 26 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert Israel, Aug 26 2025
STATUS
approved
