%I#14 2018年9月1日22:27:42
%S 1,2,9,6,25,9,49,20243,25121,27169,49225,70289243361,75441,
%电话:121529,901875169364514784122596125210892257291369,
%U 361152125016814411849363607552920931572031875260150728093645302549324984134816753721961190792442251089
%N其Dirichlet卷积与自身产生平方的序列的分子,A000290。
%H Antti Karttunen,n的表格,n=1..65537的a(n)</a>
%F a(n)=F(n)的分子,其中F(1)=1,F(n)=(1/2)*((n^2)-Sum_{d|n,d>1,d<n}F(d)*F(n/d))对于n>1。
%F a(n)=n*A318512(n)*A299149(n)/A299150(n)。
%o(PARI)
%o up_to=65537;
%o DirSqrt(v)={my(n=#v,u=vector(n));u[1]=1;对于(n=2,n,u[n]=(v[n]/v[1]-sumdiv(n,d,if(d>1&&d<n,u]*u[n/d],0))/2);u};
%o v318649_aux=DirSqrt(向量(up_to,n,(n*n)));
%o A318649(n)=分子(v318649_aux[n]);
%Y参见A000290,A318512(分母)。
%Y另请参阅A046643、A299149、A318511、A318651、A318654(给出了偶数项的位置)、A318655(2元估值)。
%K nonn公司
%O 1,2号机组
%安蒂·卡通纳,2018年8月31日
|