g: =proc(n)选项记住`如果`(n=0,1,加上(g(n-j)*add(
`如果`(d::奇数,d,0),d=numtheory[除数](j)),j=1..n)/n)
结束时间:
s: =proc(n)选项记忆;展开(`if`(n=0,1,
x*加法(s(n-j)*二项式(n-1,j-1),j=1..n))
结束时间:
S: =proc(n,k)选项记忆;系数(s(n),x,k)结束:
b: =proc(n,i)选项记忆`如果`(n=0,1,
`如果`(i=0,g(n),加上(b(n-j,i-1),j=0..n))
结束时间:
a: =n->加(S(n,j)*b(n,j),j=0..n):
seq(a(n),n=0..20);