Pentagonal Number Theorem - Example Program

Example Program

Here is a simple Python program which computes p(n), the number of partitions, using the recurrence resulting from the pentagonal number theorem.

pentagonal = lambda n : n*(3*n-1)/2 def generalised_pentagonal(n): # 0, 1, -1, 2, -2 if n < 0: return 0 if n%2 == 0: return pentagonal(n/2+1) else: return pentagonal(-(n/2+1)) pt = for n in range (1, 1000+1): r = 0 f = -1 i = 0 while 1: k = generalised_pentagonal(i) if k > n: break if i%2==0: f = -f r += f*pt i += 1 pt.append(r) print pt

Read more about this topic:  Pentagonal Number Theorem

Famous quotes containing the word program:

    According to legend, Dr. Sappington purchased his coffin several years before his death and kept it under his bed, with apples and nuts in it for his visiting grandchildren.
    —Administration in the State of Miss, U.S. public relief program (1935-1943)

    Typical of Iowa towns, whether they have 200 or 20,000 inhabitants, is the church supper, often utilized to raise money for paying off church debts. The older and more conservative members argue that the “House of the Lord” should not be made into a restaurant; nevertheless, all members contribute time and effort, and the products of their gardens and larders.
    —For the State of Iowa, U.S. public relief program (1935-1943)