%I M0616#22 2023年12月21日10:23:23

%S 1,2,3,5,6,7,9,11,12,13,15,16,17,19,21,22,23,25,26,27,29,30,31,33,35，

%电话：36、37、39、40、41、43、45、46、47、49、50、51、53、54、55、57、59、60、61、63、64、65、67，

%U 69,70,71,73,74,75,77,79,80,81,83,84,85,87,88,89,91,93,94,95,97,98,99101单位

%N光谱阵列W的第一行2-洗牌（平方码2）。

%D N.J.A.Sloane和Simon Plouffe，《整数序列百科全书》，学术出版社，1995年（包括该序列）。

%H Aviezri S.Fraenkel和Clark Kimberling，<a href=“http://dx.doi.org/10.1016/0012-365X（94）90259-3“>Generalized Wythoff arrays，shuffle and interspersions，《离散数学》126（1-3）（1994）137-149。【摘自R.J.Mathar_，2009年8月17日】

%p位数：=200:WythSpec:=proc（n，x）floor（n*x）；结束：A001951:=进程（n）WythSpec（n，sqrt（2））；结束：A001952:=进程（n）A001951（n）+2*n；end:Wsqrt2:=proc（i，j）选项记住；如果j=1，则为A001951（A00195l（i））；elif j＝2，则为A001952（A001951（i））；elif类型（j，‘add’），然后是A001951（procname（i，j-1））；其他A001952（程序名（i，j-2））；fi；end:A070071:=proc（n）选项记住；局部a；如果n=1，则为1；如果Wsqrt2（k，1）=a，则返回（a）；elif Wsqrt2（k，1）>a然后中断；fi；od：对于从1开始的k，如果Wsqrt2（k，2）=a，则返回（a）；elif Wsqrt2（k，2）>a然后中断；fi；od:od:fi；结束：序列（A007071（n），n=1..100）；#_R.J.Mathar，2009年8月17日

%t WythSpec[n_，x_]：=地板[n*x]；

%t A001951[n_]：=WythSpec[n，Sqrt[2]；

%t A001952[n]：=A001951[n]+2n；

%t WSqrt2[i_，j_]：=WSqrt2[i，j]=其中[j==1，A001951[A001951[i]]，j==2，A001 952[A00195 1[i]，奇数Q[j]，A00195 1[WSqrt2]，j-1]]，真，A001 1952[WSqrt2[i，j-2]]；

%t A007071[n_]：=A007071[n]=模块[{a，k}，如果[n==1，1，对于[a=A007071-[n-1]+1，真，a++，对于[k=1，真，k++，如果[WSqrt2[k，1]==a，返回[a]，如果[WSqrt2[k，1]>a，中断[]]]；对于[k=1，True，k++，If[WSqrt2[k，2]==a，Return[a]，If[WSqrt2[k，2]>a，Break[]]]]]；

%t表[A007071[n]，{n，1，72}]（*_Jean-François Alcover_，2023年12月20日，在R.J.Mathar_*之后）

%K nonn，简单

%O 1,2号机组

%A _N.J.A.Sloane，_Mira Bernstein_

%E更多条款摘自R.J.Mathar_，2009年8月17日