%I#11 2023年11月2日10:27:43
%S 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,8,17,18,19,20,21,22,23,24,25,26,
%电话:27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,24,49,
%U 50,51,52,53,54,55,56,57,58,59,60,61,62,63,32,65,66,67,68
%N非A322448项的N的最大除数。
%C First与A365683不同,n=64。
%n的最大因子,其素因式分解的指数都是1或素数。
%C这些除数的个数是A366991(n),它们的和是A3669962(n)。
%H Amiram Eldar,n的表,a(n)表示n=1..10000</a>
%F与a(p)=p和a(p^e)=p^A007917(e)相乘,e>=2。
%F a(n)<=n,当且仅当n不在A322448中时相等。
%F Sum_{k=1..n}a(k)~c*n^2,其中c=(1/2)*Product_{p prime}F(1/p)=0.8535795387619596052…,其中F(x)=(1-x)*(1+Sum_{k>=1}x^(2*k-s(k))),当k>=2时,s(k)=A0071917(k),s(1)=1。
%t f[p_,e_]:=p^如果[e==1,1,NextPrime[e+1,-1]];a[1]=1;a[n_]:=倍@@f@@FactorInteger[n];数组[a,100]
%o(PARI)a(n)={my(f=因子(n)
%Y参见A007917、A322448、A365683、A366989、A366991和A366992。
%K nonn,简单,多
%O 1,2号机组
%A阿米拉姆·埃尔达尔,2023年10月31日
|