%I#14 2022年10月25日02:38:01

%S 1,7,3,19,3,21,3,37,6,21,3,57,3,21,9,61,3,42,3,57,9,21,3111,6,21,10，

%电话57,3,63,3,91,9,21,9114,3,21,911，3,633,57,18,21,3183,6,42,9,57,3，

%U 70,9111,9,21,3171,3,21,18127,9,63,3,57,9,66,3222,3,21，18,57,9，63,3183单位

%N Dirichlet g.f.：（1+4/2^s+1/4^s）*zeta（s）^3。

%H Antti Karttunen，<a href=“/A145399/b145399.txt”>n的表，a（n）表示n=1..65537</a>

%H J.S.卢瑟福，<a href=“https://dx.doi.org/10.107/S0108767392000898“>导数格的枚举和对称重要性质，Act.Cryst.A48（1992），500-508。表1，对称性Pmmm。

%F From _Amiram Eldar_，2022年10月25日：（开始）：

%如果p>2，F与a（2^e）=3*e*（e+1）+1和a（p^e）=（e+1）*（e+2）/2相乘。

%F和{k=1..n}a（k）~（13/8）*n*log（n）^2+c1*n*log（n jes常数（A082633）。（完）

%p读取（“转换”）：nmax:=100：

%p L：=[1,4,0,1，seq（0，i=1..nmax）]：

%p莫比乌西（%）：

%p莫比乌西（%）：

%p莫比乌西（%）；#_R.J.Mathar，2017年9月25日

%tf[p，e]：=（e+1）*（e+2）/2；f[2，e_]：=3*e*（e+1）+1；；a[1]=1；a[n_]：=倍@@f@@FactorInteger[n]；阵列[a，100]（*_Amiram Eldar_，2022年10月25日*）

%o（PARI）t1=方向（p=2200,1/（1-X）^3）

%o t2=方向（p=2,2,1+4*X+X^2200）

%o t3=dirmul（t1，t2）

%o A145399（n）=t3[n]；\\该行由_Antti Karttune_于2018年9月27日添加

%o（PARI）a（n）={my（f=因子（n）

%Y参考A007425，A145511。

%Y参考A001620、A082633。

%K nonn，多个

%O 1,2号机组

%A _N.J.A.Sloane，2009年3月13日