/*
Euler #21 in Picat.
Problem 21
"""
Let d(n) be defined as the sum of proper divisors of n (numbers less
than n which divide evenly into n).
If d(a) = b and d(b) = a, where a /= b, then a and b are an amicable
pair and each of a and b are called amicable numbers.
For example, the proper divisors of 220 are
1, 2, 4, 5, 10, 11, 20, 22, 44, 55 and 110; therefore d(220) = 284.
The proper divisors of 284 are 1, 2, 4, 71 and 142; so d(284) = 220.
Evaluate the sum of all the amicable numbers under 10000.
"""
This Picat model was created by Hakan Kjellerstrand, hakank@gmail.com
See also my Picat page: http://www.hakank.org/picat/
*/
main => go.
go => time(euler21).
euler21 =>
S = new_map(),
foreach(A in 1..10000-1)
B = sum_divisors2(A),
C = sum_divisors2(B),
if A != B, A == C then
S.put(A, 1),
S.put(B, 1)
end
end,
% writeln(s=S),
println(sum(S.keys())),
nl.
table
sum_divisors2(N) = Sum =>
D = floor(sqrt(N)),
Sum1 = 1,
foreach(I in 2..D, N mod I == 0)
Sum1 := Sum1+I,
if I != N div I then
Sum1 := Sum1 + N div I
end
end,
Sum = Sum1.