(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日