%I #19 Apr 02 2019 14:34:42
%S 0,2,3,256,38880,1289945088,42855402240000000,
%T 605828739547255327948800000000,
%U 13263549731442762279026688000000000000000000000000000,1334793240853871268746431553848403294648071618560000000000000000000000000000000000000000000
%N a(n) is the Narumi-Katayama index of the Lucas cube Lambda(n).
%C The Lucas cube Lambda(n) can be defined as the graph whose vertices are the binary strings of length n without either two consecutive 1's or a 1 in the first and in the last position, and in which two vertices are adjacent when their Hamming distance is exactly 1.
%C The Narumi-Katayama index of a connected graph is the product of the degrees of the vertices of the graph.
%H I. Gutman and M. Ghorbani, <a href="https://doi.org/10.1016/j.aml.2011.12.018">Some properties of the Narumi-Katayama index</a>, Applied Mathematics Letters, Vol. 25, No. 10 (2012), 1435-1438.
%H S. Klavžar, M. Mollard and M. Petkovšek, <a href="https://doi.org/10.1016/j.disc.2011.03.019">The degree sequence of Fibonacci and Lucas cubes</a>, Discrete Mathematics, Vol. 311, No. 14 (2011), 1310-1322.
%e a(2) = 2 because the Lucas cube Lambda(2) is the path-tree P_3 having 2 vertices of degree 1 and 1 vertex of degree 2; consequently, the Narumi-Katayama index is 1*1*2 = 2.
%p G := (1+(1-y)*x+x^2*y^2+(1-y)*x^3*y-(1-y)^2*x^4*y)/((1-x*y)*(1-x^2*y)-x^3*y):
%p g := expand(series(G, x=0, 40)): T := (n, k) -> coeff(coeff(g, x, n), y, k):
%p a := n -> mul(k^T(n, k), k=0..n): lprint(seq(a(n), n=1..10));
%Y Cf. A307157, A307181, A307307.
%K nonn
%O 1,2
%A _Emeric Deutsch_, Mar 28 2019