%I#15 2019年9月3日20:35:19

%S 0，4，9，16，25，13，49，64，81，9121，25169，53，34256289，85361，41，58125，

%电话：529,73625173729,65841,389611024130293,74,971369365178,89，

%电话：1681,621849137106533220926524016292981852809733146113

%N a（N）=A008475（N^2）。

%H Alois P.Heinz，n的表格，n=1..10000的a（n）</a>

%F a（p^e）=p^（2e）的加法。

%e a（1）=0，因为1没有素因子；n=p^2:a（p^2）=p^2；n=6:a（6）=4+9=13；如果gcd（u，w）=1，a（u*w）=a（u）+a（w）；a（21）=a（7）+a（3）=49+9=58；关于酉素因子分解的加法。

%p a:=n->添加（i[1]^i[2]，i=ifactors（n^2）[2]）：

%p序列（a（n），n=1..60）；#_Alois P.Heinz，2019年9月3日

%t ffi[x_]：=展平[FactorInteger[x]]lf[x_]:=长度[FactorInteger[x]]ba[x_】：=表[Part[ffi[x]，2*w-1]，{w，1，lf[x]}]ep[x_'：=表[Part[ffi[x]，2*w]，{v，1，lp[x]{x]}]supo[x__]：=Apply[Plus，ba[x]^ep[x]]表[supo[w]

%o（PARI）f（n）={my（f=因子（n））；向量（向量（#f~，i，f[i，1]^f[i、2]））；}；\\A008475型

%o a（n）=f（n^2）；\\_米歇尔·马库斯，2019年9月3日

%Y参考A008475、A000142、A081403和A081404。

%K nonn公司

%O 1,2号机组

%A _Labos Elemer，2003年3月31日