发件人:亚历克·米哈伊洛夫斯日期:2006年12月17日星期日21:01:22-0600阿奇姆·弗拉门坎普的网页http://wwwhomes.uni-bielefeld.de/achim/addition_chain.html包含指向文件add25.bits.gz(约3.34MB)的链接2比特分组数0≤a(n)≤3,这样A003313(n)=地板(log_2(n))+天花板(log_2[nu(n)])+a(n)其中nu(n)是以2为基数的n的位数之和;对于从1到的n2^25=33554432,由阿希姆·弗拉门坎普和尼尔·克利夫。当我想检查n=27782的计算结果时在修改ddt.c(上述网页上的链接)时,我使用了Python(在交互模式下),今天我把它改写为Python脚本,#-----------查找.py----------从sys导入argv从结构导入解包从数学导入日志、地板、天花板f=打开(“./add25.bits”,'rb')n=整数(argv[1])-1第二周(n/4)r=拆包('B',f.read(1))[0]f.关闭()a=(r>>2*(n%4))%4n+=1a+=地板(对数(n,2))u=0当n>0时:u+=n%2n=n>>1打印int(a+ceil(log(u,2))#-------------电动势-------------它应该与未压缩的文件放在同一目录中add25.bits文件(8MB)。它可以在Windows中从cmd运行(在cd到该目录后,否则f中的路径应为以完整形式书写)作为python查找.py 27782例如,或使用从1到2^25的任何其他整数第27782页。亚历克·米哈伊洛夫斯http://mihailovs.com/Alec/