472,976 Members | 1,447 Online

# Taxicab Numbers

I'm new to Python and have been putting my mind to learning it over my
holiday break. I've been looking over the functional programming
aspects of Python and I'm stuck trying to come up with some concise
code to find Taxicab numbers (http://mathworld.wolfram.com/
TaxicabNumber.html).

"Taxicab(n), is defined as the smallest number that can be expressed
as a sum of two positive cubes in n distinct ways"

In Haskell something like this could easily be done with:
cube x = x * x * x

taxicab n = [(cube a + cube b, (a, b), (c, d))
| a <- [1..n],
b <- [(a+1)..n],
c <- [(a+1)..n],
d <- [(c+1)..n],
(cube a + cube b) == (cube c + cube d)]

Output::
*Maintaxicab 10
[]
*Maintaxicab 12
[(1729,(1,12),(9,10))]
*Maintaxicab 20
[(1729,(1,12),(9,10)),(4104,(2,16),(9,15))]

I'm looking for a succinct way of doing this in Python. I've been
toying around with filter(),map(), and reduce() but haven't gotten
anything half-way decent yet.

So, how would you implement this taxicab(n) function in Python?
Dec 27 '07 #1
3 3844 On Dec 27, 9:48*pm, rgalgon <rgal...@gmail.comwrote:
I'm new to Python and have been putting my mind to learning it over my
holiday break. I've been looking over the functional programming
aspects of Python and I'm stuck trying to come up with some concise
code to find Taxicab numbers (http://mathworld.wolfram.com/
TaxicabNumber.html).

"Taxicab(n), is defined as the smallest number that can be expressed
as a sum of two positive cubes in n distinct ways"

In Haskell something like this could easily be done with:
cube x = x * x * x

taxicab n = [(cube a + cube b, (a, b), (c, d))
* * * * * * | a <- [1..n],
* * * * * * * b <- [(a+1)..n],
* * * * * * * c <- [(a+1)..n],
* * * * * * * d <- [(c+1)..n],
* * * * * * * (cube a + cube b) == (cube c + cube d)]

Output::
*Maintaxicab 10
[]
*Maintaxicab 12
[(1729,(1,12),(9,10))]
*Maintaxicab 20
[(1729,(1,12),(9,10)),(4104,(2,16),(9,15))]

I'm looking for a succinct way of doing this in Python. I've been
toying around with filter(),map(), and reduce() but haven't gotten
anything half-way decent yet.

So, how would you implement this taxicab(n) function in Python?
Python's list comprehensions are quite similar to haskell's, so this
code can be translated almost word-for-word.

def cube(x):
return x * x * x

def taxicab(n):
return [(cube(a) + cube(b), (a, b), (c, d))
for a in range(1, n + 1)
for b in range(a + 1, n + 1)
for c in range(a + 1, n + 1)
for d in range(c + 1, n + 1)
if cube(a) + cube(b) == cube(c) + cube(d)]

for n in (10, 12, 20):
print list(taxicab(n))

--
Paul Hankin
Dec 28 '07 #2
>>>>"rgalgon" == rgalgon <rg*****@gmail.comwrites:

