%I#15 2024年4月9日08:39:08

%S 1,6,20,72,72,72，20,72，7217280480017280,7217280115200051845184，

%电话：5184960005184345600124416032000012441608294400089579520，

%电话：71663616000298598400124416082944000230400008294400019906560000

%N从2^N开始重复还原usigma的终点。

%C设PF_p（n）是p除以n的最大幂。例如，在p=2、3和5的情况下，PF_2（n）=A006519（n）、PF_3（n=A038500（n）和PF_5（n）=5^A112765（n）。

%C多诱导PF_（p1，p2，…）。。。

%对于每个n，我们定义了一个辅助序列b（k），从b（0）=2^n开始，通过b（k+1）=A034448（b（k。b（k）在某个k处终止，b（k）=1。此外，还有一个辅助并行序列c（k），由c（0）=2^n定义，递归地定义为c（k+1）=c（k。

%序列由a（n）＝C（k）定义，辅助序列C在b终止的点处。

%C序列a（n）的所有值都是5-光滑的，即A051037的成员。

%e n=5

%e b（n）：2^5->11->1

%e c（n）：2^5->2^5*3->2^3*3^2

%因此a（5）=c（2）=2^3*3^2=72。

%p PF:=proc（n，p）局部nshf，a；a:=1；nshf：=n；而（nshf mod p）=0执行nshf：=nshf/p；a:=a*p；od:a；结束时间：

%p A146891:=进程（n）局部b，a，k，t；

%磅：=[2^n]；

%p，而op（-1，b）<>1 do

%p t：=A034448（op（-1，b））；

%p b：=[操作（b），t/A006519（t）/A038500（t）/PF（t，5）]；

%日期：

%p a：=2^n；

%p代表k从2到nops（b）do

%p t：=A034448（op（k-1，b））；

%p a：=a/A006519（t）*A038500（t）*PF（t，5）；

%日期：

%p a；

%p端：

%p#_R.J.Mathar_，2009年6月24日

%t PF[n_，p_]：=p^整数指数[n，p]；

%t usigma[n_]：=如果[n==1，1，Times@@（1+Power@@@FactorInteger[n]）]；

%t A146891[n_]：=模块[{b，a，k，t}，

%tb={2^n}；

%t同时[b[[-1]]！=1,

%t t=usigma[b[[-1]]]；

%t b=附加[b，t/PF[t，2]/PF[t，3]/PF[t,5]]；

%t a=2^n；

%t对于[k=2，k<=长度[b]，k++，

%t t=usigma[b[[k-1]]；

%t a=a/PF[t，2]*PF[t，3]*PF[t，5]]；

%t a]；

%t表[A146891[n]，{n，0，32}]（*_Jean-François Alcover_，2024年4月9日，在R.J.Mathar_*之后）

%Y参考A146892、A151659。

%K nonn公司

%0、2

%A_Yasutoshi Kohmoto，2009年4月17日

%E更多条款摘自R.J.Mathar_，2009年6月24日

%E编辑：R.J.Mathar_，2009年7月2日

%E由R.J.Mathar修正的a（n）和c（k）之间关系的描述，2009年7月7日