%I#10 2023年10月27日19:16:26

%S 1,1,1,1,2,3,6,8,18,42,6420043254696038886000151802995238976，

%电话：67200285600393984163296031752004165392610560038413440，

%电话：55339200114048000205632000028021950448156800262144512037777259521294084152

%N除以1的整数k^6的数目*2!*3!*...*不！。

%H Alois P.Heinz，n表，n=1..1000时的a（n）

%e a（7）对这些整数k^6进行计数，将125411328000=A000178（6）除以：1，64，729，4096，46656，2985984，这些是k^6，表示k=1，2，3，4，6，12。

%p b:=proc（n）选项记忆；添加（i[2]*x^numtheory[pi]（i[1]），

%p i=ifactors（n）[2]）+`if`（n=1，0，b（n-1））

%p结束：

%pc:=proc（n）选项记忆；b（n）+`if`（n=1,0，c（n-1））结束：

%p a:=n->（p->mul（iquo（系数（p，x，i），6）+1，i=1.度（p））（c（n））：

%p序列（a（n），n=1..30）；#_Alois P.Heinz，2014年10月16日

%tz=40；p[n_]：=乘积[k！，{k，1，n}]；

%t f[n_]：=f[n]=系数整数[p[n]]；

%t r[m_，x_]：=r[m，x]=m*楼层[x/m]

%tu[n_]：=表[f[n][[i，1]]，{i，1，长度[f[n]}]；

%tv[n_]：=表[f[n][[i，2]]，{i，1，长度[f[n]}]；

%t t[m_，n_]：=应用[次数，1+r[m，v[n]]/m]

%t m=6；表[t[m，n]，{n，1，z}]（*A248824*）

%Y参见A000178、A248784、A248821、A24882、A248833。

%K nonn，简单

%O 1,5型

%A_Clark Kimberling_，2014年10月15日