OFFSET
0,3
PROG
(PARI) a(n) = matpermanent(matrix(n, n, i, j, sigma(gcd(i, j))));
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Dec 18 2022
STATUS
approved