Задача по теме: "Поиск делителей"

Решение:

Для нахождения всех подходящих значений напишем программу:

def prime(q):
    divs = []
    for d in range(1, round(q**0.5) + 1):
        if q % d == 0:
            divs += [d, q//d]
    divs = sorted(set(divs))[1:-1]
    return True if divs == [] else False

def s(x):
    divs = []
    for d in range(1, round(x**0.5) + 1):
        if x % d == 0:
            divs += [d, x//d]
    divs = sorted(set(divs))[1:-1]
    s1 = sum([el for el in divs if prime(el) is True])
    return s1

cnt = 0
for x in range(550001, 1000000):
    if str(s(x))[-1] == '7':
        print(x, s(x))
        cnt += 1
    if cnt == 5:
        break

Ответ: 550018 2337 550059 1197 550087 257 550099 1217 550137 117