OFFSET
2,1
LINKS
Michel Marcus, Table of n, a(n) for n = 2..500
FORMULA
For p prime, a(p) = p^2 (see A001248).
MATHEMATICA
prifacs[n_]:=If[n==1, {}, Flatten[Apply[ConstantArray, FactorInteger[n], {1}]]];
q=Table[Total[prifacs[n]]-Total[Union[prifacs[n]]], {n, 1000}];
mnrm[s_]:=If[Min@@s==1, mnrm[DeleteCases[s-1, 0]]+1, 0];
Table[Position[q, k][[1, 1]], {k, 2, mnrm[q/.(0->1)]}] (* Gus Wiseman, Feb 20 2025 *)
PROG
(PARI) sopfr(n) = my(f=factor(n)); sum(j=1, #f~, f[j, 1]*f[j, 2]);
sopf(n) = my(f=factor(n)); sum(j=1, #f~, f[j, 1]);
a(n) = {my(k = 2); while (sopfr(k) - sopf(k) != n, k++); k; }
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Dec 31 2016
STATUS
approved