%I#19 2020年12月29日09:04:13

%S 1,3,21159125710046818116741845605141469208743949193，

%电话：3339464105283430820022413244707232060357315568176324524499，

%电话：15121111593088012991212194668411179500360964780963472766997340831385476662980947182777882757416692

%N使用从不超过对角线x=y的步骤{（k，0），（0，k）|0<k<=4}，在从（0,0）到（N，N）的所有晶格路径上求和的节点总数。

%H Alois P.Heinz，n的表格，n=0..1000时的a（n）</a>

%F a（n）~c*d^n/sqrt（n），其中d=8.84734830841870961487278801886633962039798…是方程4+4*d-8*d^2-8*d*3+d^4=0和c=0.436590775682493454220000202799136319250347607927734138960545201547…-Vaclav Kotesovec_，2017年5月30日

%p b:=proc（x，y）选项记忆`如果`（y>x或y<0，0，

%p`如果`（x=0，[1$2]，加上（（p->p+[0，p[1]））(

%p b（x-j，y）+b（x，y-j）），j=1..4））

%p端：

%pa:=n->b（n$2）[2]：

%p序列（a（n），n=0..25）；

%tb[x_，y_]：=b[x，y]=如果[y>x|y<0，{0，0}，如果[x==0，{1，1}，和[Function[p，p+{0，p[[1]]}][b[x-j，y]+b[x、y-j]]，{j，1，4}]]；

%ta[n]：=b[n，n][2]]；

%t a/@范围[0，25]（*_Jean-François Alcover_，2020年12月29日，在_Alois P.Heinz_*之后）

%Y参考A175891。

%K nonn公司

%0、2

%A _Alois P.Heinz，2017年5月17日