%I#19 2017年5月28日09:14:19

%S 1,1,1,4,1,1,1,27,4,1,1,4,1,1,1,1256,1,4,4,11,1,27,4,4],1,1,1,1，

%电话3125,1,1,1,16,1,1,27,1,1,4,4,1,1256,4,1,4,14,1,271,1,1,4,1，

%U 1,446656,1,1,1,4,1,1108,1,1,4],1,1256256,1,1,1,4,1,1,1,12,27,1,4,4,1,1,13125,1,44,16,1,1,11,27,1

%如果p^b（p，N）是素数p除以N的最大幂，那么a（N）=product_{p|N}b（p、N）^b（p，N）。

%C a（1）=1（空产品）。

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

%e 45=3^2*5^1。所以a（45）=2^2*1^1=4。

%e 72=2^3*3^2。所以a（72）=3^3*2^2=108。

%t f[n_]：=块[{fi=最后@换位@因子整数@n}，Times@@（fi^fi）]；休息@阵列[f，93]（*_罗伯特·G·威尔逊v_*）

%o（PARI）a（n）=局部（v，r，i）；v=因子（n）；r=1；对于（i=1，matsize（v）[1]，r*=v[i，2]^v[i、2]）；r（Herrgesell）

%Y参考A112621、A112624。

%K nonn，多个

%O 1,4型

%A _罗伊查询，2005年12月25日

%E更多条款来自_Robert G.Wilson v _和Lambert Herrgesell（zero815（AT）googlemail.com），2005年12月27日

%E修正了起始偏移量，数据部分扩展到105个术语-Antti Karttunen，2017年5月28日