%I M2435#30 2023年3月9日08:04:53
%S 0,0,1,3,5,8,11,14,18,22,26,30,34,38,43,48,53,58,63,68,73,78,83,89,95,
%电话101107113119125131137143149155161167173179185191198,
%电话:20521221922623240247254261268275289296
%N基2中低差异序列的数量。
%D N.J.A.Sloane和Simon Plouffe,《整数序列百科全书》,学术出版社,1995年(包括该序列)。
%H Harald Niederreiter,<a href=“http://dx.doi.org/10.1016/0022-314X(88)90025-X“>低分辨率和低色散序列,《数论》30(1988),第1期,第51-70页。
%p N:=进程(b,N)
%p选项记忆;
%p局部d;
%p加(b^d*numtheory[mobius](n/d),d=numtheori[divisors](n));
%p%/n;
%p端程序:
%p M:=进程(b,n)
%p局部h;
%p如果n=0,则
%p 0;
%p其他
%p加(N(b,h),h=1..N);
%p end if;
%p端程序:
%p nMax:=进程(b,s)
%p局部n;
%从0 do到n的p
%p如果M(b,n)>s,则
%p返回n-1;
%p end if;
%p端do:
%p结束过程:
%p A005356:=进程
%p局部n,b;
%磅:=2;
%p n:=nMax(b,s);
%p n*(s-M(b,n))+加((h-1)*n(b,h),h=1..n);
%p端程序:
%p序列(A005356(n),n=1..40);#_R.J.Mathar,2016年6月9日
%tNp[b_,n_]:=Np[b,n]=和[b^d*MoebiusMu[n/d],{d,除数[n]}]/n;
%tM[b_,n_]:=如果[n==0,0,和[Np[b,h],{h,1,n}]];
%t nMax[b_,s_]:=模[{n},对于[n=0,True,n++,如果[M[b,n]>s,返回[n-1]]];
%t a[s_]:=模块[{n,b},b=2;n=最大[b,s];n*(s-M[b,n])+和[(h-1)*Np[b,h],{h,1,n}]];
%t表[a[n],{n,1,60}](*Jean-François Alcover_,2023年3月9日,在R.J.Mathar_*之后)
%Y参见A005357(基数3)、A005377(基数4)、A00.5358(基数5)。
%K nonn公司
%O 1,4型
%A _N.J.A.Sloane,西蒙·普劳夫_
%E来自Sean A.Irvine_的更多条款,2016年5月27日
|