"""███████████████ 作者、日期、时间。。。██████████████████████████████████████████████████████卡尔·海因茨·霍夫曼мм斯佩耶斯特拉ße 9м68649格罗ß-罗赫海姆мм德国мм手机+49 175 2037149（是的，我正在使用WhatsApp）邮箱：@kallimusic.de ██07.12.2021 ██10:50甲基乙酰胺██灵感来自Jon E.Schoenfield的算法。（请参见https://oeis.org/A349324) ███████████████████████████████████████████████████████████████████████████████████████████████ """导入时间从gmpy2导入issquare从sympy导入除数#██████████████ 主参数。（不要触摸其他任何东西。；-）███████████████████████████#██howfar=643000000#643000000包含至少5个解决方案。██#██#██████████████████████████████████████████████████████████████████████████████████████████ t1=time.time（）#开始计算-------------------------------------------------solutionlist=[0]#需要这个零作为排序和拆分的净空对于范围（1，int（howfar**0.5）+1）内的y：d=除数（y**4）#d是y^4的所有除数的列表对于范围（0，len（d）//2+1）：如果（d[it]+d[-（it+1）]）%2==0：a_term=（d[it]+d[-（it+1）]）//2如果a _term<=多远且不是平方（a _tern）：解决方案列表.append（a_term）打印（“计算时间=%3.3f”%（time.time（）-t1））solutionlist=排序（solutionlist）#排序和拆分------------------------------solutionlist.append（0）#这个零是排序和拆分所需的净空。计数器=0对于范围（1，len（solutionlist）-1）中的z：与计算相比性能较差。如果解决方案列表[z-1]！=解决方案列表[z]和解决方案列表[z+1]！=解决方案列表[z]：f=打开（“1_Solutions.txt”，“a”）f.write（str（解决方案列表[z]）+“\n”）f.关闭（）elif解决方案列表[z-1]！=solutionlist[z]和solutionlist[z+1]==解决方案列表[z]：计数器+=1elif solutionlist[z-1]==解决方案列表[z]和解决方案列表[z+1]==方案列表[z]：计数器+=1else:#最后出现在多个术语的末尾。计数器+=1f=打开（str（计数器）+“_Solutions.txt”，“a”）f.write（str（解决方案列表[z]）+“\n”）f.关闭（）计数器=0打印（“计算时间加上排序和拆分时间=%3.3f”%（time.time（）-t1））#如果有人对“分裂”有更好的想法。。。。。。。。。给我发邮件。（邮件地址见上文）