求亲和数

lizuyu

def a(n):
    a1=0
    c=int(n**0.5)
    for i in range(1,c+1):
        if n%i==0 and n/i!=i:
            a1=a1+i+n//i
        elif n%i==0 and n/i==i:
            a1=a1+i
    return a1-n;
b=2
while 1:
    if a(a(b))==b and a(b)<b:   
        print(a(b),b)
        break
    b+=1