登录
A057568号
n的分区数,其中n除以各部分的乘积。
30
1, 1, 1, 2, 1, 2, 1, 6, 5, 5, 1, 22, 1, 11, 23, 80, 1, 113, 1, 150, 85, 45, 1, 737, 226, 84, 809, 726, 1, 1787, 1, 4261, 735, 260, 1925, 9567, 1, 437, 1877, 16402, 1, 14630, 1, 9861, 33057, 1152, 1, 102082, 19393, 57330, 10159, 30706, 1, 207706, 47927, 200652
抵消
1,4
链接
阿洛伊斯·海因茨,n=1..1000时的n,a(n)表(Antti Karttunen提供的术语n=1..73)
例子
发件人古斯·怀斯曼2019年7月4日:(开始)
a(1)=1到a(9)=5分区如下。这些分区的Heinz数由下式给出A326149型.
(1) (2) (3) (4) (5) (6) (7) (8) (9)
(22) (321) (44) (63)
(422) (333)
(2222) (3321)
(4211) (33111)
(22211)
(结束)
MAPLE公司
b: =proc(n,i,t)选项记忆`如果`(n=0,
`如果`(t=1,1,0),`如果`(i<1,0,b(n,i-1,t)+
`如果`(i>n,0,b(n-i,min(i,n-i),t/igcd(i,t)))
结束时间:
a: =n->`if`(i素数(n),1,b(n$3)):
seq(a(n),n=1..70)#阿洛伊斯·海因茨2017年12月20日
数学
表[Length[Select[Integer Partitions[n],Divisible[Times@@#,n]&]],{n,20}](*古斯·怀斯曼2019年7月4日*)
b[n_,i_,t_]:=b[n,i,t]=如果[n==0,如果[t==1,1,0],如果[i<1,0,b[n、i-1,t]+如果[i>n,0,b[n-i,Min[i,n-i],t/GCD[i,t]]];
a[n_]:=如果[PrimeQ[n],1,b[n,n,n]];
数组[a,70](*Jean-François Alcover公司2021年5月21日之后阿洛伊斯·海因茨*)
黄体脂酮素
(方案)
;; 这是一个简单的算法,它扫描每个n的所有分区。有关fold_over_partitions_of,请参阅A000793号.
(定义(A057568美元n) (let((z(list 0)))(n 1的fold_over_partitions_of n 1*(λ(partprod))(如果(零?(模数partprodn))(set-car!z(+1(car z))))
;;Antti Karttunen公司2017年12月20日
关键词
非n
作者
勒罗伊·奎特2000年10月4日
扩展
更多术语来自詹姆斯·塞勒斯2000年10月9日
状态
经核准的