登录
e-酉Euler函数:a(1)=1,a(n)=n=乘积p(i)^e(i)的乘积uphi(e(i(A047994号).

%我#24 2022年11月29日02:41:30

%S 1,1,1,1,1,1,1,2,1,1,1,1,1,1,1,3,1,1,11,1,1,2,1,1,1,1、1,1,1,2、1,1,1和1,4,1,1,

%在1,1,1,1,1,2,1,1,1,1,1,1,1,1,3,1,1,1,1,1,2,1,1,1,1,1,1,1,1,1,1,1,1,1,2,1,1,1,1,2,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,

%U 1,1,1,2,1,1,1,1,1,1,1,1,1,3,3,1,1,11,1,1,1

%N e-酉Euler函数:a(1)=1,a(N)=乘积uphi(e(i)),N=乘积p(i)^e(i。

%C A072911的统一版本。

%C对于n=乘积p(i)^e(i)>1,a(n)是n的除数d,使得d和n是指数幺正互质,即d=乘积p(i)*f(i),其中1<=f(i。

%H Amiram Eldar,n的表,a(n)表示n=1..10000</a>

%H Nicusor Minculete和LászlóTóth,<a href=“http://ac.inf.elte.hu/Vol_035_2011/elte_annales_35_jav_vagott.pdf#page=205“>指数幺正除数,《布达佩斯科学年鉴》第35卷(2011年)第205-216页。

%F Sum_{k=1..n}a(k)~c1*n+c2*n^(1/3)+O(n ^(1/4+eps)),其中c1=A358658,c2是一个常数(参见Minculete和Tóth,2011)_阿米拉姆·埃尔达尔,2022年11月29日

%tf[p_,e_]:=p^e-1;uphi[1]=1;uphi[n_]:=倍@@f@@FactorInteger[n];fe[p,e]:=uphi[e];eufi[n_]:=倍@@fe@@@FactorInteger[n];数组[eufi,100]

%o(PARI)uphi(n)={my(f=因子(n));prod(i=1,#f~,f[i,1]^f[i、2]-1)};

%o a(n)={my(f=因子(n));prod(i=1,f~,uphi(f[i,2]))};\\_阿米拉姆·埃尔达尔,2022年11月29日

%Y参见A047994、A072911、A358658。

%K非n,简单,多

%O 1,8型

%A _Amiram Eldar,2019年1月10日