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:

    “D--n me, stranger, ef you can’t stay as long as you please, and I’ll give you plenty to eat and drink. Play away, stranger, you kin sleep on the dry spot tonight!”
    —Administration in the State of Arka, U.S. public relief program (1935-1943)

    “O Jesse had a wife, a mourner all her life
    And the children they were brave,
    But the dirty little coward that shot Mr. Howard
    He laid Jesse James in his grave.”
    —Administration in the State of Miss, U.S. public relief program (1935-1943)