[seqfan] Re: Help with A002793 needed/Karol A. Penson/2
Karol
penson at lptl.jussieu.fr
Tue Jul 6 00:01:43 CEST 2010
Amazing achievement ! , Thanks , Karol
Richard Mathar a écrit :
> Lisons http://list.seqfan.eu/pipermail/seqfan/2010-July/005230.html :
>
> kp> a(2)=4 =3+1
> kp> a(3)=20 =11+8+1
> kp> a(4)=124=50+58+15+1
> kp> a(5)=920=274+444+177+24+1,
>
> One can get an interpretation if one follows the advice
> given in A002793 to look at this as the bisection of A056952.
>
> The recurrence of A056952(n) can be looked at as a path
> with jumps/skips of length 1 or 2 totaling n, where jumps of length 1
> get weight 1, which is the factor 1 in front of a(n-1),
> and jumps of length 2 get weight floor(n/2), which is the factor in front
> of a(n-2).
>
> Obtain A056952(n) as follows: write down all compositions
> of n such that all addends are either 1 or 2, and such that
> the first addend is 1 (this last restriction comes form A056952(0)=0).
> Then associate with each composition a product of weights
> as computed above, so the n in the weight is the product of
> all floor(s/2) where s is the "partial sum" along the composition at
> the point one sees a 2. Then the terms in the triangle are the
> sums of the weights for all compositions of fixed common length.
>
> Example of n=8:
> [1, 2, 2, 2, 1], 6 = floor(3/2)*floor(5/2)*floor(7/2)
> [1, 2, 2, 1, 2], 8 = floor(3/2)*floor(5/2)*floor(8/2)
> [1, 2, 1, 2, 2], 12 = floor(3/2)*floor(6/2)*floor(8/2)
> [1, 1, 2, 2, 2], 24
> 50 = 6+8+12+24
> [1, 2, 2, 1, 1, 1], 2
> [1, 2, 1, 2, 1, 1], 3
> [1, 1, 2, 2, 1, 1], 6
> [1, 2, 1, 1, 2, 1], 3
> [1, 1, 2, 1, 2, 1], 6
> [1, 1, 1, 2, 2, 1], 6
> [1, 2, 1, 1, 1, 2], 4
> [1, 1, 2, 1, 1, 2], 8
> [1, 1, 1, 2, 1, 2], 8
> [1, 1, 1, 1, 2, 2], 12
> 58 = 2+3+6+3+6+6+4+8+8+12
> [1, 2, 1, 1, 1, 1, 1], 1
> [1, 1, 2, 1, 1, 1, 1], 2
> [1, 1, 1, 2, 1, 1, 1], 2
> [1, 1, 1, 1, 2, 1, 1], 3
> [1, 1, 1, 1, 1, 2, 1], 3
> [1, 1, 1, 1, 1, 1, 2], 4
> 15 = 1+2+2+3+3+4
> [1, 1, 1, 1, 1, 1, 1, 1], 1
> 1 =1
> 124
>
> Example of n=10:
>
> [1, 2, 2, 2, 2, 1], 24
> [1, 2, 2, 2, 1, 2], 30
> [1, 2, 2, 1, 2, 2], 40
> [1, 2, 1, 2, 2, 2], 60
> [1, 1, 2, 2, 2, 2], 120
> 274 = 24+30+40+60+120
> [1, 2, 2, 2, 1, 1, 1], 6
> [1, 2, 2, 1, 2, 1, 1], 8
> [1, 2, 1, 2, 2, 1, 1], 12
> [1, 1, 2, 2, 2, 1, 1], 24
> [1, 2, 2, 1, 1, 2, 1], 8
> [1, 2, 1, 2, 1, 2, 1], 12
> [1, 1, 2, 2, 1, 2, 1], 24
> [1, 2, 1, 1, 2, 2, 1], 12
> [1, 1, 2, 1, 2, 2, 1], 24
> [1, 1, 1, 2, 2, 2, 1], 24
> [1, 2, 2, 1, 1, 1, 2], 10
> [1, 2, 1, 2, 1, 1, 2], 15
> [1, 1, 2, 2, 1, 1, 2], 30
> [1, 2, 1, 1, 2, 1, 2], 15
> [1, 1, 2, 1, 2, 1, 2], 30
> [1, 1, 1, 2, 2, 1, 2], 30
> [1, 2, 1, 1, 1, 2, 2], 20
> [1, 1, 2, 1, 1, 2, 2], 40
> [1, 1, 1, 2, 1, 2, 2], 40
> [1, 1, 1, 1, 2, 2, 2], 60
> 444 = 6+8+12+24+8+...+60
> [1, 2, 2, 1, 1, 1, 1, 1], 2
> [1, 2, 1, 2, 1, 1, 1, 1], 3
> [1, 1, 2, 2, 1, 1, 1, 1], 6
> [1, 2, 1, 1, 2, 1, 1, 1], 3
> [1, 1, 2, 1, 2, 1, 1, 1], 6
> [1, 1, 1, 2, 2, 1, 1, 1], 6
> [1, 2, 1, 1, 1, 2, 1, 1], 4
> [1, 1, 2, 1, 1, 2, 1, 1], 8
> [1, 1, 1, 2, 1, 2, 1, 1], 8
> [1, 1, 1, 1, 2, 2, 1, 1], 12
> [1, 2, 1, 1, 1, 1, 2, 1], 4
> [1, 1, 2, 1, 1, 1, 2, 1], 8
> [1, 1, 1, 2, 1, 1, 2, 1], 8
> [1, 1, 1, 1, 2, 1, 2, 1], 12
> [1, 1, 1, 1, 1, 2, 2, 1], 12
> [1, 2, 1, 1, 1, 1, 1, 2], 5
> [1, 1, 2, 1, 1, 1, 1, 2], 10
> [1, 1, 1, 2, 1, 1, 1, 2], 10
> [1, 1, 1, 1, 2, 1, 1, 2], 15
> [1, 1, 1, 1, 1, 2, 1, 2], 15
> [1, 1, 1, 1, 1, 1, 2, 2], 20
> 177
> [1, 2, 1, 1, 1, 1, 1, 1, 1], 1
> [1, 1, 2, 1, 1, 1, 1, 1, 1], 2
> [1, 1, 1, 2, 1, 1, 1, 1, 1], 2
> [1, 1, 1, 1, 2, 1, 1, 1, 1], 3
> [1, 1, 1, 1, 1, 2, 1, 1, 1], 3
> [1, 1, 1, 1, 1, 1, 2, 1, 1], 4
> [1, 1, 1, 1, 1, 1, 1, 2, 1], 4
> [1, 1, 1, 1, 1, 1, 1, 1, 2], 5
> 24
> [1, 1, 1, 1, 1, 1, 1, 1, 1, 1], 1
> 1
> 920
> 274+444+177+24+1=920
>
> Example of n=12:
> 1764+3708+2016+416+35+1 = 7940
> Example of n=14:
> 13068+33984+23533+6560+835+48+1 = 78040
>
> Apparently the first but last term in each row is A005563.
> Apparently the first column in the triangle is A000254
> Apparently the second column in the triangle is A002538
> The last term in each row is 1.
>
> In Maple:
> n := 10; #or whatever the even index in A056952 may be
> a := 0 ;
> for c from 1 to n do
> fpar := 0 ;
> com := combinat[composition](n,c) ;
> for co in com do
> if max(op(co))<3 and op(1,co) = 1 then
> psu := 1 ;
> f := 1 ;
> for pos from 2 to nops(co) do
> psu := psu+op(pos,co) ;
> if op(pos,co) = 2 then
> f := f*floor(psu/2) ;
> end if;
> end do:
> a := a+f ;
> fpar := fpar+f ;
> print(co,f) ;
> end if;
> end do:
> print(fpar) ;
> end do:
> print(a) ;
>
>
> _______________________________________________
>
> Seqfan Mailing list - http://list.seqfan.eu/
>
More information about the SeqFan
mailing list