[seqfan] Re: extending A003241

[Max Alekseyev]

The first disagreement happens at n=4 where you give 2 while A003241(4)=4.
Why not manually construct the trees counted in A003241 and see whose
value is correct?
Max

On Wed, Sep 28, 2011 at 8:46 AM, Richard Mathar
<mathar at strw.leidenuniv.nl> wrote:
>
> In an attempt to extend http://oeis.org/A003241, I took the generating
> function of A003238, shifted by one index to generate the function P(n)
> of Harary-Robinson, then used eq (37) by Harary-Robinson
> R(x) = 2x+P^2(x)/x^2 + (1-x)*P(x)*[P(x^2)/x^2-1]/x
>   -[(P^2(x)-P(x^2))/(2x)+P(x^3)/x^2]
>   -[P^2(x^2)-P(x^4)]/(2x^3)
> and got for R(x) the coefficients
> 1,1,2,2,8,15,24,45,71,106,168,247,345,503,700,934,...
> This is not what Harary and Robinson put into table 1, column R_n,
> which is -- I think-- reproduced in A003241.
> Where is the error? Is this a typo in the paper, my misunderstanding
> of annotations found in eqs (39)-(44), or a typo in my Maple program?
>
> L := BFILETOLIST("b003238.txt") ;
> P := add( op(i,L)*x^(i+1),i=1..120) ;
> Px2 := subs(x=x^2,P) ;
> Px3 := subs(x=x^3,P) ;
> Px4 := subs(x=x^4,P) ;
> R := 2*x+P^2/x^2+(1-x)*P/x*(Px2/x^2-1)-(P^2-Px2)/2/x -Px3/x^2-(Px2^2-Px4)/2/x^3 ;
> taylor(R,x=0,40) ;
>
>
> _______________________________________________
>
> Seqfan Mailing list - http://list.seqfan.eu/
>