登录
A258116型
海因茨数按分区的递增顺序划分为不同的奇数部分。
22
1, 2, 5, 10, 11, 17, 22, 23, 31, 34, 41, 46, 47, 55, 59, 62, 67, 73, 82, 83, 85, 94, 97, 103, 109, 110, 115, 118, 127, 134, 137, 146, 149, 155, 157, 166, 167, 170, 179, 187, 191, 194, 197, 205, 206, 211, 218, 227, 230, 233, 235, 241, 253, 254, 257, 269, 274, 277, 283, 295, 298, 307, 310, 313, 314, 331, 334, 335, 341, 347
抵消
1, 2
评论
我们将分区p=[p_1,p_2,…,p_r]的Heinz数定义为乘积(p_j-th素数,j=1…r)(阿洛伊斯·海因茨在里面A215366型作为分区的“编码”)。例如,分区[1,1,2,4,10]的Heinz数是2*2*3*7*29=2436。
在Maple程序中,子程序B生成Heinz编号为n的分区。
如果将Maple程序中的350替换为更大的数字,则会获得更多术语。
参考文献
G.E.Andrews,《分割理论》,Addison-Wesley,Reading,马萨诸塞州,1976年。
G.E.Andrews、K.Eriksson,《整数分区》,剑桥大学出版社,2004年,剑桥。
链接
阿洛伊斯·海因茨,n=1..10000时的n,a(n)表
例子
170位于序列中,因为它是分区[1,3,7]的Heinz数;实际上,(第一素数)*(第三素数)x(第七素数)=2*5*17=170。
MAPLE公司
with(numtheory):B:=proc(n)local pf:pf:=op(2,ifactors(n)):[seq(seq(pi(op(1,op(i,pf))),j=1..op(2、op(i、pf){true})then DO:=`union`(DO,{q})else end if end DO:DO;
#第二个Maple项目:
a: =proc(n)选项记忆;局部k;
对于1+`if`(n=1,0,a(n-1))中的k,do
如果在map中不是false(i->i[2]=1和numtheory
[pi](i[1])::奇数,ifactors(k)[2])然后中断fi
od;k
结束时间:
seq(a(n),n=1..100); #阿洛伊斯·海因茨2016年5月10日
数学
a[n_]:=a[n]=模[{k},对于[k=1+如果[n==1,0,a[n-1]],True,k++,如果[AllTrue[FactorInteger[k],#[2]]==1&&OddQ[PrimePi[#[1]]]&],Break[]]];k];联接[{1},数组[a,100]](*Jean-François Alcover公司,2016年12月10日之后阿洛伊斯·海因茨*)
关键词
非n
作者
Emeric Deutsch公司2015年5月20日
扩展
a(1)=1由插入阿洛伊斯·海因茨,2016年5月10日
状态
经核准的