%I#22 2020年11月23日07:55:51
%S 1,3,7,12,28,31,60,91124168127360403546508744117016512418,
%电话288028213048204743683751595293609906818812493819119344,
%电话:153672256822506243842880026113931232764511815952049128792489994492202112320116281106483160797
%N最小数的除数之和(A007416)。
%C具有与A000203及其所有子序列相同的对称表示。
%H Amiram Eldar,n的表,n的a(n)=1..10000(根据A007416的b文件计算)
%F a(n)=A000203(A007416(n))。
%t(*列表minDiv[n,b]中的第d个元素是最小的数字k<=n,除数正好是d<=b,否则为零。只要不等式失败,计算就会停止。*)
%t最小刻度[n_,b_]:=
%t模块[{list=Array[0&,b],k=1,d},
%t而[k<=n,d=除数Sigma[0,k];
%t如果[d<=b&&list[[d]]==0,list[[d]]=k];
%t如果[d<=b,k++,k=n+2]];列表]
%t a256259[n,b]:=
%t映射[DivisorSigma[1,#]&,排序[Select[minDiv[n,b],#!=0 &]]]
%t a256259[100000,300](*计算序列的前60个元素*)
%t(*哈特穆特·F·W·霍夫特,2015年4月27日*)
%Y参见A000203、A005179、A007416、A007626、A053212、A063072、A081533、A237271、A237593。
%K nonn公司
%O 1,2号机组
%2015年4月20日,A_Omar E.Pol_
|