（PARI）测试（n）=（p=2,3，而（n%p==0，n/=p））；n==1；
对于（n=14000，如果（测试（n），打印1（n“，”））
（PARI）列表（lim）=我的（v=列表（），N）；对于（n=0，log（lim\1+.5）\log（3），n=3^n；而（N<=lim，listput（v，N））；N<<=1））；向量排序（Vec（v））\\查尔斯·格里特豪斯四世2011年6月28日
（PARI）列表（lim）=我的（v=列表（），N）；对于（n=0，logint（lim\=1,3），n=3^n；而（N<=lim，listput（v，N））；N<<=1））；集合（v）\\查尔斯·格里特豪斯四世2018年1月10日
（哈斯克尔）
导入数据。集合（Set、singleton、insert、deleteFindMin）
平滑：：设置整数->[Integer]
平滑s=x：平滑（插入（3*x）$插入（2*x）s’）
其中（x，s'）=删除查找最小值
a003586_list=平滑（单例1）
a003586 n=a003586_列表！！（n-1）
（鼠尾草）
定义为A003586（n）：
不返回任何（prime_divisors（n）中d的d！=2和d！=3）
@缓存函数
如果n=1：返回1
而不是A003586（k）：k+=1
返回k
（Python）
来自itertools导入计数，takewhile
定义缺陷（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）
（Python）
从sympy导入integer_log
定义平分（f，kmin=0，kmax=1）：
而f（kmax）>kmax:kmax<<=1
当kmax-kmin>1时：
kmid=kmax+kmin>>1
如果f（kmid）<=kmid：
kmax=kmid
其他：
kmin=kmid
返回kmax
定义f（x）：返回n+x-sum（（x//3**i）.bit_length（），用于范围内的i（integer_log（x，3）[0]+1））
返回二分（f，n，n）#柴华武2024年9月15日
（Python）#序列的初始段速度更快
导入heapq
从itertools导入islice
def A003586gen（）：#术语生成器
v、 oldv，h，psmooth_primes，=1，0，[1]，[2，3]
为True时：
v=heapq.heappop（h）
如果v！=旧版本：
产量v
oldv=v
对于psmooth_primes中的p：
堆堆（h，v*p）
打印（列表（islice（A003586gen（），65））#迈克尔·布拉尼基2024年9月17日
#包括<forward_list>
std:：forward_list<int>序列；
auto start_it=序列.before_begin（）；
for（int i=1；i<=阈值；i*=2）{
for（int inc=1；std:：next（start_it）！=sequence.end（）&&inc<=i；增加*=3）
++开始_ it；
auto-it=start_it；
对于（int j=1；i*j<=阈值；j*=3）{
sequence.template_after（it，i*j）；
for（int inc=1；std:：next（it）！=sequence.end（）&&inc<=i；inc*=2）
++它；
}
}
回归序列；
（Magma）[n:n在[1..4000]|素数除法（n）子集[2，3]]； //布鲁诺·贝塞利2012年9月24日
