%I#9 2013年3月2日20:31:53

%S 1,5,3,9,33,11,13,7,21,25,20,17,23,19,49,39,35,47,27,29,15,43,37,97，

%电话：55,41,51,57,85,73193111,59,61,31,89,53329145192221237,87105，

%电话：81101113481175149，69,93,77,99,79,71,67,65,45,75195157141189109,8310322912135284923345297137309281379155199159

%N A114183中第N次加倍运行的开始。

%H N.J.A.Sloane，N的表，N的A（N）=1..7328</a>

%电话：A114181

%p M：=10000；M2:=1000；

%ps1:={1}；v0：=[1]；v1：=[1]；v2:=[]；vi：=阵列（1..M2）；

%pt1:=1；r1:=1；vi[1]：=1；

%p表示n从2到M do

%p t2:=楼层（sqrt（t1））；

%p如果s1中的t2，则

%p v0：=[op（v0），2*t1]；s1:={op（s1），2*t1}；r1：=r1+1；t1:=2*t1；

%p如果t1<=M2，则vi[t1]：=n；fi；

%p其他

%p v0：=[操作（v0），t2]；s1：＝{op（s1），t2}；v1：=[op（v1），t2]；v2：=[op（v2），r1]；r1:=1；t1:=t2；

%p如果t1<=M2，则vi[t1]：=n；fi；

%p fi；

%日期：

%编号A114183：

%p[seq（v0[i]，i=1..nops（v0））]；

%第A221715页：

%p[seq（v1[i]，i=1..nops（v1））]；

%第A221716页：

%p[seq（v2[i]，i=1..nops（v2））]；

%电话：A189419

%p[seq（vi[i]，i=1..M2）]；

%Y参见A114183、A189419、A221716、A213220。

%K nonn公司

%O 1,2号机组

%A _N.J.A.Sloane，2013年1月27日