[seqfan] Re: Number of Trees

Richard Mathar mathar at strw.leidenuniv.nl
Fri Mar 6 23:53:37 CET 2009


I've placed a bfile into http://www.strw.leidenuniv.nl/~mathar/progs/b000055.txt
The Maple program is:

# maple: create b-file b000055.txt
# Richard J. Mathar, 2009-03-06
A000081 := proc(n)
	option remember ;
	local d, j;
	if n<=1 then
		n
	else
		add ( add(d*procname(d),d=numtheory[divisors](j)) *procname(n-j), j=1..n-1)/ (n-1) ;
	fi ;
end:

  	

A000055 := proc(nmax)
	local a81,n,t,a,j ;
	a81 := [seq(A000081(i),i=0..nmax)] ;
	a := [] ;
	for n from 0 to nmax do
		if n = 0 then
			t := 1+op(n+1,a81) ;
		else
			t := op(n+1,a81) ;
		fi;
		if type(n,even) then
			t := t-op(1+n/2,a81)^2/2 ;
			t := t+op(1+n/2,a81)/2 ;
		fi;
		for j from 0 to (n-1)/2 do
			t := t-op(j+1,a81)*op(n-j+1,a81) ;
		od:
		a := [op(a),t] ;
	od:
	a ;
end:

# maximum b-file elements: 1000
L := A000055(1000) ;
read("transforms3") ;
LISTTOBFILE("b000055.txt",L,0) ;




More information about the SeqFan mailing list