login
A278857
a(n) = permanent M_n where M_n is the n X n matrix m(i,j) = (i-j)^2.
2
1, 0, 1, 8, 676, 49600, 10335908, 2658757248, 1214367336000, 730771063280640, 642638269862752320, 736176718456263406080, 1122592471007868379259136, 2168016139899273930219233280, 5288852927890824307509101287680, 15889369670472598370104100032512000
OFFSET
0,4
LINKS
MAPLE
with(LinearAlgebra):
a:= n-> `if`(n=0, 1, Permanent(Matrix(n, (i, j)-> (i-j)^2))):
seq(a(n), n=0..16); # Vaclav Kotesovec, Nov 30 2016, after Alois P. Heinz
MATHEMATICA
Flatten[{1, Table[Permanent[Table[(i-j)^2, {i, 1, n}, {j, 1, n}]], {n, 1, 15}]}]
PROG
(PARI) {a(n) = matpermanent(matrix(n, n, i, j, (i-j)^2))}
for(n=0, 20, print1(a(n), ", ")) \\ Vaclav Kotesovec, Aug 12 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Nov 29 2016
STATUS
approved