%I#23 2022年11月12日18:22:29
%S 1,1,2,2,1,1,4,4,1,1,4,8,3,1,8,5,1,4,1,8,1,1,1,8,16,1,9,12,8,1,16,16,
%T 7,1,12,16,1,3,6,2,1,4,16,1,1,16,32,27,2,18,1,5,12,4,11,12,24,5,
%U 1,32,32,1,1,4,24,3,1,24,32,15,1,4,16,3,27,2,2,1,1,1,12,32,9型
%N a(1)=1;对于n>1,a(n)=a(n-1)+n除以其最大素因子。
%H Rémy Sigrist,n的表,n=1..10000的a(n)</a>
%e对于n=4,a(4)=2+4除以其最大素因子=6/3=2。
%t f[n_]:=n/因子整数[n][[-1,1]];a[1]=1;a[n]:=a[n]=f[a[n-1]+n];阵列[a,100](*_Amiram Eldar_,2020年1月16日*)
%tnxt[{n_,a_}]:={n+1,(a+n+1)/因子整数[a+n+1][[-1,1]]};NestList[nxt,{1,1},90][[All,2]](*哈维·P·戴尔,2022年11月12日*)
%o(Magma)[n eq 1 select 1 else(Self(n-1)+n)div Max(PrimeDivisors(Self-1)+n)):n in[1..85]];//_Marius A.Burtea,2020年2月17日
%o(PARI)a(n)=如果(n==1,1,my(x=a(n-1)+n);x/vecmax(系数(x)[,1]);\\_米歇尔·马库斯,2020年2月20日
%Y参考A006530(最大素因数),A208884。
%K nonn公司
%氧1,3
%A _Ali Sada,2020年1月16日
|