[seqfan] Re: checking a sequence before submission
David Newman
davidsnewman at gmail.com
Thu Dec 23 16:11:15 CET 2010
My thanks to those who checked this sequence. I'll write it up and
submit it now.
On Wed, Dec 22, 2010 at 1:26 PM, Douglas McNeil <mcneil at hku.hk> wrote:
> In very unoptimized Sage:
>
> ## testing DN's seq
> ##
>
> N = 38
> s = 1
> for n in range(1, N):
> parts = list(Partitions(n, max_slope=-1))
> max_length = max(len(part) for part in parts)
> longest = list(part for part in parts if len(part) == max_length)
> s += sum( (x**n / prod(1-x**i for i in part)) for part in longest)
>
> t = taylor(s, x, 0, N-1)
> c = [t.coeff(x,i) for i in range(N)]
>
> sage: c
> [1, 1, 2, 2, 4, 6, 7, 10, 14, 20, 24, 32, 40, 54, 69, 86, 106, 135,
> 165, 206, 256, 311, 378, 460, 555, 670, 808, 970, 1156, 1380, 1638,
> 1938, 2296, 2706, 3188, 3752, 4390, 5136]
>
>
> Doug
>
> --
> Department of Earth Sciences
> University of Hong Kong
>
> _______________________________________________
>
> Seqfan Mailing list - http://list.seqfan.eu/
>
More information about the SeqFan
mailing list