Skip to main content

Heap's Algorithm implementation with python

Heap's Algorithm is used to generate all the permutations of a given object. It was first proposed by B.R. Heap in the year 1963.

In this post, we will write a program using Heap's Algorithm which will generate permutations for a given object. Remember that the object should be iterable. Otherwise, you will get an error.

If you don't know what permutations is all about then you can visit some of the following websites:

Permutations - Wikipedia
Permutation - Wolfram math
Permutations and combinations - Math is fun  and
Permutations Khan Academy

Here goes the program:
This function will take an iterable object as input and will generate all the permutations possible with the items in the object.
For example, if the input is:heap(["A", "B", "C"]) then the output will be a generator object. You can use a for loop to iterate over this object to see the output.
But remember that if you have any duplicates in the input(Ex: heap(["A", "B", "B"])) then you will also get duplicate permutations. (Ex output: [["A", "B", "B"], ["B", "A", "B"], ["B",  "A", "B"], ["A", "B", "B"], ["B", "B", "A"], ["B", "B", "A"]])

You can eliminate duplicates by using python sets.

One of the best ways to understand this algorithm is by looking at the pseudo-code on Wikipedia Heaps Algorithm page.

We have used the nonrecursive method.

If you are not sure of what yield does, then you can have a look at this answer on Stack overflow by e-satis for the question What does the "yield" keyword do

If you use the above code and if you give a string as an input still you will get a list as output. To overcome this problem I have written a program, a slightly modified version of the above program so that we will get the output as a string.
The program is as follows:   
The above program will take a string as input and will give a generator object as output, which on further iteration will give you a list of strings.

Line 14 is just swapping of the characters a[0] , a[i] in the string.

Line 16 is again swapping of the characters a[i], a[c[i]] in the string. 

As we are doing string concatenation with more than two strings, the second program takes more time for execution than the first one.

If you have any doubt or didn't understand anything, then you can comment in the comment box below and I will be glad to help you.

If you have any suggestions or if you have found any typo then please comment in the comment box. I will be very happy to see your comment.

You can also contact me.

Thank you. Have a nice day. 😃

Popular posts from this blog

Project Euler Problem 62 solution with python

Cubic permutations ¶ The cube, 41063625 (3453), can be permuted to produce two other cubes: 56623104 (3843) and 66430125 (4053). In fact, 41063625 is the smallest cube which has exactly three permutations of its digits which are also cube. Find the smallest cube for which exactly five permutations of its digits are cube.

Project Euler Problem 67 Solution with Python

Maximum path sum II By starting at the top of the triangle below and moving to adjacent numbers on the row below, the maximum total from top to bottom is 23. 3 7 4 2 4 6 8 5 9 3 That is, 3 + 7 + 4 + 9 = 23. Find the maximum total from top to bottom in triangle.txt (right click and 'Save Link/Target As...'), a 15K text file containing a triangle with one-hundred rows.

Problem 43 Project Euler Solution with python

Sub-string divisibility The number, 1406357289, is a 0 to 9 pandigital number because it is made up of each of the digits 0 to 9 in some order, but it also has a rather interesting sub-string divisibility property. Let d 1 be the 1 st digit, d 2 be the 2 nd digit, and so on. In this way, we note the following: d 2 d 3 d 4 =406 is divisible by 2 d 3 d 4 d 5 =063 is divisible by 3 d 4 d 5 d 6 =635 is divisible by 5 d 5 d 6 d 7 =357 is divisible by 7 d 6 d 7 d 8 =572 is divisible by 11 d 7 d 8 d 9 =728 is divisible by 13 d 8 d 9 d 10 =289 is divisible by 17 Find the sum of all 0 to 9 pandigital numbers with this property. One might write a simple solution using direct if else statements for this problem. Using if else statements, the execution time may be a few seconds. But this is not a very good approach. I too had written a program with if else statement which will take each and every permutation of the 0-9 Pandigital and check for the conditions given in the qu...