（PARI）我的（v=列表（10^15），rec=0）；打印1（0，“，”）；对于（n=1，#v，如果（numdiv（v[n]）>rec，rec=numdiv；打印1（v[n]，“，”））\\宋嘉宁，2021年5月20日，见计划A054994号
（Python）
从数学导入prod
从sympy导入isprime
素数一致_1_mod_4=[5]
def素数4k_plus1（i）：#与1模4同余的第i素数
当i>=len（primes_congruent_1_mod_4）：#按需生成素数
n=素数一致_1_mod_4[-1]+4
当不为素数（n）时：n+=4
素数一致_1_mod_4.追加（n）
返回素数congruent_1_mod_4[i]
TO_DO={（1，（））}
为True时：
半径，指数=最小值（TO_DO）
屈服半径，指数
TO_DO.remove（（半径，指数））
TO_DO.update（后续项（半径、指数））
def后继项（r，指数）：
对于枚举中的i，e（指数）：
如果i==0或指数[i-1]>e：
产量（r*prime4kplus1（i），指数[：i]+（e+1，）+指数[i+1:]）
如果指数==（）或指数[-1]>0：
产量（r*prime4kplus1（len（指数）），指数+（1，））
n、 记录，半径=1,1,0
num_points=4*prod（（e+1）代表expo中的e）
如果num_points>记录：
记录=num_points
n+=1
如果n==27：断裂
打印（）
