#OEIS A368207的Python程序#Don Knuth翻译的CWEB程序#Michael S.Branicky，2023年12月19日最大值=10001th=[0代表范围内的_（最大值）]xs=[0表示_在范围内（最大值）]总和=0bfile=打开（“b368207.txt”，“w”）定义th_and_xs（）：全局th，xsm=1nn=3对于范围（1，maxn）内的n：d=米当n%d:d-=1时th[n]=dk=0对于范围（m，0，-1）中的d：如果不是n%d:如果d*d==n，则k+=d+d-1，否则d+d+d+2xs[n]=k；如果n==nn：m+=1nn+=m+m+1定义b（n）：全局th，xs，sum总和=0k、 kk=1，n-1而k<kk：d、 dd=k//第[k]条，第[kk]条如果d<dd：对于范围（d，dd）内的i：如果不是k%i：对于范围（dd，i，-1）中的j：如果不是kk%j：溶液（k//i，i，j，kk//j）k+=1kk-=1bfile.write（f“{n}{sum+xs[n]}\n”）def溶液（a、b、c、d）：全局总和如果a==b：如果c==d：总和+=1其他：总和+=2elif c==d：总和+=2其他：总和+=4bfile.write（“#computed by a Python translation of Don Knuth’s CWEB program on 19 December 2023”）；th_and_xs（）对于范围（1，maxn）中的n：b（n）bfile.write（“\n”）bfile.close（）