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:

    “Like other cities created overnight in the Outlet, Woodward acquired between noon and sunset of September 16, 1893, a population of five thousand; and that night a voluntary committee on law and order sent around the warning, “if you must shoot, shoot straight up!””
    —State of Oklahoma, U.S. public relief program (1935-1943)

    “In 1862 the congregation of the church forwarded the church bell to General Beauregard to be melted into cannon, “hoping that its gentle tones, that have so often called us to the House of God, may be transmuted into war’s resounding rhyme to repel the ruthless invader from the beautiful land God, in his goodness, has given us.””
    —Federal Writers’ Project Of The Wor, U.S. public relief program (1935-1943)