primes = 2 : [n | n <- [3..], isPrime n]
    where isPrime n = foldr (\p r-> p*p>n || (rem n p /= 0 && r)) True primes
 
sumDiv n = (divs n primes) - n where
    divs 1 _ = 1
    divs n (p:ps) =
        if p*p > n
            then n+1
            else let (c, r) = cnt n p 0 in
                if c == 0
                    then divs n ps
                    else (sum $ take (c+1) $ iterate (*p) 1) * (divs r ps)
    cnt n p c | mod n p == 0 = cnt (div n p) p (c+1)
              | otherwise = (c, n)
 
friends :: [(Integer, Integer)]
friends = filter (\(x, s) -> x < s && sumDiv s == x) $ map (\x -> (x, sumDiv x)) [1..]
 
main = mapM print $ take 37 friends

cHJpbWVzID0gMiA6IFtuIHwgbiA8LSBbMy4uXSwgaXNQcmltZSBuXQogICAgd2hlcmUgaXNQcmltZSBuID0gZm9sZHIgKFxwIHItPiBwKnA+biB8fCAocmVtIG4gcCAvPSAwICYmIHIpKSBUcnVlIHByaW1lcwoKc3VtRGl2IG4gPSAoZGl2cyBuIHByaW1lcykgLSBuIHdoZXJlCiAgICBkaXZzIDEgXyA9IDEKICAgIGRpdnMgbiAocDpwcykgPQogICAgICAgIGlmIHAqcCA+IG4KICAgICAgICAgICAgdGhlbiBuKzEKICAgICAgICAgICAgZWxzZSBsZXQgKGMsIHIpID0gY250IG4gcCAwIGluCiAgICAgICAgICAgICAgICBpZiBjID09IDAKICAgICAgICAgICAgICAgICAgICB0aGVuIGRpdnMgbiBwcwogICAgICAgICAgICAgICAgICAgIGVsc2UgKHN1bSAkIHRha2UgKGMrMSkgJCBpdGVyYXRlICgqcCkgMSkgKiAoZGl2cyByIHBzKQogICAgY250IG4gcCBjIHwgbW9kIG4gcCA9PSAwID0gY250IChkaXYgbiBwKSBwIChjKzEpCiAgICAgICAgICAgICAgfCBvdGhlcndpc2UgPSAoYywgbikKCmZyaWVuZHMgOjogWyhJbnRlZ2VyLCBJbnRlZ2VyKV0KZnJpZW5kcyA9IGZpbHRlciAoXCh4LCBzKSAtPiB4IDwgcyAmJiBzdW1EaXYgcyA9PSB4KSAkIG1hcCAoXHggLT4gKHgsIHN1bURpdiB4KSkgWzEuLl0KCm1haW4gPSBtYXBNIHByaW50ICQgdGFrZSAzNyBmcmllbmRz