OFFSET
1,1
COMMENTS
a(6) > 10^13. - Giovanni Resta, Jun 26 2015
EXAMPLE
sigma(5602 + 1) = 6048 = reverse(8406) = reverse(sigma(5602)), so 5602 is a term of the sequence.
MATHEMATICA
Do[ If[FromDigits[Reverse[IntegerDigits[DivisorSigma[1, n]]]] == DivisorSigma[1, n + 1], Print[n]], {n, 1, 10^7}]
PROG
(PARI) isok(n) = digits(sigma(n+1)) == Vecrev(digits(sigma(n))); \\ Michel Marcus, Jun 26 2015
CROSSREFS
KEYWORD
base,more,nonn
AUTHOR
Joseph L. Pe, Sep 19 2002
EXTENSIONS
a(3)-a(4) from Donovan Johnson, Feb 01 2009
a(5) from Giovanni Resta, Jun 26 2015
STATUS
approved