从数学导入日志定义IsPrime（n）：如果n<2：返回0elif n==2或n==3：返回1elif n%2==0或n%3==0：返回0其他：d、 dd=5，2而d*d<=n和n%d！=0:d、 dd=d+dd，6-dd如果d*d<=n：返回0其他：返回1定义nextp（p）：m=p%6如果m==1：d=4其他：#m==5d=2ptry，d=p+d，6-d而不是IsPrime（ptry）：ptry，d=ptry+d，6-d返回ptry文件=打开（“A166113.txt”，“r”）n=0打印（“n:A256891 A343194 A343195‘直径’A270305‘标准化平均间隙’”）而1：n=n+1p5=int（file.readline（））p6=下一个p（p5）p7=下一个p（p6）p8=下一个p（p7）p9=下一个p（p8）如果p6-p5==p9-p8：#类型1b、 c=p6-p5，p7-p6st=“类型1”elif p6-p5！=p9-p8:#类型2b、 c=p7-p5，p6-p7st=“类型2”其他：打印（“错误”）a、 d，s=p5-3*b-c，6*b+2*c，3*p5p1，p2，p3，p4，p5，p6，p7，p8，p9=a，a+b，a+2*b，a=2*b+c，a+3*b+c，a+4*b+c+4*b+2*c，a+5*b+2*c，a+6*b+2*c打印（n，“：”，a，b，c，d，s，d/（8*log（p1））