%I#8 2017年6月16日22:33:29

%S 1,4,1,16,5,1,64,23,6,1256101,30,7,11024467138,37,8,140962165，

%电话：694175,44,9,116384100553526925212,51,10,1655364670918012，

%电话：49771156249，58,11,1262144216995921402706764281387286，65,12,1

%N反对偶读取的数组：T（m，N）=长度为N的m元单词数，其循环相邻元素相差不超过3。

%C所有行都是具有常数系数的线性递归。请参阅PARI脚本以获取生成函数。

%H Andrew Howroyd，n表，n=4..1278的a（n）</a>

%e表格开始（m>=4，n>=0）：

%e 1 4 16 64 256 1024 4096 16384 65536。。。

%e 1 5 23 101 467 2165 10055 46709 216995。。。

%e 1 6 30 138 694 3526 18012 92140 471566。。。

%e 1 7 37 175 925 4977 27067 147777 808165。。。

%e 1 8 44 212 1156 6428 36338 206942 1183164。。。

%e 1 9 51 249 1387 7879 45663 267367 1575395。。。

%e 1 10 58 286 1618 9330 54994 328058 1973026。。。

%e 1 11 65 323 1849 10781 64325 388749 2371457。。。

%e 1 12 72 360 2080 12232 73656 449440 2770016。。。

%t差值=3；m0=差异+1；mmax=13；

%t传输Gf[m_，u_，t_，v_，z_]：=数组[u，m]。LinearSolve[IdentityMatrix[m]-z*数组[t，{m，m}]，数组[v，m]]

%t行Gf[d_，m_，z_]：=1+z*总和[TransferGf[m，Boole[#==k]&，Boole[绝对值[#1-#2]<=d]&，布尔值[Abs[#-k]<=d]&，z]，{k，1，m}]；

%t行[m_]：=行[m]=系数列表[RowGf[diff，m，x]+O[x]^mmax，x]；

%t t[m_/；m>=m0，n_/；n>=0]：=行[m][[n+1]]；

%t表[t[m-n，n]，{m，m0，mmax}，{n，m-m0，0，-1}]//Flatten（*Jean-François Alcover_，2017年6月16日，改编自PARI*）

%o（PARI）

%o TransferGf（m，u，t，v，z）=向量（m，i，u（i））*matsolve（matid（m）-z*矩阵（m，m，i，j，t（i，j）），向量v（m，i，v（i）））；

%o行Gf（d，m，z）=1+z*总和（k=1，m，转移Gf（m，i->如果（i==k，1,0），（i，j）->abs（i-j）<=d，j->如果（abs（j-k）<=d,1,0），z））；

%o表示（m=4，12，打印（第Gf行（3，m，x））；

%o表示（m=4，12，v=Vec（第Gf行（3，m，x）+o（x^9））；对于（n=1，长度（v），打印1（v[n]，“，”）；打印（）；）；

%Y行5-32为A124999、A125316-A125342。

%Y参见A285267、A285280、A276562。

%K nonn，表

%氧4.2

%A A Andrew Howroyd_，2017年4月15日