%I#17 2023年8月7日03:55:34

%S 0,1,1,3,3,6,12,21,38,71123234427791147745222984918659，

%电话：3541267410128644245959471166904186173823833465426452030，

%电话：12455921240764584659176690258683175029533

%N幸运数字l在2的幂之间，2^N<l<=2^（N+1）。

%e a（0）=0，因为在1（2^0）和2（2^1）之间没有幸运数字。

%e a（3）=3，因为在8（2^3）和16（2^4）之间有3个幸运数字（9，13，15）。

%o（SageMath）

%o def幸运（n）：

%o L=列表（范围（1，n+1，2））；j=1

%o而L[j]<=长度（L）-1：

%o如果（i+1）%L[j]，则L=[L[i]表示范围（长度（L））中的i=0]

%o j+=1

%o回流（L）

%o A000959=幸运（1048576）

%o定义幸运范围（a，b）：

%o幸运=[]

%A000959中的o表示l：

%o如果l>=b：

%o幸运回归

%o如果l>=a:lucky.append（l）

%o[len（lucky_range（（2^n）+1,2^（n+1）））用于范围（19）中的n）]

%Y参考A000959、A036378。

%K nonn，更多

%0、4

%A _霍克·Löffler，2019年7月28日

%E a（19）-a（30）摘自Giovanni Resta，2020年5月10日

%E a（31）-a（32），来自_Kevin P.Thompson_，2021年11月22日