%I#6 2021年8月20日00:24:04

%S 2,0,0,1,0,4,0,9,4,4,,0,32,0,8,49,0,56,0,36,16,4,0153,4,8,56,0，

%电话96,0207,8,4,16295,0,8,16187,0168,0,48120,12,0553,16,80,8,78,0，

%U 444,8323,16,4,0480,0,12216745,16144,0,60,24200,01016,0,8152,90,16216,0723472,4,0786,8,8289

%N a（1）=2；当n>1时，a（n）=A336853（n）+A347098（n）。

%｛A336853和A063524（1，0，0，0，…）的逐点和｝的C和及其狄利克雷逆。

%C第一个负项是a（720）=-6306。

%H<a href=“/index/Pri#gaps”>质数的索引条目</a>

%H<a href=“/index/Pri#prime_indices”>根据素因式分解中的索引计算出的序列的索引项</a>

%F a（1）=2，对于n>1，a（n）=-求和{d|n，1<d<n}A336853（d）*A347098（n/d）。

%F对于所有n>=1，a（A001248（n））=A001223（n）^2。

%o（PARI）

%o up_to=16384；

%o A336853（n）={my（f=因子（n））；对于（i=1，#f~，f[i，1]=下一素数（f[i、1]+1））；（因子（f）-n）；}；

%o Aux347098（n）=如果（1==n，n，A336853（n））；

%o方向逆校正（v）={my（u=向量（#v））；u[1]=（1/v[1]）；对于（n=2，#v，u[n]=（-u[1]*sumdiv（n，d，if（d<n，v[n/d]*u[d]，0））；（u）}；\\计算输入向量v中给定序列的Dirichlet逆。

%o v347098=方向反向校正（矢量（up_to，n，Aux347098（n）））；

%o A347098（n）=v347098【n】；

%o A347099（n）=如果（1==n，2，A336853（n）+A347098（n））；

%Y参见A001223、A001248、A003961、A063524、A336853、A347098。

%Y另请参阅A346250、A347097。

%K符号

%O 1,1号机组

%2021年8月19日，安蒂·卡图宁