%I#15 2019年11月18日14:59:48
%序号0,1,2,8,9,17,46,7918040488219464122
%N具有m<=10^k(且k=1,2,3,…)的订婚对(m,N)的数目,其中订婚对满足西格玛(m)=西格玛(N)=m+N+1且m<N。
%C也称为拟亲和对或约化亲和对。
%H J.O.M.Pedersen,<a href=“http://amicable.homepage.dk/tables.htm“>等分循环表
%H J.O.M.Pedersen,<a href=“http://web.archive.org/web/20140502102524/http://amicable.homepage.dk/tables.htm“>等分周期表
%H J.O.M.Pedersen,等分周期表
%H P.Pollack,<a href=“https://cs.uwaterloo.ca/journals/JIS/VOL14/Pollack/pollack3.html“>拟可数罕见,J.Int.Seq.14(2011)#11.5.2。
%e a(7)=46,因为有46对订婚对(m,n),m<=10^7
%t s[n_]:=除数Sigma[1,n]-n;BetrothedNumberQ[n]:=如果[s[s[n]-1]==n+1&&n>1,True,False];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];表[长度[Select[data,First[#]<10^k&]],{k,1,6}]
%Y参考A003502、A003503、A005276。
%K nonn,难,更多
%O 1,3
%安特金,2006年12月19日
%E a(13)摘自_Giovanni Resta_,2019年7月24日
|