编辑
经核准的
P: =proc(n)选项记忆;如果n<=3,则为1
(1+进程名(n-1)*进程名(n-2)*进程名称(n-3))/进程名(n-4)
fi;结束:seq(P(n),n=0..14)#保罗·拉瓦2018年9月27日之后罗伯特·伊斯雷尔在A018896号
Propp的三次递归:a(0) = a(1) = a(2) = a(3) = 1;对于 n个>三, a(n)=(1+a(n-1)*a(n-2)*a(n-3))/a(n-4)。
a(0)=a(1)=a;对于n > 3,a(n+2)*a(n-2)=1+a(n+1)*a。
a(-n)=a(n + 3).
一个(-n个) = 一(n个 + 三). 一(0) = a(1) = a(2) = a(3) = 1. ; 对于 n个 > 三, a(n) + 2) * a(n) - 2) =1+a(n + 1) * a(n) * a(n) - 1).
a(-n)=a(n+3)。
提出
(PARI)a=矢量(15);a[1]=a[2]=a[3]=1;a[4]=2;对于(n=5,#a,a[n]= (1+a[n-1]*a[n-2]*a[n-3])/a[n-4]);concat(1,a)\\阿尔图·阿尔坎2018年9月27日