def divisors(n):
d = set([1, n])
i = 2
while i**2 <= n:
q, r = divmod(n, i)
if r == 0:
d.add(i)
d.add(q)
i += 1
return d
def condition(s):
return s % 100 // 10 == 7
n = 500001
k = 5
while k:
s = sum(divisors(n))
if condition(s):
print(n, s)
k -= 1
n += 1
ZGVmIGRpdmlzb3JzKG4pOgogICAgZCA9IHNldChbMSwgbl0pCiAgICBpID0gMgogICAgd2hpbGUgaSoqMiA8PSBuOgogICAgICAgIHEsIHIgPSBkaXZtb2QobiwgaSkKICAgICAgICBpZiByID09IDA6CiAgICAgICAgICAgIGQuYWRkKGkpCiAgICAgICAgICAgIGQuYWRkKHEpCiAgICAgICAgaSArPSAxCiAgICByZXR1cm4gZAoKZGVmIGNvbmRpdGlvbihzKToKICAgIHJldHVybiBzICUgMTAwIC8vIDEwID09IDcKCm4gPSA1MDAwMDEKayA9IDUKCndoaWxlIGs6CiAgICBzID0gc3VtKGRpdmlzb3JzKG4pKQogICAgaWYgY29uZGl0aW9uKHMpOgogICAgICAgIHByaW50KG4sIHMpCiAgICAgICAgayAtPSAxCiAgICBuICs9IDE=