Amicable pairs: Difference between revisions
Content deleted Content added
imported>KayproKid m →{{header|S-BASIC}}: fixed missing tag |
Add Sage implementaion |
||
Line 6,025: | Line 6,025: | ||
[1, 2, 4, 8, 16, 23, 46, 47, 92, 94, 184, 188, 368, 376, 752, 1081, 2162, 4324, 8648] |
[1, 2, 4, 8, 16, 23, 46, 47, 92, 94, 184, 188, 368, 376, 752, 1081, 2162, 4324, 8648] |
||
[1, 2, 4, 8, 16, 1151, 2302, 4604, 9208] |
[1, 2, 4, 8, 16, 1151, 2302, 4604, 9208] |
||
</pre> |
|||
=={{header|Sage}}== |
|||
<syntaxhighlight lang="Sage"> |
|||
# Define the sum of proper divisors function |
|||
def sum_of_proper_divisors(n): |
|||
return sum(divisors(n)) - n |
|||
# Iterate over the desired range |
|||
for x in range(1, 20001): |
|||
y = sum_of_proper_divisors(x) |
|||
if y > x: |
|||
if x == sum_of_proper_divisors(y): |
|||
print(f"{x} {y}") |
|||
</syntaxhighlight> |
|||
{{out}} |
|||
<pre> |
|||
Share |
|||
220 284 |
|||
1184 1210 |
|||
2620 2924 |
|||
5020 5564 |
|||
6232 6368 |
|||
10744 10856 |
|||
12285 14595 |
|||
17296 18416 |
|||
</pre> |
</pre> |
||