OFFSET
0,4
COMMENTS
a(n) = A000166(n) for n <= 9.
MAPLE
a:= n-> `if`(n=0, 1, LinearAlgebra[Permanent](Matrix(n, (i, j)->
`if`(i<>j and (i-j mod n<=4 or j-i mod n<=4), 1, 0)))):
seq(a(n), n=0..15);
MATHEMATICA
a[n_] := If[n == 0, 1, Permanent[Table[If[i != j && (Mod[i - j, n] <= 4 || Mod[j - i, n] <= 4), 1, 0], {i, 1, n}, {j, 1, n}]]]; Table[an = a[n]; Print["a(", n, ") = ", an]; an, {n, 0, 15}] (* Jean-François Alcover, Jan 06 2016, adapted from Maple *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Jul 15 2015
STATUS
approved