E para cada divisor que você encontrar, na verdade você encontrou - potencialmente - dois divisores. Por exemplo, se o número for 100 e você encontra o divisor 2, você também já encontrou o divisor 50 (resultado de 100 / 2). Então basta somar ambos, economizando uma iteração do loop. Só precisa tomar cuidado para o caso de quadrados perfeitos, para não contar duas vezes o mesmo divisor (por exemplo, se o número for 25100, não podemos usar esta lógica com o divisor 510, senão ele será contado duas vezes).
import math
def soma_divisores(num):
result = 1 + num
for i in range(2, int(math.sqrt(num)) + 1):
d, r = divmod(num, i)
if r == 0: # resto zero, é divisor
result += i;
if i != d: # somar também o outro divisor encontrado
result += d;
return result