%I#6 2020年3月29日17:14:31

%编号：1,2,4126416442302545889243104281190015070181761943620532，

%电话：269562858232108330283527839293763447678503866195269254，

%电话：74578889849353499512096334100326102297142894144039145768147664152817163125183002

%N在Zeckendorf及其所有除数（A300837和A333618）的对偶Zeckenderf表示中可被1的总数整除的数。

%e 126是一个术语，因为A300837（126）=21和A333618（126）=7都是126的除数。

%t zeckDigSum[n_]：=长度[DeleteCases[NestWhileList[#-Fibonacci[Floor[Log[Sqrt[5]*#+3/2]/Log[GoldenRatio]]&，n，#>1&]，0]]；

%t zeckDivDigSum[n_]：=除数总和[n，zeckDigSum[#]&]；

%t fibTerms[n_]：=模块[{k=天花板[Log[GoldenRatio，n*Sqrt[5]]，t=n，fr={}}，而[k>1，如果[t>=Fibonacci[k]，则附加到[fr，1]；t=t-斐波纳契[k]，附录[fr，0]]；k-]；fr]；

%t dualZeckSum[n_]：=模块[{v=fibTerms[n]}，nv=长度[v]；i＝1；当[i<=nv-2时，如果[v[i]]==1&v[i+1]]==0&&v[[i+2]]==0，v[i]=0；v[[i+1]]=1；v[[i+2]]=1；如果[i>2，i-=3]]；i++]；i=位置[v，_？（#>0&）]；如果[i=={}，0，总计[v[[i[[1，1]]-1]]]]];

%t dualZeckDivDigSum[n_]：=除数总和[n，dualZekSum[#]&]；

%t选择[Range[10^4]，Divisible[#，zeckDivDigSum[#]]&&Divisible[#，dualZeckDivSum[#]]&]

%A333619和A333620的Y交叉口。

%Y参考A007895、A112310、A300837、A330711、A333617和A333618。

%K nonn，基础

%O 1,2号机组

%A _Amiram Eldar，2020年3月29日