Sum of powers of digits

number-powers.py

import math

# replace 999999999 with upper limit of your choice
for i in range(1, 999999999):
    i_as_string = str(i)
    i_as_list = []
    j = 1
    output = 0
    for digits in i_as_string:
        output += pow(int(digits),j)
        j += 1
    if output == i:
        j = 0
        for digits in i_as_string:
            j += 1
            print(digits, "^", j, sep = "", end = "")
            if j < len(i_as_string):
                print(" + ", end = "")
        print (" = ", i)