%I#46 2024年2月19日01:58:30

%S 1,1,1,2,3,7,43452453512592029160027819184831，

%电话：758397193749171922281611126403219004744354228963383975713263866432，

%电话：456995262861174715209949594648733172150425791690122432447239675853643

%N Propp的三次递归：a（0）=a（1）=a；对于n>3，a（n）=（1+a（n-1）*a（n-2）*a。

%D _James Propp，个人沟通。

%H Seiichi Manyama，<a href=“/A051786/b051786.txt”>n表，n=0..18的a（n）</a>

%F a（0）=a（1）=a；对于n>3，a（n+2）*a（n-2）=1+a（n+1）*a。

%F a（-n）=a（n+3）。

%F From _Vaclav Kotesovec_，2015年5月20日：（开始）

%F a（n）~c1^（（1+sqrt（13）-sqrt（2*sqrt）-2）/4）^n）*c2^*arccot（平方（（2*sqrt（13）-5）/3）），其中

%传真c1=0.0858378165313271469223136812741638183980800626360336156811045938771。。。

%传真c2=1.04799811587376782356890606697393352445131341037535628998042343。。。

%F c3=1.068106045469569610547194501969993896102077685059613621050203396954。。。

%传真c4=0.0530316436302789163635657674114415892838612646004303528421194603。。。

%F（结束）

%t循环表[{a[1]==a[2]==a[3]==a[4]==1，a[n]==（1+a[n-1]a[n-2]a[n-3]）/a[n-4]}，a[n]，{n，15}]（*哈维·P·达尔，2011年5月14日*）

%o（PARI）{a（n）=如果（n<0，n=3-n）；如果（n<4，1，（a（n-1）*a（n-2）*a*/

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

%o（哈斯克尔）

%o a051786 n=a051786_列表！！n个

%o a051786_list=1:1:1:

%o zipWith div（尾部$zipWith3（\u v w->1+u*v*w）

%o（删除2 a051786_list）（尾部a051786 _list和051786_ list）

%o 051786_列表

%o——Reinhard Zumkeller，2014年1月7日

%Y参考A005246、A072713。

%K nonn，很好，很容易

%0、5

%迈克尔·索莫斯，1999年12月9日

%E根据Andrew S.Plewe的建议，由N.J.A.Sloane编辑，2007年6月17日