A080193-OEIS公司

A080193号

5-平滑数，但不是3-平滑数。

5, 10, 15, 20, 25, 30, 40, 45, 50, 60, 75, 80, 90, 100, 120, 125, 135, 150, 160, 180, 200, 225, 240, 250, 270, 300, 320, 360, 375, 400, 405, 450, 480, 500, 540, 600, 625, 640, 675, 720, 750, 800, 810, 900, 960, 1000, 1080, 1125, 1200, 1215, 1250, 1280, 1350

(列表;图表;参考;听;历史;文本;内部格式)

抵消

1,1

形式为2^r*3^s*5^t且r，s>=0，t>0的数字。

也就是说，5平滑数是5的倍数-查尔斯·格里特豪斯四世2015年3月19日

链接

阿米拉姆·埃尔达尔，n=1..10000时的n，a（n）表

配方奶粉

发件人阿米拉姆·埃尔达尔2020年11月10日：（开始）

a（n）=5*A051037号（n） ●●●●。

Sum_｛n＞＝1｝1/a（n）＝3/4。（结束）

例子

15=3*5是一个术语，但18=2*3^2不是。

数学

选择[Range[1000]，FactorInteger[#][[-1，1]]==5&](*阿米拉姆·埃尔达尔2020年11月10日*）

黄体脂酮素

（PARI）{m=1440；z=[]；对于（r=0，floor（log（m）/log（2）），a=2^r；对于（s=0，loor（log（m/a）/log（3）），b=a*3^s；对于（t=1，floor，log（m/b）/log（5）），z=concat（z，b*5^t））））；z=vecsort（z）；对于（i=1，长度（z），打印1（z[i]，“，”）}

（PARI）列表（lim）=我的（v=列表（），x=1，y，z）；而（x*=5）<=lim，y=x/3；而（y*=3）<=lim，z=y/2；而（（z*=2）<=lim，listput（v，z）））；集合（v）\\查尔斯·格里特豪斯四世2015年3月19日

（Python）

从sympy导入integer_log

定义A080193号（n）：

定义平分（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）：

c=n+x

对于范围内的i（integer_log（x，5）[0]+1）：

对于范围内的j（integer_log（y:=x//5**i，3）[0]+1）：

c-=（y//3**j）.bit_length（）

返回c

返回二分（f，n，n）*5#柴华武2024年9月16日

（Python）#序列的初始段速度更快

导入heapq

从itertools导入islice

def A080193gen（）：#术语生成器

v、 oldv，h，psmooth_primes，=1,0，[1]，[2,3,5]

为True时：

v=heapq.heappop（h）

如果v！=旧版本：

产量5*v

oldv=v

对于psmooth_primes中的p：

heapq.heappush（h，v*p）

打印（列表（islice（A080193gen（），55））#迈克尔·布拉尼基2024年9月18日

交叉参考

囊性纤维变性。A051037号,A003586号.

上下文中的序列：A313735型 A067761号 A358848飞机*A313736型 A313737型 A313738型

相邻序列：A080190型 A080191号 A080192号*A080194号 A080195号 A080196号

关键词

容易的,非n,改变

作者

克劳斯·布罗克豪斯2003年2月10日

状态

经核准的