%I#14 2023年8月24日16:22:39
%S 1,2,3,4,5,6,7,4,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,12,25,26,
%电话:9,28,29,30,31,16,33,34,35,36,37,38,39,20,41,42,43,44,45,46,47,48,49,
%U 50,51,52,53,18,55,28,57,58,59,60,61,62,63,16,65,66,67,68
%N a(N)是N的最大除数,其素因式分解的指数均为2的幂(A138302)。
%H Amiram Eldar,n的表,a(n)表示n=1..10000</a>
%F与a(p^e)的乘积=p^(2^ floor(log_2(e)))。
%F a(n)=n,当且仅当n在A138302中时。
%F求和{k=1..n}a(k)~c*n^2,其中c=0.4616988732…=(1/2)*Product_{p素数}(1+Sum_{k>=1}(p^F(k)-p^(F(k-1)+1))/p^(2*k)),F(k)=2^floor(log_2(k))和F(0)=0。
%ea(27)=9,因为9=3^2是27的最大除数,在素因式分解中有一个指数,即2的幂。
%t f[p_,e_]:=p^(2^层[Log2[e]]);a[1]=1;a[n_]:=倍@@f@@FactorInteger[n];数组[a,100]
%Y参考A138302,A353898。
%Y类似序列:A000265、A007947、A008834、A055071、A350390。
%K nonn,多个
%O 1,2号机组
%2022年5月10日
|