%I#20 2018年2月27日20:01:40
%S 1,1,2,1,4,4,8,1,6,12,16,6,32,28,1,64,16128,24,96,80256,8,44192,
%电话22,80512,961024,1288448224,3020481024800,4040964008192,
%电话:240168230416384.10360204211267232768,6883216053765120
%形状为Heinz数为N的整数分区的正规半标准Young表的N个数。
%如果表的条目跨越正整数的初始区间,则表是正常的。整数分区的Heinz数(y_1,…,y_k)是素数(y_1)**质数(yk)。
%D Richard P.Stanley,《枚举组合数学》第2卷,剑桥大学出版社,1999年,第7.10章。
%H FindStat-组合统计查找器,<a href=“http://www.findstat.org/SemistandardTableaux“>半标准杨表</a>
%设b(n)=Sum_{d|n,d>1}b(n*d'/d),其中如果d=Product_i素数(s_i)^m(i),则d'=Product_i-素数(s1)^m。则a(n)=b(conj(n)),其中conj=A122111。
%e a(9)=6表aux:
%e 1 3 1 2 1 2 1 2 1 1 1 1 1 1
%e 2 4 3 4 3 3 2 2 2 3 2
%t conf[y_List]:=如果[Length[y]==0,y,表[Length[Select[y,#>=k&]],{k,1,Max[y]}];
%t conf[n_Integer]:=Times@@Prime/@conf[If[n===1,{},Join@@Cases[FactorInteger[n]//反向,{p_,k_}:>表[PrimePi[p],{k}]];
%t ssyt[n_]:=如果[n===1,1,总和[ssyt[n/q*Times@@Cases[FactorInteger[q],{p_,k_}:>如果[p===2,1,NextPrime[p,-1]^k]],{q,Rest[Divisors[n]]];
%t表[ssyt[conf[n]],{n,50}]
%Y参见A000085、A001222、A056239、A063834、A112798、A122111、A138178、A153452、A191714、A210391、A228125、A296150、A296560、A29656、A299202、A299966、A300056、A300121。
%K nonn公司
%氧1,3
%A _Gus Wiseman_,2018年2月14日