from math import sqrt
def d(n):
s = 1
t = sqrt(n)
for i in range(2, int(t)+1):
if n % i == 0: s += i + n / i
if t == int(t): s -= t
return s
limit = 10000
sum = 0
for a in range(2, limit):
d_a = d(a)
if d_a > a and d(d_a) == a: sum += d_a + a
print "Answer to PE21 = ",sum
ZnJvbSBtYXRoIGltcG9ydCBzcXJ0CiAKZGVmIGQobik6CiAgcyA9IDEKICB0ID0gc3FydChuKQogIGZvciBpIGluIHJhbmdlKDIsIGludCh0KSsxKToKICAgIGlmIG4gJSBpID09IDA6IHMgKz0gaSArIG4gLyBpCiAgaWYgdCA9PSBpbnQodCk6IHMgLT0gdAogIHJldHVybiBzCiAKbGltaXQgPSAxMDAwMApzdW0gPSAwCmZvciBhIGluIHJhbmdlKDIsIGxpbWl0KToKICBkX2EgPSBkKGEpCiAgaWYgZF9hID4gYSBhbmQgZChkX2EpID09IGE6IHN1bSArPSBkX2EgKyBhIAogCnByaW50ICJBbnN3ZXIgdG8gUEUyMSA9ICIsc3Vt