Revision history [back]
 | 1 | initial version |
Here is a possible approach:
def all_two_squares(n):
return [(abs(d[0]),abs(d[1])) for d in divisors(GaussianIntegers()(n)) if norm(d)==n]
| | 1 | initial version |
Here is a possible approach:
def all_two_squares(n):
return [(abs(d[0]),abs(d[1])) for d in divisors(GaussianIntegers()(n)) if norm(d)==n]
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.