Задача №635. Дружественные числа

def sd(k):
    s = 0
    m = int(k**0.5)
    for i in range(1, m+1):
        if k%i == 0:
            if i**2 == k:
                s += i    
            else:
                s+=i
                s+=(k//i)
    s -= k
    return s
 
a = input().split()
n,m = int(a[0]), int(a[1])
f = 0
for i in range(n, m):
    for j in range(i+1, m+1):
        if sd(i) == j and sd(j) == i:
            print(i, j)
            f = 1
if f == 0:
    print('Absend')

def fun(num):
    res = 1
    d = 2
    while d*d <= num:
        if num%d == 0:
            res += d + num//d*(d*d != num)
        d += 1
    return res 
 
m, n = map(int,input().split())
lst = set()
for i in range(m,n + 1):
    k = fun(i)
    if i < k <= n and fun(k) == i:
        lst.add((i,k))
if len(lst):
    for i in sorted(lst):
        print(*i)
else:
    print('Absent')

def primes(n):
    sieve = [True] * n 
    for i in range(3, int(n ** 0.5) + 1, 2):
        if sieve[i]:
            sieve[i * i::2 * i] = [False] * ((n - i * i - 1) // (2 * i) + 1)
    return [2] + [i for i in range(3, n, 2) if sieve[i]] 
 
 
def fact(n):
    res = []
    for d in pr:
        if d*d > n :
            break
        cnt = 0
        while not n % d:
            n //= d
            cnt += 1
        if cnt > 0:
            res.append((d, cnt))
    if n > 1:
        res.append((n, 1))
 
    return res
 
 
def divs(n):
 
    res = [1]
    for v, i in fact(n):
        res += [x*v**k for k in range(1,i+1) for x in res]
    return sum(res) - n
 
 
 
m, n = map(int,input().split())
pr = primes(n)
lst = set()
 
for i in range(m,n + 1):
    k = divs(i)
    if i < k <= n and divs(k) == i:
        lst.add((i,k))
if len(lst):
    for i in sorted(lst):
        print(*i)
else:
    print('Absent')