%I#8 2019年10月25日09:58:43

%S 1,2,3,4,5,6,8,10,12,13,14,16,18,21,22,24,26,27,30,34,36,42,45,48,55，

%电话：56,58,60,66,68,69,72,76,78,80,81,84,89,90,92,93,94,96,99102105108，

%电话：110111116120126132135140144146150152153156159162

%N Zeckendorf-已知数：可被Zeckendor表示中的项数整除的数（A007895）。

%D Andrew Ray，《关于k-Zeckendorf Niven数的自然密度》，密苏里州立大学博士论文，2005年。

%H Robert Israel，n的表，n的a（n）=1..10000</a>

%H Helen G.Grundman，<a href=“https://www.fq.math.ca/Papers1/45-3/grundman.pdf“>连续Zeckendorf-Niven和懒惰Fibonacci-Niven数，《斐波纳契季刊》，第45卷，第3期（2007年），第272-276页。

%H Andrew Ray和Curtis Cooper，<a href=“http://cs.ucmo.edu/~cnc8851/articles/kzeckniven.pdf“>关于k-Zeckendorf Niven数的自然密度。

%由于A007895（12）=3，且3是12的除数，因此e12是顺序中的。

%p fib:=组合：-fibonacci：

%pφ：=1/2+平方英尺（5）/2：

%p fibapp:=n->phi^n/sqrt（5）：

%p invfib:=进程（x:：posint）

%p局部q，n；

%pq:=evalf（（ln（x+1/2）+ln（5）/2）/ln（φ））；

%p n：=楼层（q）；

%p如果fib（n）<=x，则

%p，而fib（n+1）<=x do

%p n：=n+1

%p端do

%p其他

%当fib（n）>x do时为p

%p n：=n-1

%p端do

%p end if；

%第n页

%p端程序：

%p zeck：=进程（x）局部n；

%p如果x=0，则为0

%p其他

%p n：=invfib（x）；

%p F[n]+zeck（x-fib（n））；

%功率因数

%p端程序：

%p过滤器：=n->n mod nops（zeck（n））=0：

%p选择（过滤器，[1..200]）；#_罗伯特·伊斯雷尔，2019年10月25日

%tz[n_]：=长度[DeleteCases[NestWhileList[#-Fibonacci[Floor[Log[Sqrt[5]*#+3/2]/Log[GoldenRatio]]&，n，#>1&]，0]]；aQ[n_]：=可除[n，z[n]]；选择[范围[1000]，aQ]（*在A007895*处的_Alonso del Arte_后面）

%Y参考A005349、A007895。

%K nonn公司

%O 1,2号机组

%A _Amiram Eldar，2019年10月7日