%I#8 2017年12月4日12:09:53
%S 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,21,
%电话:26,27,28,29,30,31,32,33,34,35,36,37,34,38,39,13,40,41,42,43,44,45,46,
%U 47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,54,71,72,73,67,74,75,76,52,77,78,79,80,81,82,83,84,85,86,73
%N过滤器序列组合A003557(N)和A023900(N)。
%C对于所有i,j:a(i)=a(j)=>A295887(i)=A295887(j)。
%H Antti Karttunen,n的表,n=1..65537的a(n)</a>
%F序列a(n)=(1/2)*(2+((A003557(n)+A295876(n))^2)-A003557(m)-3*A295875(n)的限制生长序列变换。
%o(PARI)
%o分配(2^30);
%o up_to=65537;
%o rgs_transform(invec)={my(om=Map(),outvec=vector(length(invec)),u=1);对于(i=1,长度(invesc),如果(mapisdefined(om,invec[i]),my(pp=mapget(om,invec[i];
%o写入to_b文件(start_offset,vec,bfilename)={对于(n=1,长度(vec),写入(bfilename,(n+start_offset)-1,“”,vec[n]);}
%o A003557(n)={my(f=系数(n));对于(i=1,#f~,f[i,2]=最大值(0,f[i,2]-1));系数回退(f);};
%o A023900(n)=总和(n,d,d*moebius(d));\\此功能来自2014年9月9日的_Charles R Greathouse IV_
%o v295876=rgs_transform(向量(up_to,n,A023900(n)))
%o A295876(n)=v295876[n];
%o提交的6(n)=(1/2)*(2+((A003557(n)+A295876(n;
%o写入to_b文件(1,rgs_transform(向量(up_to,n,Anotsubmitted6(n)),“b295886.txt”);
%Y参见A003557、A023900、A295876、A295887。
%K nonn公司
%O 1,2号机组
%A _Antti Karttunen_,2017年12月3日
|