OFFSET
1,2
COMMENTS
The terms are conjectural and based on a search for solutions to esigma(x) - x = k for k in the range of the data section and x < 10^12 (esigma(x) - x = A051377(x) - x = A126164(x) is the sum of aliquot exponential divisors of x). - Amiram Eldar, Jan 22 2020
MATHEMATICA
fun[p_, e_] := DivisorSum[e, p^# &]; esigma[1] = 1; esigma[n_] := esigma[n] = Times @@ fun @@@ FactorInteger[n]; untouchableQ[n_] := Catch[ Do[ If[n == esigma[k]-k, Throw[True]], {k, 0, (n+1)^2}]] === Null; Reap[ Table[ If[ untouchableQ[n], Sow[n]], {n, 1, 130}]][[2, 1]] (* after Jean-François Alcover at A005114 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Feb 20 2019
EXTENSIONS
Data corrected by Amiram Eldar, Jan 22 2020
STATUS
approved