%I#18 2020年2月14日16:37:27
%S 1,2,5,10,11,17,22,23,31,34,41,46,47,55,59,62,67,73,82,83,85,94,97,
%电话:103109110115118127134137146149155157166167170179187,
%电话:191194197205206211218227230235235241253254257269277283295308307313143331347
%N将海因茨数按递增顺序划分为不同的奇数部分。
%C我们将分区的Heinz数p=[p_1,p_2,…,p_r]定义为乘积(p_j-th素数,j=1…r)(_Alois p.Heinz_在A215366中使用的概念是分区的“编码”)。例如,分区[1,1,2,4,10]的Heinz数是2*2*3*7*29=2436。
%C在Maple程序中,子程序B生成Heinz数为n的分区。
%C如果将Maple程序中的350替换为更大的数字,则会获得更多的项。
%D G.E.Andrews,《分割理论》,Addison-Wesley,Reading,马萨诸塞州,1976年。
%D G.E.Andrews,K.Eriksson,《整数分区》,剑桥大学出版社,2004年,剑桥。
%H Alois P.Heinz,n的表格,n=1..10000的a(n)</a>
%e 170位于序列中,因为它是分区[1,3,7]的Heinz数;实际上,(第一素数)*(第三素数)x(第七素数)=2*5*17=170。
%p with(numtheory):B:=proc(n)local pf:pf:=op(2,ifactors(n)):[seq(seq(pi(op(1,op(i,pf={true})然后DO:=`union`(DO,{q})else end if end DO:DO;
%p#第二个Maple程序:
%p a:=proc(n)选项记忆;局部k;
%p表示从1+`if`(n=1,0,a(n-1))do中的k
%p如果在映射中不为假(i->i[2]=1和数字理论
%p[pi](i[1])::奇数,ifactors(k)[2])然后中断fi
%p od;k个
%p端:
%p序列(a(n),n=1..100);#_Alois P.Heinz,2016年5月10日
%t a[n_]:=a[n]=模[{k},对于[k=1+如果[n==1,0,a[n-1]],真,k++,如果[AllTrue[FactorInteger[k],#[2]]==1&&OddQ[PrimePi[#[1]]]&],中断[]]];k] ;加入[{1},数组[a,100]](*_Jean-François Alcover_,2016年12月10日,继_Alois P.Heinz_*之后)
%Y参见A215366、A258117。
%K nonn公司
%O 1,2号机组
%德国电子报,2015年5月20日
%E a(1)=1由_Alois P.Heinz插入,2016年5月10日
|