python亲密数设计

'''
亲密数
　　　　（如果a的所有正因子和等于b，b的所有正因子和等于a，因子包括1但不包括本身，且a不等于b，则称a，b为亲密数对。一般通过叠代编程求出相应的亲密数对）
'''
n = 3000
def fun01(n):
    s = 0
    for i in range(1,n):
        if n % i == 0:
            s += i
    return s
for i in range(1,n):
    a = fun01(i)
    if i == fun01(a) and i < a:
        print(i,a)