cp's OEIS Frontend

A278857 a(n) = permanent M_n where M_n is the n X n matrix m(i,j) = (i-j)^2.

%I A278857 #12 Aug 14 2021 01:59:27
%S A278857 1,0,1,8,676,49600,10335908,2658757248,1214367336000,730771063280640,
%T A278857 642638269862752320,736176718456263406080,1122592471007868379259136,
%U A278857 2168016139899273930219233280,5288852927890824307509101287680,15889369670472598370104100032512000
%N A278857 a(n) = permanent M_n where M_n is the n X n matrix m(i,j) = (i-j)^2.
%H A278857 Vaclav Kotesovec, <a href="/A278857/b278857.txt">Table of n, a(n) for n = 0..37</a>
%p A278857 with(LinearAlgebra):
%p A278857 a:= n-> `if`(n=0, 1, Permanent(Matrix(n, (i, j)-> (i-j)^2))):
%p A278857 seq(a(n), n=0..16);  # _Vaclav Kotesovec_, Nov 30 2016, after _Alois P. Heinz_
%t A278857 Flatten[{1, Table[Permanent[Table[(i-j)^2, {i, 1, n}, {j, 1, n}]], {n, 1, 15}]}]
%o A278857 (PARI) {a(n) = matpermanent(matrix(n, n, i, j, (i-j)^2))}
%o A278857 for(n=0, 20, print1(a(n), ", ")) \\ _Vaclav Kotesovec_, Aug 12 2021
%Y A278857 Cf. A085750, A085807, A204249, A278845, A278847, A278858.
%K A278857 nonn
%O A278857 0,4
%A A278857 _Vaclav Kotesovec_, Nov 29 2016