%I#17 2023年10月22日16:35:57
%S 1,3,4,7,13,25,33,48,61,79,91119133157181212230269289331363,
%电话:39942348351455659665268275478684979591198910901128,
%电话:1188124413341376147215161600167817501798192219792072207221442242
%如果σ(N)>=a(N-1),则N a(1)=1,a(N)=σ。
%C a(n)=σ(n)仅适用于n=1,2,3,4吗?
%H Harvey P.Dale,n的表格,n=1..1000的a(n)</a>
%F a(n)=A024916(n)-8和A160664(n)-9,对于n>3_格雷格·德累斯顿,2020年2月23日
%tnxt[{n_,a_}]:=模[{s=DivisorSigma[1,n+1]},{n+1,If[s>=a,s,a+s]}];NestList[nxt,{1,1},60][[All,2]](*哈维·P·戴尔,2022年8月4日*)
%o(PARI){v=向量(60);v[1]=σ(1);对于(i=2,60,如果(σ(i)<v[i-1],v[i]=v[i-2]+σ
%o(Python)
%o从数学导入isqrt
%o定义A092406(n):返回(-(s:=isqrt(n))**2*(s+1)+总和((q:=n//k)*((k<<1)+q+1),对于范围(1,s+1)>>1中的k)-8,如果n>3其他(1,如果n<2其他n+1)#_Chai Wah Wu_,2023年10月22日
%Y参考A0000203、A024916、A160664。
%K nonn公司
%O 1,2号机组
%A _Jon Perry_,2004年3月22日