39.py

import math
import time
start = time.time()
p = 120
l = []
maxl = 0
maxp = 0
for p in range(1001):
    l.append([])
    for a in range(1, int(p/2)):
        for b in range(a+1, p):
            c = math.sqrt(a**2+b**2)
            s = a+b+c
            if s > p:
                break
            elif s == p:
                l[p].append((a, b, int(c)))
                break
    print(p, " ", l[p])
    if len(l[p]) > maxl:
        maxl = len(l[p])
        maxp = p
print(maxp, l[maxp])
print(str(time.time()-start)[:7]+' s')