OFFSET
0,3
COMMENTS
The matrix M(n) is nonsingular only for n in {1, 5, 6, 7, 8, 9}.
EXAMPLE
a(4) = 29376:
[1, 2, 3, 4]
[2, 4, 6, 8]
[3, 6, 9, 21]
[4, 8, 21, 61]
MAPLE
rev:= proc(n) local L, i;
L:= convert(n, base, 10);
add(L[-i]*10^(i-1), i=1..nops(L))
end proc:
f:= proc(n) LinearAlgebra:-Permanent(Matrix(n, n, (i, j) -> rev(i*j))) end proc:
f(0):= 1:
map(f, [$0..20]); # Robert Israel, Feb 16 2026
MATHEMATICA
a[n_]:=Permanent[Table[IntegerReverse[i*j], {i, n}, {j, n}]]; Join[{1}, Array[a, 15]]
PROG
(PARI) a(n) = matpermanent(matrix(n, n, i, j, fromdigits(Vecrev(digits(i*j))))); \\ Michel Marcus, Feb 16 2026
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Stefano Spezia, Feb 16 2026
STATUS
approved
