从sympy导入integer_ntroot
A375750型,x个 k个 =设置(), 2[], 1
而len(A375750型) <100060:
对于范围（1，x）中的y：
c=（b:=k*4）+1
mid=integer_ntroot（（c_ttpo2:=c**2）//2，2）[0]
当b>mid时：
一 秒 =x个整数_n根(c（c）_ttpo2型 - b条**2 - 年**,2)
b条 一 =2 * x个 * 年秒[0]
如果s[1]==真并且gcd（a，b，c）==1:
A375750型.append（k）
c（c） = x个**2 + 年**2 打破
如果 最小值(一,b条) % 2 = -=1 和 最大值(一,b条) % 2 == 0:
如果gcd（a，b，c）==1且（c-1）%4==0：
A375750型.添加（（c-1）//4）
x个 k个 += 1
打印(已排序(列表(A375750型))[:60]) #卡尔·海因茨·霍夫曼2024年9月17日
