%I#27 2021年5月20日10:55:16
%S 1,1,2,4,9,18,40,80168340698139628445688114294846072,
%电话9214418469636939273361479232295986059197201184269623685473,
%电话:47376634947539401895195763790391527581029001516205800
%N在变换T下左移,其中Ta是DCONV a。
%H Reinhard Zumkeller,n表,n=1..1000的a(n)</a>
%H N.J.A.Sloane,转换</a>
%F From _Benoit Cloitre,2004年8月29日:(开始)
%F a(n+1)=和{d|n}a(d)*a(n/d),a(1)=1。
%F a(素数(k)+1)=2*a(素素(k));
%F a(n)渐近于c*2^n,其中c=0.353030198…(结束)
%计算公式:A(x)=和{n>=1}A(n)*x^n=x*(1+Sum_{i>=1}和{j>=1}A(i)*A(j)*x*(i*j))_伊利亚·古特科夫斯基(Ilya Gutkovskiy),2019年5月1日[修改人:伊莉亚·古特科夫斯基(_Ilya Gutcovskiy),2018年5月9日]
%F a(n+1)=和{k=1..n}a(gcd(n,k))*a(n/gcd(n,k))/phi(n/gccd(n、k)),其中phi=A000010_理查德·奥尔勒顿(Richard L.Ollerton),2021年5月19日
%p带有(numtheory);EIGENbyDIRCONV:=proc(upto_n)局部n,a,j,i,s,m;a:=[1];对于i从1到upto_n做s:=0;m:=转换(除数(i),集合);n:=nops(m);对于从1到n的j,dos:=s+(a[m[j]]*a[m[(n-j)+1]]);od;a:=[操作(a),s];od;返回(a);结束;
%t dc[b_,c_]:=模[{p},p[n_]:=p[n]=和[b[d]*c[n/d],{d,如果[n<0,{},除数[n]]}];p] ;A[n_,k_]:=模[{f,b,t},b[1]=dc[f,f];对于[t=2,t<=k,t++,b[t]=dc[b[t-1],b[t-2]];f=函数[m,如果[m==1,1,b[k][m-1]];f[n]];a[n_]:=a[n,1];阵列[a,40](*_Jean-François Alcover_,2017年3月20日,A144324*之后)
%o(哈斯克尔)
%o导入数据。功能(打开)
%o a038044 n=a038044_列表!!(n-1)
%o a038044_list=1:f 1[1]其中
%o f x ys=y:f(x+1)(y:ys)其中
%o y=总和$zipWith((*)`on`a038044)divs$reverse divs
%o其中divs=a027750_row x
%o——_Reinhard Zumkeller_,2014年1月21日
%奇数项的Y位置由A003095给出。其他自进化序列:A000108、A007460-A007464、A025192、A061922、A062177。
%Y列k=1(A144324和A144823)。-_Alois P.Heinz,2012年11月4日
%Y参考A038040。
%Y参考A000010。
%K nonn,本征
%氧1,3
%基督教G.鲍尔_
|