登录
A370352型
存在n个或多个数k的最小数,其中k-a(n)=sopfr(k)+sopfr。
4
1, 1, 6, 22, 46, 526, 838, 838, 5667, 5667, 20158, 32127, 56697, 82617, 174718, 174718, 314492, 314492, 415789, 498957, 1142398, 1713598, 1713598, 2280067, 2280067, 4324316, 4324316, 5847653, 6918908, 6918908, 6918908, 9979197, 15855829, 24023995, 28274398, 28274398, 28274398, 28274398
抵消
1,3
链接
Michael S.Branicky,n=1..51时的n,a(n)表
例子
a(1)=a(2)=1,因为1是具有两个数字(k=1,6)的最小数字,使得1-1=0=sopfr(1)+sopfr。
a(3)=6,因为6是具有三个数字(k=20,21,26)的最小数字,使得20-6=14=sopfr(20)+sopfr。
黄体脂酮素
(Python)
来自sympy导入因子
从itertools导入计数,islice
从集合导入计数器
kcount,kmax=计数器(),0
def-sopfr(n):返回和(p的p*e,factorint(n).items()中的e)
定义f(n):
全球kcount,kmax
目标=n+sopfr(n)
对于射程内的k(kmax+1,2*目标+5):
k计数[k-sopfr(k)]+=1
kmax+=1
return kcount[目标]
def agen():#术语生成器
adict,n=dict(),1
对于计数(1)中的m:
v=f(米)
如果v不在adict:adict[v]=m中
对于范围(n,v+1)内的i:产量m;n+=1
打印(列表(islice(agen(),16))#迈克尔·布拉尼基2024年2月17日
关键词
非n
作者
斯科特·R·香农2024年2月16日
扩展
a(21)及以上迈克尔·布拉尼基2024年2月16日
状态
经核准的