%I#9 2019年9月16日06:38:57
%S 1,4,8,6,12,14,16,9,18,20,24,26,28,22,39,15,36,35,40,38,57,34,48,49,
%电话51,44,46,33,60,77,64,32,75,56,81,68,76,58100,55,84111,88,62125,70,
%U 96,91,98,95134,72108,82141,80140,92120156124,94121,52152145
%N Shotgun(或Schrotschuss)数:递归极限B(k)=T[k](B(k-1)),其中B(1)=(1,2,3,4,5,…),T[k'是对所有正整数i的k(2i-1)和k(2i)项进行置换的变换。
%C序列是无素数的。
%H K.Strassburger,猎枪编号图</a>
%e B(1)=1,2,3,4,5,6,7,8,9,10,11,12,13,14,。。。
%e B(2)=1,4,3,2,5,8,7,6,9,12,11,10,13,16,。。。
%e B(3)=1,4,8,2,5,3,7,6,10,12,11,9,13,16,。。。
%e B(4)=1,4,8,6,5,3,7,2,10,12,11,14,13,16,。。。
%t最大值=66;b[1,j_]:=j;b[k_,j_]:=b[k,j]=b[k-1,j];Do[b[k,2j*k-k]=b[k-1,2j*k];b[k,2j*k]=b[k-1,2j*k-k],{k,2,max},{j,1,max}];a[n]:=b[最大值,n];表[a[n],{n,1,max}](*Jean-François Alcover_,2012年10月11日*)
%o(SageMath)
%o定义分配(s,k):
%o如果不是k.divides:返回0
%o返回(-1)^(s//k)*k
%o定义A064494(n):
%o s=n
%o表示范围(n,1,-1)中的k:
%o s-=分配(s,k)
%o返回s
%o打印([A064494(n)代表n in(1..66)])#_Peter Luschny_,2019年9月16日
%Y参见A064728、A064590、A064627。
%很好,不
%O 1,2号机组
%A Klaus Strassburger(斯特拉斯(AT)ddfi.uni-duesseldorf.de),2001年10月16日
|