%I#29 2023年9月17日01:36:57

%S-1,0,0,0,1,2,3,3,5,6,7,9,10,12,14,15,18,20,22,25,27,30,33,35,39,42，

%电话：45,49,52,56,60,63,68,72,76,81,85,90,95,99105110115121126132，

%电话：138143150156162169175182189195203210217225232

%N热带版本的Somos-4序列A006720。

%C给定广义Somos-4序列的变量s（1）、s（2）、s n）是Z.-Michael Somos_中所有n的不可约多项式，2023年9月16日

%H G.C.Greubel，n的表格，n=1..5000的a（n）</a>

%H A.Fordy和A.Hone，<A href=“http://arxiv.org/abs/1207.6072“>来自簇映射的离散可积系统和泊松代数，arXiv:1207.6072[nlin.SI]，2012，参见示例3.6。

%H A.P.Fordy，<A href=“http://arxiv.org/abs/1403.8061“>周期性簇突变和相关可积映射，arXiv预印本arXiv:1403.8061[math-ph]，2014年。

%H<a href=“/index/Rec#order_10”>带常系数线性递归的索引条目，签名（2，-1,0,0,0-0,0,1，-2,1）。

%F From _Michael Somos，2012年12月27日：（开始）

%传真：x*（x^6-x^5+x^4-x^2+2*x-1）/（（1-x）^2*（1-x^8））。

%F a（2-n）=a（n）。（结束）

%F第二个差异具有周期8。-_迈克尔·索莫斯（Michael Somos），2012年12月27日

%F a（n）=A236294（n-5）=Z.-Michael Somos中所有n的最大值（a（n-1）+a（n-3），2*a（n-2））-a（n-4），2023年9月16日

%e G.f.=-x+x^5+x^6+2*x^7+3*x^8+3*x^9+5*x^10+6*x^11+7*x^12+。。。

%p A118825x:=进程（n）

%p系数（（1-2*x+x^2）/（x^4+1），x=0，n）；

%p端程序：

%p A056594：=程序（n）

%p系数（1/（x^2+1），x=0，n）；

%p端程序：

%p A220838：=程序（n）

%p-9/32-1/8*n+1/16*n^2+1/32*（-1）^n；

%p%+A118825x（n）/4-A056594（n+3）/8；

%p端程序：

%p序列（A220838（n），n=0..80）；#_R.J.Mathar，2013年1月30日

%t线性递归[{2，-1，0，0，0,0，0,1，-2，1}，{-1，0、0，0、1，2，3，3，5}，62]（*Jean-François Alcover_，2017年11月26日*）

%t a[n_]：=使用[{m=n-1}，楼层[m^2/16]-布尔[Mod[m，8]==0]]；（*迈克尔·索莫斯，2023年9月16日*）

%o（PARI）｛a（n）=如果（n<1，n=2-n）；极系数（x*（x^6-x^5+x^4-x^2+2*x-1）/（（1-x）^2*（1-x^8））+x*o（x^n），n）｝/*_Michael Somos，2012年12月27日*/

%o（PARI）{a（n）=n-；n^2\16-！（n%8）}；/*_迈克尔·索莫斯（Michael Somos），2023年9月16日*/

%o（岩浆）m:=25；R<x>：=PowerSeriesRing（整数（），m）；系数（R！（x*（x^6-x^5+x^4-x^2+2*x-1）/（（1-x）^2*（1-x^8）））；//_G.C.Greubel，2018年8月10日

%Y参考A006720，A236294。

%K符号，简单

%O 1,7型

%A _N.J.A.Sloane，2012年12月23日