OFFSET
1,1
MATHEMATICA
psi[n_] := n * Times @@ (1 + 1/FactorInteger[n][[;; , 1]]); psi[1] = 1; t[n_] := psi[n] - n; seq[max_] := Module[{s = {}, n}, Do[n = t[m]; If[n > m && t[n] == m, AppendTo[s, n - m]], {m, 1, max}]; s]; seq[200000] (* Amiram Eldar, Sep 05 2025 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
S. I. Dimitrov, Sep 04 2025
STATUS
approved
