OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..2500
MATHEMATICA
f[p_, e_] := p^e + e; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; Position[Partition[Array[s, 10^5], 2, 1], _?(SameQ @@ # &)] // Flatten
PROG
(PARI) s(n) = {my(f = factor(n)); prod(i = 1, #f~, f[i, 1]^f[i, 2] + f[i, 2]); }
lista(nmax) = {my(s1 = s(1), s2); for(n=2, nmax, s2=s(n); if(s1 == s2, print1(n-1, ", ")); s1 = s2); }
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Feb 10 2023
STATUS
approved