rgalgonI'm new to Python and have been putting my mind to learning it
rgalgonover my holiday break. I've been looking over the functional
rgalgonprogramming aspects of Python and I'm stuck trying to come up with
rgalgonsome concise code to find Taxicab numbers
rgalgon(http://mathworld.wolfram.com/ TaxicabNumber.html).

rgalgon"Taxicab(n), is defined as the smallest number that can be
rgalgonexpressed as a sum of two positive cubes in n distinct ways"

rgalgonIn Haskell something like this could easily be done with:
rgalgoncube x = x * x * x

rgalgontaxicab n = [(cube a + cube b, (a, b), (c, d))
rgalgon| a <- [1..n],
rgalgonb <- [(a+1)..n],
rgalgonc <- [(a+1)..n],
rgalgond <- [(c+1)..n],
rgalgon(cube a + cube b) == (cube c + cube d)]

rgalgonI'm looking for a succinct way of doing this in Python. I've been
rgalgontoying around with filter(),map(), and reduce() but haven't gotten
rgalgonanything half-way decent yet.

rgalgonSo, how would you implement this taxicab(n) function in Python?

To give a fairly literal translation of your Haskell, you could just use
this:

def cube(n): return n ** 3

def taxicab(n):
n += 1
return [(cube(a) + cube(b), (a, b), (c, d))
for a in xrange(1, n)
for b in xrange(a + 1, n)
for c in xrange(a + 1, n)
for d in xrange(c + 1, n)
if cube(a) + cube(b) == cube(c) + cube(d)]

If you print the return value of this you get the same results as your

if __name__ == '__main__':
import sys
print taxicab(int(sys.argv))

Then I see

\$ time taxicab.py 20
[(1729, (1, 12), (9, 10)), (4104, (2, 16), (9, 15))]

real 0m0.098s
user 0m0.046s
sys 0m0.046s
Terry
Dec 28 '07 #3
I hate having to correct myself...

On Fri, 28 Dec 2007 02:31:50 +0000, Steven D'Aprano wrote:
You seem to be treating n as the number to be split into the sum of
cubes,
Not quite... in your code, n appears to be the largest number to test.
That is, if n is twelve, you test up to 12 cubed.
but that's not what n is supposed to be. n is the number of
distinct sums.
That remains true.

--
Steven.
Dec 28 '07 #4

This thread has been closed and replies have been disabled. Please start a new discussion.

### Similar topics

 4 by: August1 | last post by: A handful of articles have been posted requesting information on how to use these functions in addition to the time() function as the seed to generate unique groups (sets) of numbers - each group... 16 by: StenKoll | last post by: Help needed in order to create a register of stocks in a company. In accordance with local laws I need to give each individual share a number. I have accomplished this by establishing three tables... 12 by: Jim Michaels | last post by: I need to generate 2 random numbers in rapid sequence from either PHP or mysql. I have not been able to do either. I get the same number back several times from PHP's mt_rand() and from mysql's... 7 by: Gerard Flanagan | last post by: All would anyone happen to have code to generate Cutter Numbers: eg. http://www1.kfupm.edu.sa/library/cod-web/Cutter-numbers.htm or is anyone looking for something to do?-) (I'm under... 10 by: LuTHieR | last post by: Hi, I'm reading a string of numbers from a file (using Borland C++ Builder 6), and I'm doing it like this: first I use FileRead to store all the data in the file to a char* variable (appropriately... 7 by: newstips6706 | last post by: 1, 2, 3, 5, 7... PRIME Numbers ________________________________ Definitions What is a PRIME Number ? 17 by: Ron | last post by: I want to write a program that will accept a number in a textbox for example 23578 and then in a label will display the sum of the odd and even number like this... the textbox containsthe number... 24 by: pereges | last post by: I need to generate two uniform random numbers between 0 and 1 in C ? How to do it ? I looked into rand function where you need to #define RAND_MAX as 1 but will this rand function give me ... 0 by: lllomh | last post by: Define the method first this.state = { buttonBackgroundColor: 'green', isBlinking: false, // A new status is added to identify whether the button is blinking or not } autoStart=()=>{ 0 by: Aliciasmith | last post by: In an age dominated by smartphones, having a mobile app for your business is no longer an option; it's a necessity. Whether you're a startup or an established enterprise, finding the right mobile app... 0 by: tracyyun | last post by: Hello everyone, I have a question and would like some advice on network connectivity. I have one computer connected to my router via WiFi, but I have two other computers that I want to be able to... 2 by: giovanniandrean | last post by: The energy model is structured as follows and uses excel sheets to give input data: 1-Utility.py contains all the functions needed to calculate the variables and other minor things (mentions... 3 by: NeoPa | last post by: Introduction For this article I'll be using a very simple database which has Form (clsForm) & Report (clsReport) classes that simply handle making the calling Form invisible until the Form, or all... 1 by: Teri B | last post by: Hi, I have created a sub-form Roles. In my course form the user selects the roles assigned to the course. 0ne-to-many. One course many roles. Then I created a report based on the Course form and... 0 by: isladogs | last post by: The next Access Europe meeting will be on Wednesday 1 Nov 2023 starting at 18:00 UK time (6PM UTC) and finishing at about 19:15 (7.15PM) Please note that the UK and Europe revert to winter time on... 0 by: NeoPa | last post by: Introduction For this article I'll be focusing on the Report (clsReport) class. This simply handles making the calling Form invisible until all of the Reports opened by it have been closed, when it... 4 by: GKJR | last post by: Does anyone have a recommendation to build a standalone application to replace an Access database? I have my bookkeeping software I developed in Access that I would like to make available to other...