%I#14 2016年2月19日06:40:39
%S 1,1,2,2,4,5,9,11,21,28,50,681231733104417891147204429995351,
%电话:793814143211383768656729101144153085273077415407741301,
%电话:113237320218313100128553778285851907615225373234914134200737864979069116239502
%N根树的森林数量,其中花费N美元,一棵N节点树花费2n-1美元。
%H Alois P.Heinz,n的表格,n=1..1000的a(n)</a>
%H N.J.A.Sloane,转换</a>
%{c(x)}的F欧拉变换,其中{c(x)}具有g.F.B(x^2)/x=x+x^3+2*x^5+4*x^7+9*x^9+。。。和B(x)=x+x^2+2*x^3+4*x^4+9*x^5+。。。是A000081的g.f。
%p带有(numtheory):
%p b:=proc(n)选项记忆;局部d,j`如果`(n<=1,n,
%p(加(加(d*b(d),d=除数(j))*b(n-j),j=1..n-1))/(n-1)
%p端:
%pc:=proc(n)局部r`如果`(irem(n,2,'r')=0,0,b(r+1))结束:
%p a:=proc(n)选项记忆`如果`(n=0,1,
%p加(加(d*c(d),d=除数(j))*a(n-j),j=1..n)/n)
%p端:
%p序列(a(n),n=1..50);#_Alois P.Heinz_,2013年5月16日
%tb[n_]:=b[n]=如果[n<=1,n,(和[Sum[d*b[d],{d,除数[j]}]*b[n-j],{j,1,n-1}])/(n-1)];c[n_]:=如果[Mod[n,2]==0,0,b[n/2//天花板]];a[n_]:=a[n]=如果[n==0,1,Sum[Sum[d*c[d],{d,Divisors[j]}]*a[n-j],{j,1,n}]/n];表[a[n],{n,1,50}](*_Jean-François Alcover_,2016年2月19日,在_Alois P.Heinz_*之后)
%K nonn,简单
%氧1,3
%基督教G.鲍尔_
|