s[n_]:=除数Sigma[1,n]-n;BetrothedNumberQ[n_]:=如果[s[s[n]-1]==n+1&&n>1,对,错];BetrothedPairList[k_]:=(anlist=Select[Range[k],Betrothed NumberQ[#]&];prlist=Table[Sort[{anlist[[n]],s[anlist[n]]-1}],{n,1,Length[anlist]}];Union[prlist,prlist]);data=BetrothedPairList[10^6];表[Length[Select[data,First[#]<10^k&]],{k,1,6}]