（PARI）lista（nmax）={my（list=list（），rmax=0，e2，e3，r）；对于（n=1，nmax，e2=evaluation（n，2）；e3=evaluance（n，3）；r=if（2^e2*3^e3==n，（2^（e2+1）-1）*（3^（e3+1）-1）/（2*n），0）；如果（r>rmax，rmax=r；listput（list，n））））；Vec（list）}；
（Python）
从分数导入分数
从sympy导入多重性为v
来自itertools导入计数，takewhile
def f（n）：返回分数（（2**（v（2，n）+1）-1）*（3**（v（3，n）+1）-1）//2，n
定义平滑3（lim）：
pows2=列表（takewhile（λx：x<lim，（计数（0）中i的2**i））
pows3=列表（takewhile（λx：x<lim，（计数（0）中的i为3**i））
返回已排序的（如果c*d<=lim，则返回pows2中c的c*d，返回pows3中d的c）
定义缺陷（lim）：
数据，记录，记录=平滑3（lim），[]，-1
对于zip中的argv、v（数据、地图（f、数据））：
如果v>记录：记录=v；记录.附加（argv）
退货记录
