[seqfan] Partition numbers

Mats Granvik mgranvik at abo.fi
Sat Aug 20 12:19:19 CEST 2011


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





More information about the SeqFan mailing list