Untitled

def list_squared(m, n):
    solution = []
    
    while m < n:
        divisors_sum = 0
        for i in range(1, m+1):
            if m % i == 0:
                divisors_sum += i**2
        if (divisors_sum**(1/2)).is_integer():
            solution.append([m, divisors_sum])
        m+=1
    
    return solution