OFFSET
1,2
COMMENTS
A066466 is a subsequence of this sequence.
LINKS
Dumitru Damian, Table of n, a(n) for n = 1..235 (terms up to 10^9)
EXAMPLE
n=317. Anti-divisors: 2, 3, 5, 127, 211.
Sum of the reversals of the anti-divisors: 2+3+5+721+112=843.
Sigma*(317)=348 and its reversal is 843.
n=1473. Anti-divisors: 2, 5, 6, 7, 19, 31, 95, 155, 421, 589, 982.
Sum of the reversals of the anti-divisors:
2+5+6+7+91+13+59+551+124+985+289=2132.
Sigma*(1473)=2312 and its reversal is 2132.
MAPLE
isA203616:=proc(j) local a, b, c; a:=0; b:=0;
for c from 2 to j-1 do
if abs((j mod c)-c/2)<1 then a:=a+A004086(c); b:=b+c; fi;
od;
evalb(A004086(b)=a) end: # simplified by M. F. Hasler, Jan 29 2012
for n to 10^7 do if isA203616(n) then lprint(n) fi od: # simplified by M. F. Hasler, Jan 29 2012
PROG
(Python)
from itertools import count, islice
from sympy.ntheory.factor_ import antidivisors
def a203616():
isa = lambda n: str(sum((a:=antidivisors(n))))[::-1]==str(sum(map(int, (str(_)[::-1] for _ in a))))
yield from (n for n in count(1) if isa(n))
a203616_list = [*islice(a203616(), 20)] # Dumitru Damian, Feb 12 2024
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Paolo P. Lava, Jan 20 2012
EXTENSIONS
a(22)-a(40) from Dumitru Damian, Feb 12 2024
STATUS
approved