（PARI）/*使用矩阵幂法生成：*/a（n）=局部（a=Mat（1），B）；对于（m=1，n+2，B=矩阵（m，m）；对于（i=1，m，对于（j=1，i，如果（j==i，B[i，j]=1，B[i，j]=（A^（fibonacci（i-1）））[i-1，j]）； ));A=B）；返回（（（A）[n+2,2]））
（PARI）/*使用部分和方法生成：*/a（n）=局部（a=向量（n+1），p）；A[1]＝1；对于（j=1，n，p=fibonacci（n+3）-fibonacci（n-j+3）-j；A=Vec（（Polrev（A）+x*O（x^p））/（1-x））；答[#A]
