%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日