[seqfan] Re: Partition numbers

[franktaw at netscape.net]

I'm not sure I entirely follow your program (I don't use Mathematica), 
but take a look at A008284 and A026820, and see whether either of them 
matches what you are doing.

Franklin T. Adams-Watters

-----Original Message-----
From: Mats Granvik <mgranvik at abo.fi>

With the Mathematica program:

Clear[nn, t, n, k, p];
t[n_, 1] = 0;
t[1, k_] = 0;
t[1, 1] = 1;
t[n_, k_] :=
   t[n, k] = If[n >= k, t[n - k + 1, k - 1] + t[n - k + 1, k], 0];
(*nn is size of matrix*)
nn = 351;
(*Partitions numbers*)
p = DeleteCases[Table[t[n, Floor[1/2 + Sqrt[2*nn]]], {n, 1, nn}], 0]

I get the partition numbers:
{1, 1, 2, 3, 5, 7, 11, 15, 22, 30, 42, 56, 77, 101, 135, 176, 231, \
297, 385, 490, 627, 792, 1002, 1255, 1575, 1958}

Mats



_______________________________________________

Seqfan Mailing list - http://list.seqfan.eu/