(PARI){a(n)=如果(n<1,0,sumdiv(n,d,moebius(n/d)*(斐波那契(d-1)+斐波那契(d+1))^6)/n)}
(PARI){Lucas(n)=斐波那契(n-1)+斐波那奇(n+1)}
{a(n)=局部(F=exp(sum(m=1,n,Lucas(m)^7*x^m/m)+x*O(x^n));如果(n==1,1,polcoeff(F*prod(k=1,n-1,(1-Lucas(k)*x^k+(-1)^k*x^(2*k)+x*O(x^n))^a(k)),n)/Lucas(n)}