%I#12 2017年3月30日21:58:46
%S 0,1,2,3,5,6,10,10,19,15,19,21,24,28,39,33,53,44,49,53,60,61,69,72,79,
%电话82,92,93117108115140121174146205155233217267192295,
%电话:20922522223824926727029929033631337332841134745138048640053442244741460460511479650455529602553579577626612681632747665796695
%N当从k=N^2开始并重复应用将k替换为k-d(k)的映射时,达到零所需的步长数,其中d(k)是k的除数(A000005)。
%H Antti Karttunen,n表,n=0..5792的a(n)</a>
%F a(n)=A155043(A000290(n))=A165043(n^2)。
%t f[n_]:=长度[NestWhileList[#-除数Sigma[0,#]&,n^2,#!=0&]]-1;f/@范围[0,85](*_Ivan N.Ianakiev,2015年9月25日*)
%o(PARI)
%o分配((2^31)+(2^30));
%o上传=2^25;
%o v155043=矢量(uplim);
%o v155043[1]=1;v155043[2]=1;
%o表示(i=3,uplim,v155043[i]=1+v155043[i-numdiv(i)];如果(!(i%65536),打印1(i,“,”););
%o A155043=n->如果(!n,n,v155043[n]);
%o A261088=n->A155043(n^2);
%o表示(n=0,5792,写入(“b261088.txt”,n,“”,A261088(n));
%o(方案)(定义(A261088 n)(A155043(A000290 n))
%Y参见A000005、A000290、A155043、A049820、A259934、A261085。
%Y另请参阅A260732、A261222。
%K nonn公司
%0、3
%2015年9月23日,安蒂·卡图宁
|