OFFSET
0,1
COMMENTS
a(n) is the least solution k to A346377(k) = n.
EXAMPLE
MAPLE
f:= proc(n) local t; add(t[1]*t[2], t=ifactors(n)[2])+n end proc:
N:= 10^6: # for terms <= N
V:= Vector(N):
for n from 1 to N do
v:= f(n);
if v <= N then V[v]:= V[v]+1 fi
od:
m:= max(V):
A:= Array(0..m):
for i from 1 to N do
if A[V[i]] = 0 then A[V[i]]:= i fi
od:
convert(A, list);
MATHEMATICA
f[1] = 1; f[n_] := n + Plus @@ Times @@@ FactorInteger[n]; m = 10^7; v = Table[0, {m}]; Do[i = f[n]; If[i <= m, v[[i]]++], {n, 1, m}]; TakeWhile[Table[ FirstPosition[v, k][[1]], {k, 0, Max[v]}], NumericQ] (* Amiram Eldar, Jul 14 2021 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Jul 14 2021
STATUS
approved