|
| |
| |
|
|
|
1, 4, 12, 47, 170, 750, 3255, 16010, 81199, 448156, 2579626, 15913058, 102488024, 698976419, 4976098729, 37195337408, 289517846210, 2352125666883, 19841666995265, 173888579505200, 1577888354510786, 14820132616197925, 143746389756336173, 1438846957477988926
(列表;图表;参考文献;听;历史;文本;内部格式)
|
|
|
|
抵消
|
1、2
|
|
|
评论
|
用多重数划分n个多集的方法n的某些划分。
大小为n的强正规多集的多集划分数，其中，如果有限多集覆盖具有弱递减重数的正整数的初始区间，则它是强正规的。（弱）正常版本为2005年2月. -古斯·怀斯曼2019年12月31日
|
|
|
链接
|
|
|
|
例子
|
a（3）=12，因为A025487号有3个因素，即8、12、30；f（8）=3，f（12）=4，f（30）=5和3+4+5=12。
强正规多集的a（1）=1到a（3）=12个多集划分：
{{1}} {{1,1}} {{1,1,1}}
{{1,2}} {{1,1,2}}
{{1},{1}} {{1,2,3}}
{{1},{2}} {{1},{1,1}}
{{1},{1,2}}
{{1},{2,3}}
{{2},{1,1}}
{{2},{1,3}}
{{3},{1,2}}
{{1},{1},{1}}
{{1},{1},{2}}
{{1},{2},{3}}
（结束）
|
|
|
MAPLE公司
|
带有（数字理论）：
g： =proc（n，k）选项记忆；
`if`（n>k，0，1）+`if`（isprime（n），0，
加法（`if`（d>k，0，g（n/d，d）），d=除数（n）减去{1，n}）
结束时间：
b： =程序（n，i，l）
`如果`（n=0，g（mul（ithprime（t）^l[t]，t=1..nops（l））$2），
`如果`（i<1，0，加（b（n-i*j，i-1，[l[]，i$j]），j=0..n/i））
结束时间：
a： =n->b（n$2，[]）：
|
|
|
数学
|
g[n_，k_]：=g[n，k]=如果[n>k，0，1]+如果[PrimeQ[n]，0，Sum[If[d>k，0,g[n/d，d]]，{d，除数[n]~补~{1，n}}]]；b[n_，i_，l]：=如果[n==0，g[p=乘积[Prime[t]^l[[t]]，{t，1，Length[l]}]，p]，如果[i<1，0，Sum[b[n-i*j，i-1，Join[l，Array[i&，j]]]，}，0，n/i}]]；a[n_]：=b[n，n，{}]；表格[打印[an=a[n]]；an，{n，1，13}](*Jean-François Alcover公司2013年12月12日之后阿洛伊斯·海因茨*)
|
|
|
黄体脂酮素
|
（Python）
从sympy.core.cache导入缓存
从辛导入除数、isprime、prime
从运算符导入mul
@缓存
定义g（n，k）：
return（0如果n>k else 1）+（0如果isprime（n）else sum（g（n//d，d）对于除数（n）中的d[1:-1]如果d<=k））
@缓存
定义b（n，i，l）：
如果n==0：
p=减少（mul，（素数（t+1）**l[t]表示范围内的t（len（l）））
返回g（p，p）
其他：
如果i＜1，则返回0，否则求和（[b（n-i*j，i-1，l+[i]*j）对于范围（n//i+1）]中的j）
定义a（n）：
返回b（n，n，[]）
对于范围（1，11）中的n：打印（a（n））#因德拉尼尔·戈什2017年8月19日，在Maple代码之后
（PARI）
EulerT（v）={Vec（exp（x*Ser（dirmul（v，vector（#v，n，1/n）））-1，-#v）}
D（p，n）={my（v=向量（n））；对于（i=1，#p，v[p[i]]++）；my（u=EulerT（v））；Vec（1/prod（k=1，n，1-u[k]*x^k+O（x*x^n））-1，-n）/prod（i=1，#v，i^v[i]*v[i]！）}
序列（n）={my（s=0）；对于部分（p=n，s+=D（p，n））；s}\\安德鲁·霍罗伊德2020年12月30日
|
|
|
交叉参考
|
|
|
|
关键词
|
非n,美好的
|
|
|
作者
|
|
|
|
扩展
|
|
|
|
状态
|
经核准的
|
| |
|
|
