%I#17 2022年9月8日08:46:24
%S 1,3,5,16,11,43,17,88,48,95,31320,41145157486,59554,67696243,
%电话:265,8322042183474581062109191127266844749352350444157,
%电话:5655774780179302119381998697211145905162538823252624622939
%N a(1)=1;a(n)=Sum_{d|n,d<n}素(n/d)*a(d)。
%F G.F.A(x)满足:A(x”)=x+Sum_{k>=2}素数(k)*A(x^k)。
%t a[n_]:=如果[n==1,n,和[If[d<n,素数[n/d]a[d],0],{d,除数[n]}]];表[a[n],{n,1,55}]
%t nmax=55;A[_]=0;Do[A[x_]=x+和[Prime[k]A[x^k],{k,2,nmax}]+O[x]^(nmax+1)//正常,nmax+1];系数列表[A[x],x]//静止
%o(PARI)seq(n)={my(v=向量(n));v[1]=1;对于(n=2,#v,v[n]=sumdiv(n,d,v[d]*prime(n/d)));v}\\_Andrew Howroyd_,2019年9月5日
%o(岩浆)溶胶:=[1];对于[2.60]中的n,执行追加(~sol,&+[NthPrime(Floor(n/d))*sol[d]:Set中的d(Divisors(n));结束;溶胶;//_Marius A.Burtea,2019年9月5日
%Y参考A000040、A007445、A008966(a(n)的奇偶校验)、A034696。
%K nonn公司
%O 1,2号机组
%A_Ilya Gutkovskiy_,2019年9月5日
|