r=range(-50,51) r=sorted([a*a+b*b+c*c for a in r for b in r for c in r]) for d in range(51):print r.index(d*d+1)
Note that non-ascii characters in the above source code will be escaped (such as \x9f).
download
return to the top page