%I#40 2021年12月12日22:11:09

%S 0,1,2,5,8,9,10,21,32,33,34,37,40,41,42,85128129130136137，

%电话：138149160161162165168169170341512513514517520521522，

%电话：533544545546955525535545976406416426456489650661

%N a（N）是k的最小值，A020986（k）=N。

%C Brillhart和Morton导出了k的最大值的ω函数。这个序列似乎是由一个类似的α函数给出的。

%H Michael Day，n的表，a（n）表示n=1..10000</a>

%H J.Brillhart和P.Morton，<a href=“http://www.maa.org/programs/maa-awards/writing-awards/a-case-study-in-mathematical-research-the-golay-rudin-shapiro-sequence“>数学研究案例研究：Golay-Rudin-Shapiro序列，美国数学月刊，103（1996）854-869。

%H Kevin Ryde，<a href=“http://user42.tuxfamily.org/alternate/index.html“>交替折纸曲线的迭代，见索引GRScumulFirstN。

%F a（2*n-1）-a（2*n-2）=（2^（2*k+1）+1）/3和a（2*n）-a。

%F重申上述公式，a（n+1）-a（n）=A007583（A050605（n-1））=A276391，术语重复_约翰·基思，2021年3月4日

%o（PARI）

%α（n）={

%o如果（n<2，返回（最大值（0，n-1））；

%o本地（nm1=n-1，

%o mi=m=天花板（nm1/2），

%o r=楼层（log（m）/log（2）），

%o i，fi，α=0，a）；

%o对于步骤（i=1,2*r+1，2，

%o mi/=2；

%o fi=（1+2^i）\3；

%oα+=fi*地板（0.5+mi）；

%o）；

%oα*=2；

%o如果（nm1%2，\\调整为偶数n

%o a=系数（2*m）[1,2]-1；

%oα-=（1+2^（1+2*a））\3；

%o）；

%o返回（α）；

%o}（o）

%o（J）

%o注：。向量上的J函数

%o注：。注意大参数的舍入错误

%o注：。确定到~1e8

%o字母=：3:0

%o n=.<：年

%如果.+/ntlo=，则为o。n>0 do。

%o n=。ntlo#n

%o m=.>.-：n个

%o r=.<。2^.米

%o f=.<。3%~2+2^2*>：i.>./>：第页

%o z=。0

%o mi=。米

%o代表i.i.f做。

%o z=。z+（i{f）*<0.05+mi=.mi%2

%o结束。

%o nzer=。(+/ @: (0=>./\)@:|.)“1@：#：m

%o ntlo#^：_1 z-（2|n）*<.-：nzer{f

%o其他。

%o ntlo公司

%o结束。

%o）

%o注：。例如alphav 1 3 5 100 2 8 33

%Y参见A020985、A020991、A020986。

%K nonn公司

%O 1,3

%2012年5月22日，迈克尔日

%E小修编辑：N.J.A.Sloane，2012年6月6日