%I#13 2021年11月5日21:12:25

%S 1,3,4,1,6,12,8,9,2,18,12,4,14,24,24,18,6,20,6,32,36,24,8,8，

%电话30,72,32,3,48,54,48,2,38,60,56,54,42,96,44,12,72,48,4,2,6,72,14，

%U 54,12,72,72,80,90,60,24,62,96,16,1,84144,68,18,96144,72,18,74114,8,20,96168,80,6,2126,84,32

%N a（N）=gcd（A003959（N），A034448（N）。

%C这不是乘法。第一点，其中a（m*n）=a（m）*a（n）不适用于互素m，n为1444=2^2*19^2，其中a（1444）=10！=1*2=a（4）*a（361）。有关这些位置的列表，请参见A348740。

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

%F a（n）=gcd（A003959（n），A034448（n））。

%F a（n）=gcd（A003959（n），A348732（n））=gcd（A034448（n），A348732（n））。

%F a（n）=A003959（n）/A348734（n）=A034448（n）/A348735（n）。

%t f1[p_，e_]：=（p+1）^e；f2[p_，e_]：=p^e+1；a[1]=1；a[n_]：=GCD[Times@@f1@@@（f=FactorInteger[n]），Times@@f2@@@f]；阵列[a，100]（*_Amiram Eldar_，2021年11月5日*）

%o（PARI）

%o A003959（n）=｛my（f=因子（n））；对于（i=1，#f~，f[i，1]++）；因子返回（f）；｝；

%o A034448（n）={my（f=因子（n））；prod（k=1，#f~，1+（f[k，1]^f[k，2]））；}；

%o A348733（n）=gcd（A003959（n），A034448（n））；

%Y参考A003959、A034448、A348732、A34873、A34875、A348740。

%Y也可参见A344695、A348047、A348503、A348946中类似的几乎是乘法序列。

%K nonn公司

%O 1,2号机组

%2021年11月5日，安提·卡图内