%2018年11月13日03:17:21

%S 1,1,2,3,3,4,4,7,7,6,6,12,7,8,12,16,9,15,10,18,16,12,12,12,17,14,22，

%电话：24,15,30,16,34,24,18,24,48,19,20,28,48,21,40,22,36,45,24,24,27,32,37，

%U 36,42,27,54,36,40,30,30,96,31,32,60,78,42,60,34,54

%N a（N）=psi_c（N），其中Product_{k>1}1/（1-1/k^s）^A007897（k）=和{k>0}psi_c（k）/k^s。

%D F.V.Weinstein，《斐波那契分区》，预印本，1995年。

%H F.V.Weinstein，<a href=“https://arxiv.org/abs/math/0307150“>斐波纳契分区注释，arXiv:math/0307150[math.NT]，2003-2018。

%e G.f.=x+x ^2+2*x ^3+3*x ^4+3*x^5+4*x ^6+4*x^7+7*x ^8+7*x^9+。。。

%t dircon[v_，w_]：=模块[{lv=长度[v]，lw=长度[w]，fv，fw}，fv[n]：=如果[n<=lv，v[[n]]，0]；fw[n_]：=如果[n<=lw，w[[n]]，0]；表[DirichletConvolve[fv[n]，fw[n]、n、m]、{m、Min[lv、lw]}]；

%t a[n_]：=模[{a，v，w，m}，如果[n<1，0，v=表[Boole[k==1]，{k，n}]；对于[k=2，k<=n，k++，m=长度[整数位数[n，k]]-1；A=乘积[{p，e}=pe；如果[p==2，如果[e<3，e，2^（e-2）+2]，1+p^（e-1）（p-1）/2]，{pe，因子整数[k]}]；A=（1-x）^-A+x O[x]^m//正常；w=表[0，{n}]；对于[i=0，i<=m，i++，w[[k^i]]=系数[A，x，i]]；v=dircon[v，w]]；v[[n]]]；

%t阵列[a，68]（*Jean-François Alcover_，2018年11月12日，来自PARI*）

%o（PARI）{a（n）=我的（a，v，w，m，p，e）；如果（n<1，0，v=向量（n，k，k==1）；对于1）/2））；a=（1-x）^-a+x*o（x^m）；w=向量（n）；对于（i=0，m，w[k^i]=波尔科夫（a，i）；v=dirmul（v，w））；v[n]）}；/*_迈克尔·索莫斯（Michael Somos），2014年5月26日*/

%Y参考A007896、A007897。

%K nonn公司

%O 1,3

%费利克斯·温斯坦（wain（AT）ana.unibe.ch）

%E米歇尔·马库斯的新定义，2014年5月12日

%E定义由N.J.A.Sloane编辑，2014年5月26日