fork download
  1. from math import sqrt
  2.  
  3. def d(n):
  4. s = 1
  5. t = sqrt(n)
  6. for i in range(2, int(t)+1):
  7. if n % i == 0: s += i + n / i
  8. if t == int(t): s -= t
  9. return s
  10.  
  11. limit = 10000
  12. sum = 0
  13. for a in range(2, limit):
  14. d_a = d(a)
  15. if d_a > a and d(d_a) == a: sum += d_a + a
  16.  
  17. print "Answer to PE21 = ",sum
Success #stdin #stdout 0.28s 6516KB
stdin
Standard input is empty
stdout
Answer to PE21 =  31626