%I#48 2022年9月8日08:46:08

%S 0,1,12202080117812038411154440199743601131233611957283461，

%电话：11084934960191793804840108621050274118793835590881，

%电话：106437544333680184160409410152010429793134197921180458407386358101102013289607062600

%N数字m，使3*m+1和8*m+1都是正方形。

%C当然，所有项都是三角数。

%C数字m，使k*m+1和8*m+1都是正方形：

%Ck=1:A006454；

%C k=3：该序列；

%Ck=4:A029549；

%Ck=5:0，3，231，4560，333336，6575751。。。

%C k=6:A200999；

%Ck=7:A157879。

%C数字m，使3*m+1和k*m+1都是正方形：

%Ck=1:A045899；

%Ck=2:A045502；

%Ck=4:A059989；

%C k=5:A159683；

%Ck=6:8*A029546；

%Ck=7:A160695；

%C k=8：这个序列。

%H Bruno Berselli，n的表格，n=1..500时的a（n）</a>

%H<a href=“/index/Rec#order_05”>具有常系数的线性重复出现的索引条目，签名（1,98，-98，-1,1）。

%传真：x^2*（1+20*x+x^2）/（（1-x）*（1-10*x+x2）*（1+10*x+2））。

%F a（n）=a（n-1）+98*a（n-2）-98*a（n-3）-a（n-4）+a（n-5）。

%四边形的F G.F:

%F a（4k+1）:40*x*（52+3*x）/（（1-x）*（1-9602*x+x^2））；

%Fα（4k+2）：（1+2178*x+21*x^2）/（（1-x）*（1-9602*x+x2））；

%（4k+3）：（21+2178*x+x^2）/（（1-x）*（1-9602*x+x^2））；

%Fα（4k+4）:40*（3+52*x）/（1-x）*（1-9602*x+x^2））。

%t线性递归[{1，98，-98，-1，1}，{0，1，21，120，2080}，20]（*或*）系数列表[x（1+20 x+x^2）/（1-x）（1-10 x+x2）（1+10 x+x^2）），{x，0，20}]，x]

%o（PARI）a=矢量（20）；a[1]=0；a[2]=1；a[3]=21；a[4]=120；a[5]=2080；对于（i=6，#a，a[i]=a[i-1]+98*a[i-2]-98*a[i3]-a[i-4]+a[i-5]）；一

%o（Maxima）a[1]：0$a[2]：1$a[3]:21$a[4]:120$a[5]：2080$a[n]：=a[n-1]+98*a[n-2]-98*a[n-3]-a[n-4]+a[n-5]$makelist（a[n，n，1，20）；

%o（岩浆）I:=[0,1,211202080]；[n le 5选择I[n]else Self（n-1）+98*Self，n-2）-98*Selve，n-3，n-4，+Selve（n-5）：[1..20]]中的n；

%Y参考A000217。

%Y参见A006454、A029546、A0295402、A045502、A045%、A0899、A059989、A157879、A159683、A160695、A200999。

%K nonn，简单

%氧1,3

%2014年7月15日，A _布鲁诺·贝塞利

%E 2016年3月3日，Bbruno Berselli_将偏移量从0更改为1并采用公式