%I#20 2021年4月21日11:25:22
%S 1,1,2,1,2,2,6,2,1,20,8,2,60,2,10450,1,2726,21403321,14,25896,
%电话572,1626426264,2394406,2,1226020,2051886961584,2,222044895,
%电话:38740,220959503,2676478164163,26684908631206529681847044980,2127370696620376153361803937344
%N N到N的奇数除数的组成数(有序分区)。
%H Alois P.Heinz,n的表格,n=0..2000的a(n)</a>
%H<a href=“/index/Com#comp”>与合成相关的序列的索引条目</a>
%F a(n)=[x^n]1/(1-和{d|n,d正奇数}x^d)。
%如果n是2的幂,则F a(n)=1。
%如果n是奇数素数,则F a(n)=2。
%e a(10)=8,因为10有4个除数{1,2,5,10},其中2个是奇{1,5},所以我们有[5,5],[5,1,1,1],[1,5,1。
%p(数字理论):
%p a:=proc(n)选项记忆;局部b,l;
%pl,b:=select(x->is(x::奇数),除数(n)),
%p过程(m)选项记忆`如果`(m=0,1,
%p加(`if`(j>m,0,b(m-j)),j=l)
%p端;b(n)
%p端:
%p序列(a(n),n=0..60);#_Alois P.Heinz,2017年3月30日
%t表[d=除数[n];系数[级数[1/(1-和[Boole[Mod[d[[k]],2]==1]x^d[[k]],{k,长度[d]}]),{x,0,n}],x,n],{n,0,57}]
%o(Python)
%o来自sympy导入除数
%o从sympy.core.cache导入缓存
%o@缓存
%o定义a(n):
%o l=[x代表除数(n)中的x,如果x%2]
%o@缓存
%o def b(m):如果m==0,则返回1,否则求和(如果j<=m,则b(m-j)用于l中的j)
%o返回b(n)
%o打印([a(n)代表范围(61)中的n)]#_Indranil Ghosh,2017年8月1日,在Maple代码之后
%Y参见A000045、A005408、A032021、A100346、A171565。
%K nonn公司
%0、4
%A _Ilya Gutkovskiy_,2017年3月27日
|