Euler #346

seen, value = [1], 1
for base in range(2, 10**6 + 1):
    value, exponent = 1 + base + base**2, 3
    while value <= 10**12:
        seen.append(value)
        value += base**exponent
        exponent += 1
print(sum(set(seen)))