亲密数就是在完数的基础上更进一步

for i in range(2,3001):

    s1=s2=1

    for j in range(2,i//2+1):

        if i%j==0:

            s1+=j

    for k in range(2,s1//2+1):

        if s1%k==0:

            s2+=k

    if i==s2 and i!=s1 and i<s1:#小于就会不等于我真是多此一举

        print('(%d,%d)'%(i,s1),end='')