OFFSET
1,1
COMMENTS
EXAMPLE
Anti-divisors of 1656 are 7, 11, 16, 43, 48, 77, 144, 301, 368, 473, 1104 and their sum is 2592. Then, 711164348771443013684731104 / 2592 = 274368961717377705896887.
MAPLE
P:=proc(q) local a, b, k, n; for n from 3 to q do a:=0: b:=0:
for k from 2 to n-1 do if abs((n mod k)-k/2)<1 then
a:=a*10^length(k)+k: b:=b+k: fi; od; if frac(a/b)=0 then print(n);
fi; od; end: P(10^6);
PROG
(Python)
from itertools import islice, count
from sympy.ntheory.factor_ import antidivisors
def A308533gen(): # generator of terms
for n in count(3):
a = antidivisors(n)
if int(''.join(str(s) for s in a)) % sum(a) == 0:
yield n
A308533_list = list(islice(A308533gen(), 22)) # Chai Wah Wu, Dec 08 2021
CROSSREFS
KEYWORD
base,nonn,more
AUTHOR
Paolo P. Lava, Jun 06 2019
EXTENSIONS
a(22)-a(30) from Chai Wah Wu, Dec 08 2021
a(31)-a(32) from Chai Wah Wu, Dec 14 2021
STATUS
approved