OFFSET
1,1
EXAMPLE
sigma(1302) - sigma(reverse(1302)) = sigma(1302) - sigma(2031) = 3072 - 2712 = 360 = phi(1302).
MATHEMATICA
r = {}; Do[ If[ DivisorSigma[1, n] - DivisorSigma[1, FromDigits[ Reverse[ IntegerDigits[n]]]] == EulerPhi[n], r = Append[r, n]], {n, 1, 10^7}]; r
Select[Range[81*10^5], DivisorSigma[1, #]-DivisorSigma[ 1, IntegerReverse[ #]] == EulerPhi[#]&] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Aug 08 2020 *)
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Joseph L. Pe, Jun 09 2002
EXTENSIONS
Edited and extended by Robert G. Wilson v, Jun 11 2002
STATUS
approved