题目39：如果p是直角三角形{a,b,c}的周长，1000以下的p中哪一个具有最多的解？

芒果加黄桃

1. # encoding:utf-8
   
2. # 寻找周长1000以下，能够成直角三角形最多的数
   
3. from time import time
   
4. def euler039():
   
5.     maxcount = 0
   
6.     result = []
   
7.     for p in range(12, 1001):
   
8.         count = 0
   
9.         tmp = []
   
10.         for i in range(int(p / (1 + 2 ** 0.5)) - 1, int(p / 2) + 1):
    
11.             for j in range(int(i / (2 ** 0.5)) - 1, i):
    
12.                 if i ** 2 == j ** 2 + (p - i - j) ** 2:
    
13.                     count += 1
    
14.                     tmp.append((p, i, j, (p - i - j)))
    
15.                     break
    
16.         if count > maxcount:
    
17.             result = tmp
    
18.             maxcount = count
    
19.     print(maxcount, result)
    
20. if __name__ == '__main__':
    
21.     start = time()
    
22.     euler039()
    
23.     print('cost %.6f sec' % (time() - start))
    

复制代码

8 [(840, 348, 252, 240), (840, 350, 280, 210), (840, 357, 315, 168), (840, 364, 336, 140), (840, 370, 350, 120), (840, 375, 360, 105), (840, 394, 390, 56), (840, 401, 399, 40)]
cost 5.641026 sec