from math import ceil
def dsum(num):
res = 1
for k in range(2, int(num**.5) + 1):
if num % k == 0:
res += k + num//k
return res
my_list = [0]*101
my_set = set()
for i in range(1, 101):
for j in range(1000*i - 999, 1000*i+1):
if j not in my_set:
ssum1 = dsum(j)
if ssum1 < 100001:
ssum2 = dsum(ssum1)
if ssum2 == j and ssum1 != ssum2:
m = ceil(ssum1/1000)
my_list[i] += 1
my_list[m] += 1
my_set |= {ssum1, ssum2}
for i in range(1, 11):
print(f'({i*1000-999}, {i*1000}): {my_list[i]}')
ZnJvbSBtYXRoIGltcG9ydCBjZWlsCgpkZWYgZHN1bShudW0pOgogICAgcmVzID0gMQogICAgZm9yIGsgaW4gcmFuZ2UoMiwgaW50KG51bSoqLjUpICsgMSk6CiAgICAgICAgaWYgbnVtICUgayA9PSAwOgogICAgICAgICAgICByZXMgKz0gayArIG51bS8vawogICAgcmV0dXJuIHJlcwoKCm15X2xpc3QgPSBbMF0qMTAxCm15X3NldCA9IHNldCgpCmZvciBpIGluIHJhbmdlKDEsIDEwMSk6CgogICAgZm9yIGogaW4gcmFuZ2UoMTAwMCppIC0gOTk5LCAxMDAwKmkrMSk6CiAgICAgICAgaWYgaiBub3QgaW4gbXlfc2V0OgogICAgICAgICAgICBzc3VtMSA9IGRzdW0oaikKICAgICAgICAgICAgCiAgICAgICAgICAgIGlmIHNzdW0xIDwgMTAwMDAxOgogICAgICAgICAgICAgICAgc3N1bTIgPSBkc3VtKHNzdW0xKQogICAgICAgICAgICAgICAgaWYgc3N1bTIgPT0gaiBhbmQgc3N1bTEgIT0gc3N1bTI6CiAgICAgICAgICAgICAgICAgICAgbSA9IGNlaWwoc3N1bTEvMTAwMCkKICAgICAgICAgICAgICAgICAgICBteV9saXN0W2ldICs9IDEKICAgICAgICAgICAgICAgICAgICBteV9saXN0W21dICs9IDEKICAgICAgICAgICAgICAgICAgICBteV9zZXQgfD0ge3NzdW0xLCBzc3VtMn0KCmZvciBpIGluIHJhbmdlKDEsIDExKToKICAgIHByaW50KGYnKHtpKjEwMDAtOTk5fSwge2kqMTAwMH0pOiB7bXlfbGlzdFtpXX0nKQ==
(1, 1000): 2
(1001, 2000): 2
(2001, 3000): 2
(3001, 4000): 0
(4001, 5000): 0
(5001, 6000): 2
(6001, 7000): 2
(7001, 8000): 0
(8001, 9000): 0
(9001, 10000): 0