OFFSET
1,2
COMMENTS
See the comments in A382883.
MATHEMATICA
V[n_, e_] := If[e == 1, 1, IntegerExponent[n, e]]; f[n_] := f[n] = -DivisorSum[n, V[n, #] * f[#] &, # < n &]; f[1] = 1; Select[Range[100], f[#] != 0 &] (* Amiram Eldar, Apr 29 2025 *)
PROG
(SageMath)
def A383106List(upto): return [n for n in srange(1, upto) if A382883(n) != 0]
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Apr 16 2025
STATUS
approved
