Python program to find the sum of the even divisors of a number

Hello everyone, in this Python tutorial, we are going to write a python program to find the sum of the even divisors of a number. Suppose we have a number as 12. Its divisors are 1, 2, 3, 4, 6, 12. The even divisors are 2, 4, 6, 12 and their sum is 24 (2+4+6+12). Let’s have a look at steps we need to follow to write a Python program to find such a sum.

Find the sum of the even divisors of a number in Python

Let’s say our input number is n and its prime factors are p1, p2, …pk. Let a1, a2, …ak be the highest powers of corresponding prime factor. Now we can represent n as below:

n= p1^a1 * p2^a2 * ....pk^ak

We can easily infer from the above that the sum of all the divisors is :

sum of all factors=
(1 + p1 + p1^2 + ...p1^a1) *
(1 + p2 + p2^2 + ...p2^a2) *
.......................... *
(1 + pk + pk^2 + ...pk^ak)

If the number is odd, there will be no even divisors. Hence the sum will be zero. If the number has 2 as one of its prime factors, then the sum can be written as:

sum of all factors= 
(1 + 2 + 2^2 + ...2^a1) * 
(1 + p2 + p2^2 + ...p2^a2) * 
.......................... * 
(1 + pk + pk^2 + ...pk^ak)

Now in order to find the sum of even divisors we need to include only the cases when there is 2 in the factor. Therefore we exclude 1 (2^0) from the above from the second and the sum becomes:

sum of even factors= 
(2 + 2^2 + ...2^a1) * 
(1 + p2 + p2^2 + ...p2^a2) * 
.......................... * 
(1 + pk + pk^2 + ...pk^ak)

Let’s understand the above equation with an example.

Let’s say we have n=12.

The prime factors of n are 2, 3 and we can represent it as 2^2 * 3^1.

Now the sum of even divisors is (2 + 2^2) (1 + 3) = 24.

We can use property of Geometric Progression to find the sum (1 + p + p^2 + …p^a) and (2 + 2^2 + …2^a).

1 + p + p^2 + ...p^a = 1 (p^(a+1) - 1)/(p - 1) = (p^(a+1) - 1)/(p - 1)

Now we will implement it with an example program in Python. See the code.

def EvenDivisorsSum(n):
    
    if n%2 != 0:
        return 0
    if n == 2:
        return 2
    result = 1
    lim = n
    for i in range(2, lim):
        a = 0
        sum_i = 0
        while(n % i == 0):
            n = n // i
            a = a + 1
        sum_i = (i**(a+1)-1)/(i-1)
        if(i==2):
            sum_i = sum_i - 1
        result = result * sum_i
    return result
    
        
        
print(EvenDivisorsSum(30))

The output of the above code is:

48.0

Thank you.

Also, read: Sum of all the factors of a number in Python