%I#20 2022年3月29日04:10:12
%S 1,3,4,5,7,8,10,11,13,15,16,18,19,21,22,23,25,26,28,30,31,33,35,
%电话:36、37、39、39、41、42、44、46、47、49、50、51、53、54、56、57、59、59、60、61、62、64、66,
%U 68、70、71、73、73、75、76、77、79、80、82、84、85、87、88、89、90、91、91、93、94、96、98
%N a(N)=b(N+2)+b(N),其中a(0)=1,其中b(N)=A005229(N)对于N>2。
%C应用于斐波那契序列(A000045)的类似定义导致卢卡斯序列(A0000 32)。
%H G.C.Greubel,n的表格,n=0..5000时的a(n)</a>
%对于n>=1,F a(n)=A005229(n+2)+A005228(n)。
%p b:=proc(n)选项记忆;如果n<=2,则1其他b(b(n-2))+b(n-b(n-2)):fi:结束:seq(b(n),n=1..75):a[0]:=1:对于从1到70的n,执行a[n]:=b(n;
%t M[n_]:=M[n]=如果[n<3,1-布尔[n==0],M[M[n-2]]+M[n-M[n-2]]];
%t L[n_]:=L[n]=如果[n==1,1,M[n-1]+M[n+1]];
%t表[L[n],{n,100}](*由_G.C.Greubel_修改,2022年3月28日*)
%o(圣人)
%o@CachedFunction
%o定义b(n):#A005229
%o如果(n<3):返回1
%o else:返回b(b(n-2))+b(n-b(n-3))
%o定义A116591(n):返回b(n+2)+b(n)-bool(n==0)
%o[A116591(n)代表n in(0..100)]#_G.C.Greubel_,2022年3月28日
%Y参考A000032、A005185、A005229。
%K nonn公司
%0、2
%A _Roger L.Bagula,2006年3月27日
%E编辑:N.J.A.Sloane_,2006年4月15日