OFFSET
1,1
COMMENTS
The palindromic primes (A002385) (whose sum of proper divisors is 1) are not included in this sequence.
LINKS
Donovan Johnson, Table of n, a(n) for n = 1..260 (terms < 10^14)
EXAMPLE
10801 = 1543*7*1 and 1543 + 7 + 1 = 1551 (a palindrome). So, 10801 is a member of this sequence.
PROG
(Python)
from sympy import divisors
def rev(n):
..r = ''
..for i in str(n):
....r = i + r
..return int(r)
for i in range(10**6):
..s = sum(divisors(i))-i
..if rev(i) == i and s > 10 and rev(s) == s:
....print(i, end=', ')
# simplified by Derek Orr, Oct 16 2014
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Derek Orr, Sep 19 2013
STATUS
approved