%I#44 2022年11月19日04:44:17

%S 1,1,1,1,3,1,1,5,1,3,1,7,5,1,1,9,1,5,7,11,1,31,13,1,7,1,15,1,1，

%电话：11,17,7,9,1,19,13,5,12,1,11,5,23,1,3,1,25,17,13,1,27,11,7,19,29,1，

%U 15,1,31,7,1,13,33,1,17,23,35,1,9,1,37,25,19,11,39,1,5,1,41,1,21

%N如果N=p_1^e_1*…*p_k^e_k，p_1<…<p_k素数，则a（n）=n/p_1^e_1，其中a（1）=1。

%C与A067029一起用于定义与a（p^e）=f（e）相乘的序列，作为形式的递归：a（1）=1，对于n>1，a（n）=f_Antti Karttunen_，2017年5月29日

%H Reinhard Zumkeller，n的表，n=1..10000的a（n）</a>

%F a（n）=n/A028233（n）。

%F A0001221（a（n））=A0001221（n）-1；A001222（a（n））=A001222-A067029（n）.-_Reinhard Zumkeller_，2006年5月13日

%F求和{k=1..n}a（k）~c*n^2，其中c=（1/2）*Sum_{k>=0}A005867（k）/（素数（k+1）*（素数）+1）*A002110（k））=0.114813….-_阿米拉姆·埃尔达尔，2022年11月19日

%t a[n_]：=n/Power@@First[FactorInteger[n]]；表[a[n]，{n，1，84}]（*Jean-François Alcover_，2012年6月12日*）

%o（哈斯克尔）

%o a028234 n=n`div`a028233 n---Reinhard Zumkeller_，2013年3月27日

%o（PARI）a（n）={my（f=因子（n））；如果（#f~，f[1,1]=1）；因子回复（f）；}\\米歇尔·马库斯，2016年2月11日

%o（Python）

%o来自sympy进口保理商

%o定义a（n）：

%o f=因子（n）

%o如果n==1则返回1，否则n/（min（f）**f[min（f）]）#_Indranil Ghosh，2017年5月12日

%o（方案）（定义（A028234 n）（/n（A028232 n））；；还需要A020639和A028233的代码_Antti Karttunen_，2017年5月29日

%o（GAP）a:=列表（列表（列表）（列表（[1..10^3]，因子），已收集），i->i[1]），j->j[1]^j[2]）；；

%o A028234:=列表（[1..长度（a）]，i->i/a[i]）；#_Muniru A Asiru_，2018年1月27日

%Y参见A001221、A001222、A002110、A005867、A020639、A028233、A067029。

%K nonn，很好，很容易

%O 1,6型

%电弧勒布朗（_M）_

%E编辑后的名称包含a（1）=1，作者：Franklin T.Adams-Waters_，2018年1月27日